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ABSTRACT 

FIRST MEASUREMENT OF THE RATIO OF BRANCHING FRACTIONS 
B{K ^ A>-I7^)/S(A6 A+TT-) AT CDF II 
Shin-Shan Yu 
Nigel Lockyer 

We present the first measurement of the ratio of branching fractions i3(A5 — >■ 
fi^Vfj^) /B{Ai) A+TT^) based on 171.5 pb^-*^ of pp colhsions at ^/s = 1.96 TeV 
taken with the CDF-II detector. In addition, we present measurements of B{B — > 
D*+^-Vf,)/B(B" D*+7i-) and B(b" D+ii-V^)/B(^ D+n-), which 
serve as control samples to understand the data and Monte Carlo used for the A;, 
analysis. We find the relative branching fractions of the control samples to be: 

^^^(B^ ^D*+ -) " ^^'^ ^ ^'^ ^ ^^^^^^ ^ ^^^^ ^ ^'^ (UBR), 

and 

B(^ ^ D+fl-T7^^ 
B(^ D+TT-) 



9.8 ± 1.0 {stat) ± 0.6 {syst) ± 0.8 {BR) ± 0.9 {UBR), 



which are consistent with the ratios obtained by the Particle Data Group at the 
0.7 and 1.1 a level, respectively. Finally, we obtain the relative Af, branching 
fraction to be: 

^i^l ^ ^if^^^ = 20-0 ± 3.0 {stat) ± 1.2 {syst) {BR) ± 0.5 {UBR). 

The uncertainties of the three relative branching fractions are from statistics, 
CDF internal systematics, external measured branching ratios and unmeasured 
branching ratios, respectively. We present a method to derive B{Ah —>■ A+tt^) 
using previous CDF measurements and obtain 

B{A, ^ A+TT-) = (o.41 ± 0.19 {stat syst) ^Jj'^g (Pt)) %, 

where the last uncertainty is due to the measured A;, Pt spectrum. Combining 
B{Ah — i> A+TT^) with our result, we determine the exclusive semileptonic branch- 
ing fraction for the A;,; 

B{A, ^ A+/i-l7^) = (8.1 ± 1.2 {stat) ^ {syst) ± 4.3 {B{A, ^ A+Tr"))) %. 
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Chapter 1 
Introduction 



In this dissertation, we measure the properties of the lowest-mass beauty baryon, 
Aft. Baryons are the bound states of three quarks. Protons and neutrons, con- 
stituents of atomic nuclei, are the most common baryons. Other types of baryons 
can be produced and studied in the high-energy collider environment. Three- 
body dynamics makes baryons composed of low mass quarks difficult to study. 
On the other hand, baryons with one heavy quark simplify the theoretical treat- 
ment of baryon structure, since the heavy quark can be treated the same way as 
the nucleus in the atom. The A;, is composed of u, d, and b quarks, where the 
b quark is much heavier than the other two. Although, it is accessible, little is 
known about A;,. In 1991, UAl [1] reconstructed 9 ± 1 A;, — > J/^'A candidates. 
In 1996, ALEPH and DELPHI reconstructed the decay A;, A+n~ and found 
only 3-4 candidates [2,3]. ALEPH measured a A^ mass of 5614 ± 21 MeV/c^ 
while DELPHI measured 5668 ± 18 MeV/c^, about 2 a higher. Subsequently, 
CDF-I observed 20 A5 J/^A events [4], confirmed the existence of A5 unam- 
biguously and made a more precise measurement of A^ mass, 5621 ± 5 MeV/c^. 
A recent CDF-II measurement by Korn [5] yields 5619. 7± 1.7 MeV/c^, which will 
significantly improve the current world average, 5624±9 MeV/c^, and resolve the 
discrepancy of ALEPH and DELPHI. 

Several experiments have also measured the product of a fragmentation frac- 
tion and a branching ratio, such as: /Aji3(Aft — > J/^'A) [4], and f\^B{Ab — >■ 
A+/i^X) [6,7]. However, branching ratios derived from measurements rely on 
the knowledge of the A;, fragmentation fraction (/aJ, which is defined as the 
probability for a b quark to hadronize into A5. Assuming that the A;, dominates 
the production of the beauty baryons, i.e. = /baryom and applying the world 
average /baryon compiled by the Particle Data Group (PDG) [8], we obtain the 
branching ratios: 

B{Ab J/*A) = (4.7 ± 2.8) x 10-^ 

and 

B{Ab A+ii-X) = (9.2 ± 2.1) X 10"^ 

The uncertainties on the above branching ratios arc large, about 60% and 20%. 
Several other decays, have been searched for, but either only 1 candidate was 
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A, A, 




Figure 1.1: Feynman diagram of — > A+/i i/^ (left) and — > A+tt (right) 
decays. 



observed (eg: A5 A^ai(1260)~) or no candidates were found (eg: A;, pK~) 
and an upper limit was set. In addition to the mass and branching ratios, the 
Ab lifetime is an important physics quantity. However, the world average lifetime 
ratio, T(Ab)/T(S°), is 0.797±0.052, in disagreement with the range of theoretical 
predictions: between 0.9 and 1.0 [9]. The properties listed above are all that is 
known about A;, and its decays, which motivates us to measure the branching 
ratios of Af,. 

Currently, the Fermilab Tevatron is the only facility which produces a large 
sample of A;,. The Heavy Quark Effective Theory (HQET) has been used to 
predict the mass, lifetime, and decay rates of the A5. Studying A5 at the Tevatron 
through various measurements allows us to test HQET in different aspects. We 
present a measurement of the relative branching fractions of A;, — > A+/x~I7^ to 
Afe — > A+7r~. Figure 1.1 shows that these two decays have very similar Feynman 
diagrams: a b quark decays into a c quark via a virtual W boson exchange, 
and W decay into a muon and an anti- neutrino or an u and a d quark. The 
advantage of measuring the ratio of branching fractions is that several systematic 
uncertainties cancel, such as those from the trigger and reconstruction efficiencies. 
To understand the A;, measurement, we perform a similar analysis on the better 
understood decays: we measure B{B D*+ II y^/B{B — > D*+7r ) and 
^ D+ii-V^)/B(B^ ^ D+TT-). 

The analysis strategy is as follows: The number of signal events observed in 
the data (N) is the product of the number of produced B hadrons {Nb), the 
branching ratio (i3), the detector acceptance and reconstruction efficiency (e) ob- 
tained from a Monte Carlo (MC) program, i.e. N — Nb ■ B ■ e. Therefore, the A;, 
(S°) relative branching fraction is the yield ratio divided by the efficiency ratio 
since the numbers of B hadrons cancel. For the semileptonic mode, several back- 
grounds exhibit a similar signature to the real signal. We estimate the amount of 
these backgrounds (iVbg) and subtract them from the observed yield in the data. 
The formula for extracting the relative branching ratio is then expressed as: 



B, 



semi ^ 



inclusive semi 



had 



had 



A^bg X Chad 
-) X 
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In this dissertation, we measure the following relative branching fractions: 

-> D*+^-V^) /B(^ D*+TT-), where D*+ L)°7r+, D° ^ K-7r+, 
B(B° D+fi-V^)/B(B^ D+TT-), where D+ K-n+n+, 
B{Ab ^ A+/x-I7^)/i3(Ab ^ A+TT-), where A+ ^ pK-7r+. 

In Chapter 2, a brief overview of the Standard Model and HQET is presented. 
A general description of the Tevatron, the CDF-II detector and trigger is found 
in Chapter 3. Chapter 4 details our data sample and event selection. Chapter 5 
focuses on how the signal yields are extracted. The MC simulations for the 
acceptance and efficiency are described and compared with data in Chapter 6. 
Chapter 7 involves the estimate of the backgrounds present in the semileptonic 
signal. In Chapter 8, the systematic uncertainties are first discussed, then the 
results of the relative branching fractions are summarized. Through out the whole 
dissertation, the analyses of the A^ and the control modes are presented in 
parallel. The charge conjugates of the A;, and decays are also included in the 
reconstruction. 
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Chapter 2 

Theoretical Background 



This chapter first gives a brief overview of the fundamental particles and inter- 
actions in Section 2.1. Then a general idea of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix and the Heavy Quark Effective Theory (HQET) are discussed in 
Sections 2.2-2.3. 

2.1 Fundamental Particles and Interactions 

The "Standard Model" [10, 11] is the accepted theory that describes particles 
with no internal structure and their interactions with matter. In the Standard 
Model, the basic constituents of matter are six flavors of spin-1/2 quarks: the 
down-type quarks {d, s, h) and the up- type quarks (u, c, t), and six kinds of spin- 
1/2 leptons: the charged leptons, e, /x and r, and the neutrinos ^fi, and z/^-. 
The magnitude of the electron's electric charge is denoted as "e" . The down-type 
quarks carry electric charge — |e and the up- type quarks carry electric charge +|e. 
The charged leptons have electric charge -le while the neutrinos have zero electric 
charge. The masses of the quarks and leptons exhibit a hierarchy. The u, d and 
s quarks are much lighter (< 200 MeV/c^) than the c, h and t quarks (~ 1100, 
4500 and 175000 MeV/c^). This hierarchy is not understood. The three charged 
leptons also have progressively increasing masses: ~ 0.51(e), 106(/x), 1777 (r) 
MeV/c^. For the neutrinos, currently, only upper limits exist for their masses. 
Separate neutrino types can undergo transitions into one another if at most one 
type of neutrino has zero mass. The current best estimates require three-flavor 
mixing to explain the full range of results from the solar neutrino [12, 13], and 
the atmospheric neutrino experiments [14-16]. 

The Standard Model describes the following three types of interaction among 
quarks and leptons: electromagnetic, weak and strong interactions. The gravi- 
tational interaction is not described by the Standard Model. The gravitational 
force dominates in the large mass scale, such as a galaxy, but has little influence 
on the scale of quarks and leptons. Therefore, it is usually ignored in the funda- 
mental particle interactions. The quarks and leptons interact via the exchange 
of the gauge bosons. The Lagrangian of each interaction is invariant under a 
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transformation that corresponds to a symmetry group. The Standard Model is a 
theory based on the symmetry group SU{3) x SU{2) x U{1). Both SU{n) and 
U{n) groups are Lie groups, i.e. any element in the group can be represented by 
m fundamental elements or generators [17]: 

m 

E^eMH&iFi), (2.1) 

i 

where Fi is the i^^ generator and 9i is the "rotation" angle corresponding to each 
generator. Elements of the SU {n) groups are represented by n x n imitary matri- 
ces, U^U = 1, with det U = +1 and have n^ — 1 generators. The theory introduces 
r? — \ gauge bosons, analogous to the rotation angle in Equation 2.1. They form 
a scalar product with the r? — \ generators and make the Lagrangian invariant. 
The t^(l) group is a one dimensional unitary group with single generator, where 
the elements are specified by a continuous parameter, and expressed as e*^. 

The U{X) group describes the electromagnetic interaction among quarks and 
the charged leptons, via the exchange of a massless spin-1 photon. The electro- 
magnetic interaction binds the electrons and atomic nuclei together and forms 
atoms. The SU{2) group describes the weak interaction experienced by all the 
fundamental particles, where the gauge bosons are massive spin-1 and 
The masses of and Z° are about 80 and 91 GeV/c^. A well known weak 
interaction process is the neutron /3-decay: n{udd) —>■ p{uud)ei7g. The right- and 
left-handed fundamental particles transform differently under SU{2). The right- 
handed quarks and leptons do not couple to and arc singlets under SU{2). 
While the left-handed quarks and leptons are doublets under SU (2) and classified 
into three generations: 




The weak interaction allows transitions between quarks of different flavors. The 
transitions within the same generation are more favored than those across the gen- 
erations and the coupling strength is given by the CKM matrix (see Section 2.2). 
The coupling strength is the same for all leptons. 

The SU{?>) group describes the strong color interaction among quarks, me- 
diated via the exchange of eight massless spin-1 gluons. The quarks carry three 
possible "chromoelectric charges": red, green and blue (i?, G, -B), which are anal- 
ogous to the electric charge in the electromagnetic interaction. The eight gluons 
are associated with the color combinations: 

RB, RG, BR, BG, GR, GB, {RR - GG) /^/2 and {RR + GG- 2BB) / ^6 

The strong interaction binds the quarks together to form a colorless state, qq or 
qqq {RR or RGB). The qq bound state is referred to as "meson" and the qqq 



5 



bound state is referred to as "baryon" . For example, a bd bound state is a 

meson and a udb bound state is a beauty baryon, A;,. Both mesons and baryons 
are called "hadrons". Just as the residual electric field outside of the neutral 
atoms causes them to combine into molecules, the residual color field outside of 
the protons and neutrons forms nuclei. 

Each fundamental particle has an associated antiparticle, i.e. of which the 
electric charge, color charge and fiavor are reversed, but the mass and the spin 
are the same. In addition, the Standard Model introduces a neutral spin-0 Higgs 
boson, if", to accommodate the masses of the gauge bosons, quarks and leptons. 
The Higgs boson has not been discovered, yet. The search for the Higgs boson 
remains an important goal of several running and future high energy experiments. 



2.2 CKM Matrix 

The strong interaction conserves the fiavor of quarks, and only transitions of the 
same fiavor quark will take place (eg: charmness conserved decay, D*^ — > D+vr*'), 
while the flavor-changing decays are allowed in the electroweak interaction (eg: 
beauty to charm decay, — > D*+7r~). The Cabibbo - Kobayashi - Maskawa 
(CKM) matrix [18, 19] in Equation 2.3 describes the coupling in the weak in- 
teraction between different flavors of quarks. For instance, Vch describes the 
electroweak coupling strength of the b quarks to the c quarks. The CKM ma- 
trix represents a unitary transformation from the flavor (mass) eigenstates to the 
weak interaction eigenstates. 



Vc 



CKM 



Vud 




Vub 


Vcd 




Vc, 


Vu 


Vts 


Vt, 



(2.3) 



A standard CKM matrix parametrization proposed by Chau,et al. [20-23], which 
is similar to Kobayashi - Maskawa's original parametrization [19], has four free 
parameters: three mixing angles between any two generations, ^12, 6*23, 6*13 and a 
phase, 5i3; 



V 



CKM 



/ C12C13 

"•§12023 — Ci2<S23Sl3e 
V •S12S23 - Ci2C23Sl3e 



iSi3 



( 



C12 



ii5i3 



-S12C23 - Ci2S23Sl3e 
\ S12S23 - Ci2C23Sl3e*'^" 



S12C13 
C12C23 ~ 'Si2>S23Si3e 
~Ci2S23 ~ Si2C2ZSi^e 



i5i3 



S12 

C12C23 - Sl2S23Sl3e 
~Ci2S23 ~ Si2C2ZSi^e 



ii5i3 



si3e-*^« \ 
S23C13 (2.4) 
C23C13 / 



sise 
S23 
C23 



(2.5) 



with c. 



cos 9 ij and sij 



sinOij where i, j labels the generations. The matrix 
elements with simpler forms appearing in the first row and third column, have 
been measured directly in decay processes. C13 is known to be very close to unity: 
|ci3 — 1| ~ 10~^, and this gives an approximation in Equation 2.5. In the Standard 
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Model, the complex phase, ^13, is the origin of the Charge-Parity (CP) violation 
in the weak interaction. We refer the reader to the Tevatron Run II B Physics 
Workshop Report [9] for a more detailed description of the weak CP violation 
mechanisms. 

Using the world average experimental results of the weak decays as the input, 
and assuming that only three generations exist, with the unitarity, a 90% confi- 
dence limit can be placed on the amplitude of the matrix elements (eg: 



/ 0.9739 ~ 0.9751 0.221 ~ 0.227 0.0029 ~ 0.0045 \ 

0.221 ~ 0.227 0.9730 ~ 0.9744 0.039 ~ 0.044 
^ 0.0048 ~ 0.014 0.037 ~ 0.043 0.9990 ~ 0.9992 y 



(2.6) 



As seen in Equation 2.6, quark transitions within the same generation are fa- 
vored over the transitions across generations. The latter are called "Cabibbo 
suppressed" decays. The goal of several analyses at CDF II, together with the 
experiments BELLE [24], BABAR [25] and KTeV [26], etc. , is to make precise mea- 
surements of many matrix elements. For example, a measurement of the — 
oscillation frequency can infer the CKM matrix element Vtd- While V^d can be 
obtained by comparing the nuclear /3-decay n peV^ to muon decay and V^s 
can be measured from the decay: K 7rei7g. By making measurements of all 
the CKM elements, we can determine whether the CKM matrix is unitary. If 
the matrix is not unitary, this would be a signature for additional physics not 
currently described by the Standard Model. 



2.3 Heavy Quark Effective Theory 

This dissertation presents a measurement of the relative decay rates. The 
transition amplitude {M.) that describes the decay rate of a S hadron into some 
final state /, can be derived by drawing all the possible Feynman diagrams at 
the quark level and summing up all the contributions. The underlying weak 
interaction is simple but the strong interaction that binds the quarks into hadrons 
introduces complications. When the quarks or gluons travel over a distance of 
1/Aqcd or longer, the coupling constant of the strong interaction (a^) diverges, 
so perturbation theory breaks down and the nonperturbative effect takes over. 
For the energy scale of our concern, Aqcd is around 200 MeV/c^. 

One theoretical tool, the Operator Product Expansion (OPE) [27], separates 
the perturbative from nonperturbative physics and M. becomes: 



M = 



VckmEC,-(/|0,-[S) 



i + e>( 



ml . 



(2.7) 



where j indicates the contribution from the j*^ Feynman diagram, is the 
Fermi coupling constant, Vckm is the CKM matrix element in Equation 2.3. 
The Wilson coefficients [27], Cj, act as effective coupling constants and contain 



7 



the physics at short distance. The Wilson coefficients can be calculated using 
perturbation theory and are model independent. The {f\Oj\B) is usually referred 
to as the hadronic matrix element, where Oj is a local operator. The hadronic 
matrix elements contain the long distance physics and can only be evaluated 
using nonperturbative methods. Contributions from the higher order operators 
are suppressed by a power of m^/M^, where and Mw are the masses of 
the b quark and the W boson. The Heavy Quark Effective Theory (HQET) 
significantly simplifies the form of the hadronic matrix element. This section 
gives a review of the HQET and shows how the decay rates may be derived 
using the HQET with other theoretical assumptions. 

The HQET stems from the Standard Model and describes the hadrons con- 
taining a 6 or c quark. The concept of "heavy" is relative. In the HQET, the 
masses of the "heavy" c, b and t quarks are much larger than QCD energy scale, 
while the masses of the "light" u, d and s quarks are much smaller than Aqcd- 
In the limit of mc,6,t ^ ^qcd, a new type of symmetry, "spin-flavor heavy quark 
symmetry" arises. The momentum transfer between the heavy quark and the 
light quarks in the hadron system is of the order of Aqcd- Ot cquivalently speak- 
ing, the typical size of a hadron system is of the order of Ag^^,. The change 
in the heavy quark velocity is then ~ Kqcd/^Qi which vanishes when mq is 
infinitely large. The velocity of the heavy quark is, therefore, almost unaffected 
by the strong interaction, i.e. the quark-quark interaction terms disappear in the 
Lagrangian. The only strong interaction of a static heavy quark is with gluons 
via choromoelectric charge. This quark-gluon interaction is spin-independent. 
Consequently, the light quark system knows nothing about the spin, mass and 
flavor of the "nucleus" , i.e. a B hadron at rest is identical to a charm hadron at 
rest regardless of their spin orientations. 

The "heavy quark symmetry" implies that we can relate properties of the 
beauty hadrons to those of the charm hadrons. For example, Aglietti [28] derived 
a formula to estimate the Af, mass: Ma, - \{Md + SM^,*) = Ma^ - \ {Mb + 3Mb* ) , 
which gives Ma^ ~ 5630 MeV/c^, in good agreement with the world average, 
5624±9 MeV/c^. An analogy can be found in atomic systems, where the isotopes 
with different nuclei have nearly the same chemical properties. When performing 
a calculation of the B or charm hadron mass, decay rate or lifetime, we could 
start from the limit of m^b,* S> ^qcd- Then the correction terms arc added 
in expansion of the power of l/mq, where ttiq is the mass of the heavy quark. 
The l/mq corrections take into account flnite mass effects and are different for 
quarks of different masses. A more complete description of HQET may be found 
in Manohar [29], Godfrey [30] and Isgur [31]. 

The focus of this analysis, examining A^ to Ac baryon decay, is best suited 
to treatment using HQET since both the initial- and the flnal- state hadrons 
contain a heavy quark. In addition, the light quark system in a A^ baryon is 
in a spin-0 state; the sub-leading corrections have a simpler form than those for 
the mesons [32]. This analysis concerns the relative branching fractions of A;, — >■ 
A+/i~I7^ to A;, — > A+7r~, where the leading order Feynman diagrams are shown in 
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Figure 1.1. Using the Soft CoUinear Effective Tfieory (SCET), Leibovich, Ligeti, 

Stewart, and Wise [33] relate the decay rate of A5 ^^^r" to A5 A+/i~z7^. The 
SCET assumes that in the A5 — A+vr" decay, the pion mass can be neglected and 
rrib — rric = Et^^ Aqcd- Therefore, the ud quark pair acquires large momentum, 
remains close together and acts as a color dipole (singlet). Within the A^, the 
h quark and the diquark [ud) form a color dipole since the A^ is colorless; the 
same hold for the c quark and the diquark within the Ac- The color dipole- 
dipole interaction is weaker than the color monopole-monopole interaction, which 
means the pion interacts weakly with the rest of the system. In the end, the 
Ab — > A+7r~ decay factorizes into two subprocesses: the hadronization of the 
Ac and, completely decoupled, the hadronization of the tt. The hadronic matrix 
element is then expressed as (Ac|c7'^(l— 7^)6|A6)(7r|(i7^(l— 7^)m|0). The first term 
is common to both hadronic and semileptonic decays and can be inferred from 

the differential decay rate of the semileptonic mode, ^ at maximal 

recoil, where w is the scalar product of the A;, and A^ four- velocities, v and v': 

w = v-v' (2.8) 
= K, + ^L-?')/(2^A6^Ac)- (2-9) 

Here is the four-momentum transfer in the decay. Maximal recoil refers to the 
kinematic configuration when the charged lepton and the neutrino momenta are 
parallel, or equivalently, when is at its minimum, mf, which is approximately 
zero. At this configuration, semileptonic decay. A;, — > A+/i~I7^, has the same 
kinematics as the hadronic decay, A^ — > A+tt", since the (f of the hadronic decay 

is always m^, which is neglected in the SCET, and = 0. The '''^^^^^^^ ^ '^"^ can 
be constructed from six form factors, which are functions of w. The form factors 
can be considered as the Fourier transformation of the weak charge distribution 
and describe the interaction between the h and the c quarks. The second term of 
the hadronic matrix element is the pion decay constant, /^r. The value of /^r is 
131 MeV and was extracted from the decay width of 7r~ — > /I'Vi^. 
In a more exact form, the A5 A+tt" decay rate is 

mi ta \ dw I 

where Ci and C2 are the first two terms of the Wilson coefficients, the higher 
order terms are suppressed, and 

ta = mAe/rriAf (2.11) 

The Wmax corresponds to = mf ~ ~ 0. In the limit of uiq 3> Aqcd, the 
six form factors that describe the semileptonic differential decay rate are reduced 
to one universal function, the Isgur-Wise function {({w)) [31,34], and 

dTjAb A+/X-I7 ) GlmijKbP 3 r-2 — Tk, 2 . 2^2. n 

— — — -2-^ r\\/w^ — l[Qw + owr]^ — 4rA — 8r\w \Q [w). 

dw 2477'^ 

(2.12) 
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Note that Equation 2.12 proposes an alternative way to measure \ Vcb\ using A{, — > 
A+/i~z7^ decay. 

Although HQET reduces the form factors to the Isgur-Wise function, it can 
not predict the functional form of C(^)- One functional form easy for calculation 
was suggested by Isgur and Wise [31]: 

((w) = e-^'^"'-^^ (2.13) 

where the slope has to be calculated using other theoretical assumptions. 
Assuming the number of colors (Nc) in the baryon is infinitely large, Jenk- 
ins, et al. [35] derived = 1.3. Using the QCD sum rules, Huang, et al. [36] 
calculated = 1.35 ± 0.12. A recent DELPHI measurement [37] gives = 
2.03 ± 0A6{stat)^^'J^^{syst), consistent with the numbers from Huang and Jenk- 
ins, et al. Combining Equations 2.10-2.13 and the slope value from Jenkins, 
Leibovich, et al. predict B{Ah A+tt") = 0.45% and B{Ab A+p-Jy^) = 6.6%. 
However, the correction to the large Nc limit is of order 1/Nc, which is 30% in 
the case of baryons {Nc = 3). The uncertainty from the QCD sum rule is about 
10%. 

To summarize, measuring the relative branching fractions of A^ — > A+/i~l7^ to 
Ah — > A+7r~ allows us to obtain a ratio free from several experimental systematic 
uncertainties. With the external input of i3(Ab — * A+vr"), we could infer the 
i3(Af, A+/i~z7^) and vice versa. Finally, the absolute branching ratios of A;, 
A+7r~ and A;, — > A+/x~I7^ increase our knowledge of the A^ baryon and test the 
validity of Equations 2.10 and 2.12 which are derived from the HQET and the 
SCET. 
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Chapter 3 

The CDF-II Detector and Trigger 



The Fermilab Tevatron is currently the highest energy accelerator in the world. 
Protons and anti-protons (pp) are accelerated in its ~ 6 km (4 mile) circumfer- 
ence to be brought into collision with a center of mass energy of approximately 
1.96 TcV. The collisions take place at the center of two detectors: the Collider 
Detector at Fermilab (CDF-II) and DO-II. These two detectors are about 120° 
away from each other on the ring as indicated in Figure 3.1. The "luminos- 
ity" is a measure of collision rate normalized by the collision cross section, in 
unit of 1/sec-cm^. The common dimension for the time integrated luminosity 
is "barn~^", which is 10^'^/cm^. During 1992-1995, the predecessors of CDF-II 
and DO-II, CDF-I and DO-I had collected data with a time-integrated luminos- 
ity of ~110 pb~^ (inverse pico-barn) and pubhshed more than 100 papers. This 
analysis uses data collected by the CDF-II experiment. 

Both the accelerator and the collider detectors underwent major upgrades 
between 1997 and 2001. The main goals of these upgrades were to increase the 
luminosity of the accelerator, and to collect data samples with an integrated lumi- 
nosity of 2 fb"^ (inverse femto-barn) or more. The upgraded Tevatron accelerates 
36 bunches of p and p, whereas the previous accelerator operated with only 6x6. 
Consequently, the time between collisions (or beam crossings) has decreased from 
3.5 fis to 396 ns for the current collider. The new collider configuration required 
extensive detector upgrades at CDF-II to accommodate the shorter bunch spac- 
ings. In Section 3.1, we give an overview of how the proton and anti- proton beams 
are accelerated to their final center of mass energy of 1.96 TeV, and collided. We 
then describe in Sections 3.2 -3.5 the components of the CDF-II detector, and 
trigger, which are used to measure the properties of the particles produced in the 
pp collisions. 

3.1 pp acceleration and collisions 

In order to create the world's most energetic particle beams, Fermilab uses a series 
of accelerators. The diagram in Figure 3.1 shows the paths taken by protons and 
anti-protons from initial acceleration to collision in the Tevatron. The first stage 
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Figure 3.1: Layout of the Fermilab accelerator complex. The proton (solid arrow) 
is accelerated at the Cockcroft- Walton, Linac, Booster, Main Injector and finally 
at the Tcvatron. The anti-proton (dashed arrow) from the anti-proton source is 
first accelerated at the Main Injector and then at the Tevatron. 
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of acceleration is in the Cockcroft- Walton pre-accelerator [38] , where H~ ions are 
created from the ionization of the hydrogen gas and accelerated to a kinetic energy 
of 750 keV. The ions enter a hnear accelerator (Linac) [39], approximately 
500 feet long, where they are accelerated to 400 MeV. The acceleration in the 
Linac is done by a series of "kicks" from Radio Frequency (RF) cavities. The 
oscillating electric field of the RF cavities groups the ions into bunches. Before 
entering the next stage, a carbon foil removes the electrons from the H~ ions at 
injection, leaving only the protons. The 400 MeV protons are then injected into 
the Booster, a 74.5 m-diameter circular synchrotron. The protons travel around 
the Booster about 20,000 times to a final energy of 8 GeV. 

Protons are then extracted from the Booster into the Main Injector [40], 
where the protons are accelerated from 8 GeV to 150 GeV before the injection 
into the Tevatron. The Main Injector also produces 120 GeV protons, where the 
protons collide with a nickel target, and produce a wide spectrum of secondary 
particles, including anti-protons. In the collisions, about 20 anti-protons are 
produced per one million protons. The anti-protons are collected, focused, and 
then stored in the Accumulator ring. Once a sufficient number of anti-protons 
are produced, they are sent to the Main Injector and accelerated to 150 GeV. 
Finally, both the protons and anti-protons are injected into the Tevatron. The 
Tevatron, the last stage of Fermilab's accelerator chain, receives 150 GeV protons 
and anti-protons from the Main Injector and accelerates them to 980 GeV. The 
protons and anti-protons travel around the Tevatron in opposite directions. The 
beams are brought to collision at the center of the two detectors, CDF-II and 
DO-II. 

We use the term "luminosity" to quantify the beam particle density and the 
crossing rate. The luminosity in units of cm~^s~^ can be expressed as: 

where / is the revolution frequency, Nb is the number of bunches, -A^p/p are the 
number of protons/anti-protons per bunch, and ap/p are the RMS beam sizes 
at the interaction point. F is a form factor which corrects for the bunch shape 
and depends on the ratio of ai, the bunch length, to (3*, the beta function, at 
the interaction point. The beta function is a measure of the beam width, and is 
proportional to the beam's x and y extent in phase space. Figure 3.2 shows the 
peak luminosities for the stores used in this analysis. The collision products are 
recorded in the CDF-II and DO-II detectors. 



3.2 The CDF-II Detector 

The Collider Detector at Fermilab (CDF) is a general purpose, azimuthally and 
forward-backward symmetric apparatus, designed to study pp collisions at the 
Tevatron. Figure 3.3 shows the detector and the different sub-systems in a solid 
cutaway view. 
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Collider Run II Peak Luminosity 
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Figure 3.2: Peak luminosities for stores collided between April 2001 and February 
2004. This analysis uses the data collected from February 2002 to September 
2003. 



Standard Coordinates in CDF-II 

Because of its barrel-like detector shape, CDF-II uses a cylindrical coordinate 
system (r, 0, z) with the origin at the center of the detector. The z axis is along 
the direction of the proton beam. The r indicates the radial distance from the 
origin and is the azimuthal angle. The r-0 plane is called the transverse plane, 
as it is perpendicular to the beam line. The polar angle, 6, is the angle relative 
to the z axis. An alternative way of expressing 9, pseudorapidity (rj), is defined 
as: 

T] = -\ntan{9/2). (3.2) 

The coverage of each CDF-II detector sub-system will be described using combi- 
nations of 1], r, (f) and z. 

Overview 

The CDF-II detector consists of five main detector systems: tracking, particle 
identification (for e, K, p and vr), calorimetry, muon identification and luminosity 
measurement. 

The innermost system of the detector is the integrated tracking system: a 
silicon microstrip system and an open-cell wire drift chamber, the Central Outer 



14 



5^ 




MinipluQ' 




Calorimeter \ 



Figure 3.3: The CDF-II detector with quadrant cut to expose the different sub- 
detectors. 

Tracker (COT) that surrounds the sihcon detector. The tracking system is de- 
signed to measure the momentum and the trajectory of charged particles. Re- 
constructed particle trajectories are referred to as "tracks". Multiple-track re- 
construction allows us to identify a vertex where either the pp interaction took 
place (primary vertex) or the decay of a long-lived particle took place (secondary 
or displaced vertex). 

The silicon microstrip detector consists of three sub-detectors in a barrel ge- 
ometry that extends from the radius of r= 1.35 cm to r= 28 cm and covers the 
track reconstruction in the range of |?7| < 2. Closest to the beam pipe is a single- 
sided, radiation tolerant silicon strip detector, Layer 00 (LOO), with sensors at 
r=1.35 cm and r=1.62 cm. LOO is followed by five concentric layers of double- 
sided silicon sensors (SVX-II) from r=2.45 cm to 10.6 cm. The outermost silicon 
detector is the Intermediate Silicon Layers (ISL), from r=20 cm to r=28 cm. LOO 
only provides r — (j) measurements, while the SVX-II and ISL provide both r-(f) 
and z measurements. 

Surrounding the silicon detector is the COT, which covers the radius from 
40 cm to 137 cm and \7]\ < 1. Figure 3.4 and Figure 3.5 give a r-z view of the 
CDF-II tracker and the silicon tracking system, respectively. 

Immediately outside the COT is the Time of Flight system (TOF), which 
consists of 216 scintillator bars, roughly 300 cm in length and with a cross-section 
of 4 X 4 cm^. The bars are arranged into a barrel around the COT outer cylinder. 
The TOF system is designed for the particle identification of charged particles 
with momentum below 2 GeV/c. Both the tracking system and the TOF are 
contained in a superconducting solenoid, 1.5 m in radius, 4.8 m in length, that 
generates a 1.4 Tesla magnetic field parallel to the beam axis. The solenoid is 
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Figure 3.4: A diagram of the CDF-II tracker layout showing the different subde- 
tector systems. 
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Figure 3.5: Coverage of the different sihcon subdetector systems projected into 
the r-z plane. The r and z axes have different scales. 
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Figure 3.6: Locations of the CDF-II Muon System in the rj-ct) view. 

surrounded by the electromagnetic and hadronic calorimeters, which measure the 
energy of particles that shower when interacting with matter. The coverage of 
the calorimeters is \ri\ < 3. The electromagnetic calorimeter is a lead/scintillator 
sampling device and measures the energy of the electrons and photons. The 
hadronic calorimeter is an iron/scintillator device and measures the energy of the 
hadrons, e.g.: pions. 

The calorimeters are surrounded by the muon detector system. Muons inter- 
act with matter primarily through ionization. As a result, if a muon is created 
in the collision and has enough momentum, it will pass through the tracking 
system, TOF, the solenoid and the calorimeters with minimal interaction with 
the detector material. Muon detectors are, therefore, placed radially outside the 
calorimeters. The CDF-11 detector has four muon systems covering different 77 
and (f) regions: the Central Muon Detector (CMU), the Central Muon Upgrade 
Detector (CMP), the Central Muon Extension Detector (CMX), and the Interme- 
diate Muon Detector (IMU). Figure 3.6 shows the coverage of each muon detector 
in the 'q-cf) view. 

At the extreme forward region of the CDF-II detector, 3.75 < l^^l < 4.75, 
two modules of Cherenkov Luminosity Counters (CLC) are placed pointing to 
the center of the interaction region to record the number of pp interactions. The 
number of particles recorded in the CLC modules is combined with an accep- 
tance of the CLC {A) and the inelastic pp cross section ((Tj„) to determine the 
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instantaneous luminosity using the following equation: 



L = (3.3) 

where L is the instantaneous luminosity, fsc is the rate of bunch crossings in 
the Tevatron and // is the average recorded number of pp interactions per bunch 
crossing. 

This analysis uses SVX-II, COT, CMU and the trigger, which will be described 
in detail in the following sections. More detailed information about each sub- 
detector and the trigger may be found in Bishai [41]. We refer the reader to Af- 
folder and Hill [42,43] for a full documentation about the ISL and LOO, Acosta [44] 
about the TOF, Balka, Bertolucci and Kuhlmann [45-47] about the calorimeters, 
Artikov [48] about the CMP, the CMX and the IMU, and Acosta [49, 50] about 
the CLC. 



3.3 Tracking Systems 

3.3.1 Silicon Vertex Detectors II 

Silicon tracking detectors are used to obtain precise position measurements of 
the path of a charged particle. They present some advantages over the gas filled 
drift chamber. The typical electron-hole creation energy of the sihcon is about 
3 eV, while the ionization energies arc about 10-15 eV for the drift chamber gas 
(Argon or Ethane). The increased number of electrons gives better energy and 
position resolutions and signal to noise ratio. The fundamental component of 
the sihcon detector is a reverse-biased p-n junction. The reverse-biased voltage 
increases the gap between the conduction band and the valence band across the 
p-n junction and reduces the current from the thermal excitation. By segmenting 
the p or n side of the junction into "strips" and reading out the charge deposition 
separately on every strip, we obtain sensitivity to the position of the charged 
particle. More information about the principles of silicon detector may be found 
in Knoll [51]. 

The CDF SVX-II is composed of double-sided silicon sensors. Each side mea- 
sures either the (f) or z coordinates. The strips on the side run axially, while 
the strips on the z side run either perpendicular to the axial strips, or are tilted 
by 1.2°. On the (j) measurement side, 65 /xm pitch strips of p-type silicon are 
implanted near the surface of a mild n-type (high purity) silicon bulk. On the z 
measurement side, strips of n-type silicon are implanted on the same high purity 
bulk. 

Four silicon sensors are assembled on a ladder. The readout electronics are 
mounted directly to the surface of the silicon sensor at each end of the ladder, 
as shown in Figure 3.7. The ladders are arranged in three, 29 cm long barrels 
in an approximately cylindrically symmetric configuration. The barrels arc po- 
sitioned end-to-end along the beam axis as shown in Figure 3.8, and segmented 
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Tabic 3.1: Relevant parameters for the layout of the sensors of different SVX-II 
layers. 



Property 


Layer 


Layer 1 


Layer 2 


Layer 3 


Layer 4 


number of strips 




256 


384 


640 


768 


869 


number of Z strips 




256 


576 


640 


512 


869 


stereo angle 


(°) 


90 


90 


+1.2 


90 


-1.2 


strip pitch 


(/xm) 


60 


62 


60 


60 


65 


Z strip pitch 


(/xm) 


141 


125.5 


60 


141 


65 


active width 


(mm) 


15.30 


23.75 


38.34 


46.02 


58.18 


active length 


(mm) 


72.43 


72.43 


72.38 


72.43 


72.43 



azimuthally into 12 wedges (30° each). Each wedge contains 5 layers of silicon 
sensors. Table 3.1 gives the mechanical dimensions of SVX-IL The "impact 
parameter" is defined as the distance from the closest approach of the particle 
trajectory to the primary vertex. Without ISL and LOO, SVX-II reaches a impact 
parameter resolution of ~ 50 /xm, which includes the 30 \xm. contribution from 
the transverse size of the beam line. 

3.3.2 Central Outer Tracker 

The Central Outer Tracker (COT) is a cyhndrical, open-cell drift chamber. The 
active volume of the COT begins from r=43.4 cm to r=132.3 cm and spans 
310 cm in the z direction. The COT contains 96 sense wire layers in the radial 
direction, which are grouped into eight superlayers, as shown in Figure 3.9. The 
maximum drift distance is approximately the same for all superlayers. 

Each superlayer is divided in into "cells". Figure 3.10 shows the transverse 
view of three COT cells. In each cell, 12 sense wires and 17 potential wires are 
sandwiched by two grounded field (cathode) sheets. The potential wires help 
shape the electric field near the sense wires. The distance between the wires and 
each field sheet is 0.88 cm. The sense wires alternate with the potential wires at 
a pitch of 0.3556 cm. Each end of the cell is closed with two potential wires, the 
first at the same pitch of 0.3556 cm, and the second at the pitch of 0.1778 cm. 
The wires are made of 40 /xm of gold plated tungsten and the field sheets are 
made of 6.35 /xm Mylar with vapor-deposited gold on each side. The entire COT 
contains 2,520 cells and 30,240 sense wires. The wires in superlayer 1, 3, 5 and 7 
run along the z direction ("axial"). The wires in the other superlayers are strung 
at a small angle (2°) with respect to the z direction ("stereo"). 

The COT chamber is filled with an Argon-Ethane gas mixture and Isopropyl 
alcohol (49.5:49.5:1). The voltages on the sense wires are 2600-3000 volts and 
1000-2000 volts on the potential wires. The field sheets are grounded. These 
voltages have been optimized using Garfield simulation software [52] in order to 
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achieve a gas gain of 2 • 10^ and minimize the drift field deviations from wire to 
wire. When a charged particle passes through, the gas is ionized. The apphed 
electric field in the cell allows us to collect the electrons from the ionization and 
to track the passage of the charged particle. Electrons drift towards the sense 
wires with a velocity ~ 50 /xm/ns. Due to the magnetic field in which the COT is 
immersed, electrons drift at a Lorentz angle of ~ 35°. Therefore, each cell is tilted 
by 35° with respect to the radial direction to compensate for this effect. After 
the tilt, the electrons drift approximately perpendicular to the radial direction. 
The maximum electron drift time is approximately 177 ns, much shorter than 
the beam crossings of 396 ns. The electric potential in a cylindrical system 
grows logarithmically with decreasing radius. As a result, a limited avalanche via 
secondary ionization is initiated when the electron drifts close to the wire surface. 
This effect provides us a gain large enough to be read out by the electronics with 
high signal to noise. 

Signals on the sense wires are read out by the ASDQ (Amplifier, Shaper, 
Discriminator with charge encoding) chip, which was developed by Bokhari and 
Newcomer [53]. The ASDQ provides input protection, amplification, pulse shap- 
ing, baseline restoration, discrimination and charge measurement. The analog 
signal arrives at the ASDQ and the output is a digital pulse. The leading edge 
gives the arrival time information and the pulse width is related to the amount 
of charge collected by the wire. After calibrating the width variations due to the 
COT geometry, path length of the particle, gas gain difference for the 96 wires, 
the digital width is related to the ionization energy loss dE/dx, used for particle 
identification. A detailed description of the calibrations performed by the author 
and Donega, Giagu, Tonelli is found in Yu [54] and Donega [55]. 

The ASDQ pulse is then sent through ~ 105 cm of micro-coaxial cable, via 
repeater cards to Time to Digital Converter (TDC) boards in the collision hall. 
Hit times are later processed by pattern recognition (tracking) software to form 
helical tracks. The hit resolution of the COT is about 140 /^m. The transverse 
momentum (Pt) resolution, ap^/Pr is about 0.15% • Pt- Table 3.2 lists the COT 
parameters. A full documentation about the COT may be found in Affolder [42]. 

3.3.3 Track Reconstruction 
Definition of Track Parameters 

Charged particles moving through a homogeneous solcnoidal magnetic field in 
the —z direction follow helical trajectories. To uniquely parameterize a helix 
in three dimensions, five parameters are needed: C, cot^, do, 0o ^-nd Zq. The 
projection of the helix is a circle on the r-0 plane. C is the signed curvature 
of the circle, defined as C = where p is the radius of the circle and the 
charge of the particle (Q) determines the sign of C. The positive charged tracks 
curve counterclockwise in the r-cp plane when looking into the —z direction and 
the negative charged tracks bend clockwise. The transverse momentum, Pt, is 
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Figure 3.9: Layout of wire planes on a COT endplate. 




Figure 3.10: Transverse view of three COT cells. Shown are cathode field sheets 
(solid lines), 12 sense wires (circles with cross) and 17 potential wires (marked 
by open circles) in each cell. 
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Table 3.2: COT Parameters. 



Parameter 


Valnc 


Gas 


Ar/Et/Isopropyl(49.5:49.5:l) 


Max. Drift distance 


0.88 cm 


Max. Drift Time 


177 ns 


Lorentz Angle 


~ 31° (35 °cell tilt) 


Drift Field 


1.9 kV/cm 


Radiation Lengths 


1.7% 


Total sense wires 


30,240 


Number of cells per SL 


168,192,240,288, 




336,384,432,480 


Stereo Angle 


+2°, 0°, -2°, 0° 




+2°, 0°, -2°, 0° 



related to C, the magnetic field {Bmagnet), and charge of the particle: 

P _ 1A9898 ■ 10-^ ■ Bmagnet 

^T-Q ^ , I3.4j 

The 6 is the angle between the z axis and the momentum of the particle. 
Therefore, cot 6 is Pz/Pt, where Pz is the z component of the particle momentum. 
The last three parameters, do, (po, and zq, arc the r, and z cylindrical coordinates 
at the point of closest approach of the helix to the beam line. See Figure 3.11 for 
the definition of do and 0o- ^^o is a signed variable; 



do^Q-{^xl + yi-p), (3.5) 

where {xq, yo) is the center of the helix circle in the r-0 plane. Figure 3.12 illus- 
trates the sign definition of do. 

For decaying particles, we define the displacement L^y, 

L.y = d- Pt (3.6) 

where d is the displacement of the decay vertex in the transverse plane, and Pt 
is the unit vector in the direction of Pt- 

Pattern Recognition Algorithms 

The track reconstruction begins using only the COT information. The first step is 
to look for a circular path in the axial superlayers of the COT. The algorithm looks 
for 4 or more hits in each axial superlayer to form a straight line, or "segments" . 
The hits on the segment are reconstructed using the time difference between when 
the ionization occurs, to (the collision time plus the time of flight of the charged 
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Figure 3.12: Tracks of particles with positive/negative charge and posi- 
tive/negative impact parameters. 
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particle) and when the signal is picked up by the wire (the leading edge time 
of the digital pulse from the TDC). The global time offset, readout time of the 
wires and cables, electronic channel pedestals, charged-based time slewing and 
non-uniform drift velocities are corrected before using the time difference (or drift 
time) in the tracking. 

Once segments are found, there are two approaches to track finding. One 
approach is to link together the segments which are consistent with lying tangent 
to a common circle. The other approach is to constrain its circular fit to the 
beamline, and then add hits which are consistent with this path. Once a circular 
path is found in the r — (p plane, segments and hits in the stereo superlayers are 
added depending on their proximity to the circular fit. This results in a three- 
dimensional track fit. Typically, if one algorithm fails to reconstruct a track, the 
other algorithm will not. 

Once a track is reconstructed in the COT, it is extrapolated into the SVX-II. 
A three-dimensional "road" is formed around the extrapolated track, based on 
the estimated errors on the track parameters. Starting from the outermost layer, 
and working inward, silicon clusters found inside the road are added to the track. 
As a cluster is added, a new track fit is performed, which modifies the error 
matrix for the track parameters and produces a narrower road. In the first pass 
of this algorithm, r — clusters are added. In the second pass, stereo clusters are 
added to the track. If there is more than one track with different combinations of 
SVX hits associated with the same COT track, the track with maximum number 
of SVX hits is chosen. 

The track reconstruction efficiency in the COT is ~ 95% for tracks which 
pass through all 8 superlayers {Pt > 400 MeV/c ) and ~ 98% for tracks with 
Pt > 10 GeV/c. The SVX track reconstruction efficiency with the COT tracks 
in the denominator is about 93% for the tracks with at least 3 SVX r-cf) hits. A 
complete description of the COT and the SVX tracking is found in Hays [56]. 

3.4 Central Muon Detector 

The Central Muon Detector (CMU) is embedded in the central hadron calorime- 
ter wedges at r—3A7 cm as shown in Figure 3.13 and covers r] < 0.6. The detector 
is segmented in (j) into 12.6° wedges, while the calorimeter is segmented into 15° 
wedges. This leaves a 2.4° gap between each wedge. In addition, there is a gap 
between the cast and west chambers at r]=0. The detector is further segmented 
in (j) into three 4.2° modules. There are 72 modules at the east and west ends of 
the detector, which gives 144 modules in total. Each module has 4 layers of 4 
rectangular drift cells. The dimension of the drift cell is 6.35 cm (x) x 2.68 cm 
(y) X 226.1 cm (z). Each cell has a 50/um stainless steel sense wire in the center. 

The first and the third layers have small offset in the direction from the 
second and the fourth layers, which also means: for each (p module, two sense 
wires from the alternating layers ( 1 and 3 or 2 and 4) he on a radial fine. The 
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Figure 3.13: Location of the Central Muon Detector (CMU). 



other two sense wires he on a line with a 2 mm offset from the first two. The 
ambiguity as to which side of the sense wire (in 0) a track passes is resolved 
by determining which two layers of sense wires are hit first. The sense wires of 
alternating cells in the same layer are connected together so to enable readout 
at only one end of the chamber. Signals from the sense wires are discriminated 
and passed on to the same type of TDCs as used by the COT. Short tracks re- 
constructed using the TDC and ADC information are referred to as "stubs". The 
muon stubs are matched to the reconstructed tracks to form a muon candidate. 
A is computed using the distance between the track and the stubs, the differ- 
ence in the direction of the track and the stub, and the covariance matrix of the 
track. Full documentation about the CMU is found in Ascoli [57]. Table 3.3 lists 
the parameters of the CMU, where the pion interaction lengths and the multiple 
scattering are computed at a reference angle of 6* = 90°. 



3.5 Triggers 

The triggering systems play an important role in the pp collider for two reasons. 
First, the collision rate is about 2.5 MHz, which is much higher than the rate 
at which data can be stored on tape, 50 Hz. Second, the total hadronic cross- 
section (including the elastic, inelastic, and diffractive processes) is about 75 mb 
and the bb cross-section is about 1000 times smaller, 0.1 mb. Extracting the most 
interesting physics events from the large number of events reduces the cost and 
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Table 3.3: Parameters of the CMU. 



Parameter 


Value 


7] coverage 


< 0.6 


Drift tube cross-section 


2.68 X 6.35 cm 


Drift tube length 


226.1 cm 


Max drift time 


800 ns 


Total drift tubes 


2304 


Pion interaction lengths 


5.5 


Min detectable muon Pt 


1.4 GeV/c 


Multiple scattering resolution 


12cm/p 



time to reconstruct data. 

The goal of the CDF-II triggering system is to be dead-time-less so that the 
system is quick enough to make a decision for every single bunch crossing be- 
fore the next bunch crossing occurs. Each level of the trigger must reduce the 
background to a low enough level so that the rate to the next level is not satu- 
rated. Each level of the trigger is given an amount of time to make a decision 
about accepting or rejecting an event which depends on the complexity of the 
reconstruction. At the first level (Level 1), a trigger decision is made based 
only on a subset of the detector and quick pattern recognition or simple count- 
ing algorithms. The second level of the trigger (Level 2) does a limited event 
reconstruction. The third level of the trigger (Level 3) uses the full detector in- 
formation to fully reconstruct events in a processor farm. The decision time for 
Level 1, 2 and 3 is about 5.5 fis, 20 fis and 1 s, respectively. The event accept 
rate for Level 1, 2 and 3 is 40 kHz, 300 Hz and 50 Hz. The delay necessary to 
make a trigger decision is achieved by storing detector readout information in a 
storage pipeline, as shown in Figure 3.14. 

A trigger path is a well defined sequence of Level 1, Level 2 and Level 3 
triggers. An event coming from one trigger path satisfies the trigger requirements 
at each level. A well defined trigger path eliminates volunteer events. A volunteer 
event is an event which passes a higher level trigger requirement but passes a 
lower level trigger from a different path. At CDF II, there are about 100 trigger 
paths. The trigger path used in this analysis is one of the SVT trigger paths, 
B_CHARM Scenario A. Figure 3.15 shows a general diagram of the SVT trigger 
path. 

The strategy of the SVT trigger path is as follows. At Level 1, the eXtremely 
Fast Tracker (XFT) measures the track Pt and angle 0. By cutting on Pt and 
(f), most of the inelastic background will be rejected. The Extrapolation Unit 
(XTRP) selects the XFT tracks above a certain P^ threshold and sends signal to 
the Level 1 Track Trigger (Level 1 Track). The Level 1 Track Trigger counts the 
number of tracks from the XTRP, if more than 6 tracks are found, an automatic 
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Figure 3.14: Diagram of the CDF-II trigger system. 
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Figure 3.15: Diagram of the two track trigger path at Level 1 and 2. 
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Level 1 accept is generated. Otherwise, depending on the trigger requirements, 
Level 1 Track Trigger accepts or rejects the event. If a Level- 1 accept is received, 
the XFT track information is sent to the Level 2 Silicon Vertex Trigger (SVT). 
At Level 2, SVT uses the SVX-II information to obtain impact parameters of the 
tracks, do. Requiring non-zero impact parameters of tracks will require that they 
come from decays of long-lived particles: charmed and bottom hadrons. The 
requirements of each level of trigger will be described in detail in Section 3.5.4. 
The trigger components used in this trigger path, XFT, SVT, and Level 3, will be 
discussed in the following text. We refer the reader to the CDF Run II Technical 
Design Report [58] for a full descriptions of the trigger hardware. 

3.5.1 The extremely Fast Tracker (XFT) 

The Level 1 trigger decision of the "two track" trigger path is based on the 
information from the eXtremely Fast Tracker (XFT) . This device is designed to 
measure the momentum of the charged particle using the hit information of the 
4 COT axial layers. Instead of using the TDC information and a drift model to 
find a track segment as described in Section 3.3.3, the XFT uses a fast binary-like 
algorithm. 

Each hit on the wire is classified as "prompt" if the drift time ranges from 
to 66 ns and as "delayed" if the drift time ranges from 67 to 220 ns. Four 

neighboring COT cells are grouped together when searching for a segment in a 
given superlayer. A track segment in each axial superlayer is found by comparing 
the hit patterns to a list of pre-loaded patterns. The hit pattern varies depending 
on the combination of delayed and prompt hits, and the track angle through the 
cell or the track Pt. The algorithm allows two missed hits (2-miss) in each 
segment for the beginning period of the data used for this analysis and tightens 
the requirement to one missed hit (1-miss) since October, 2002. The data used 
in this analysis in the 2-miss period is about 26.4 pb~^, and 124.5 pb~^ for the 
1-miss period. 

Once a segment is found in a superlayer, it it marked as a "pixel". XFT 
compares the pixels in all 4 layers to a list of pixel patterns in a A(J) = 1.25° 
window corresponding to a valid track with Px > 1.5 GeV/c (~ 2400 roads). The 
algorithm returns the best track. Figures 3.16-3.17 extracted from Thomson [59] 
show an example of the hit and the track pattern for a track with Pt — 1.5 
GeV/c. Finally the XFT reports the track Pt and 06, the angle of the transverse 
momentum at the sixth superlayer of the COT, which is located 106 cm radially 
from the beamline. 

The XFT efficiency is > 96% for tracks which pass through all 4 axial layers. 
The momentum resolution, gp^/Pt is about 2% per GeV/c. The angular resolu- 
tion at the sixth superlayer, cr^g is about 5 mR. More detailed information about 
XFT may be found in Thomson [59]. 
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Figure 3.16: Close view of a track with Pt = 1.5 GeV/c in a cell of the 4*^^ COT 
axial layer from Thomson [59]. A collection of the prompt hits (marked by open 
circles) and the delayed hits (marked by open diamonds) is an example of the 
XFT hit pattern. 
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Figure 3.17: Close view of a track with Pt — 1.5 GeV/c traversing 4 COT axial 
layers from Thomson [59]. Several possible track patterns are displayed. 
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3.5.2 The Silicon Vertex Tracker (SVT) 

At Level 2, the Silicon Vertex Tracker (SVT) combines the Level 1 track infor- 
mation, computes the track do, and improves the measurements of Bp and 00 . 
The SVT is a new type of Level 2 trigger optimized for B physics . The heavy 
flavor particles, such as beauty and charm hadrons, decay at positions displaced 
from the primary vertices. Therefore, their daughter tracks tend to have larger 
impact parameters. The SVT cuts on the minimum track impact parameter so 
we are able to collect large sample rich in heavy flavor. The SVT also cuts on 
the maximum track impact parameter and removes background due to the long 
lived Ks, A or the secondary tracks which come from the particles interaction 
with the beam pipes. 

As mentioned in Section 3.3.1, the SVX-II is segmented into 12 wedges in 
and three mechanical barrels in z. The SVT makes use of this symmetry and does 
tracking separately for each wedge and barrel. The XFT track is extrapolated 
into the SVX-II, forming a "road". Clusters of charge on the inner four r — <p 
layers of the given wedge have to be found inside this road. Since June 2003, 
the requirement is loosened to ask for hits from any four r-0 layers out the five 
SVXII layers. This period of data corresponds to about 1/3 of the total integrated 
luminosity used for this analysis. The SVT checks if one of the roads in the list 
of pre-loaded patterns is present in the data. The found roads are fed into a 
linearized fitter which returns the measurements of pt, 0o and do for the track. 
The Level 1 trigger conditions are confirmed with the improved measurements of 
Pt and 0o- An event passes Level 2 selection if there is a track pair reconstructed 
in the SVT and additional requirements on the do, Ly.y, and scalar sum Pt of the 
track pair depending on the trigger path. 

Figure 3.18 shows the SVT track impact parameter resolution for tracks with 
> 2 GeV/c. The width of the Gaussian fit for the distribution in Figure 
3.18 is 47 //m, which is a combination of the intrinsic SVT impact parameter 
resolution, and the transverse size of the beam line: ajit = cfsvt ® o-heam, where 
o'beam is about 30 /im. Therefore, the intrinsic SVT resolution is about 35 /xm. 
Full documentation about SVT is found in Ashmanskas [60]. 

3.5.3 Level 3 Trigger 

The third level of the trigger system is implemented as a Personal Computer 
(PC) farm. The input rate of the Level 3 is roughly 300 Hz. With roughly 300 
CPUs, and one event per CPU, this allocates approximately 1 second to do event 
reconstruction and reach a trigger decision. 

Figure 3.19 shows the principle of the Level-3 farm. The detector readout from 
the Level 2 buffers is received via an Asynchronous Transfer Mode (ATM) switch 
and distributed to 16 "converter" node PC's (CV). The main task of these nodes 
is to assemble all the pieces of the same event from the different sub-detector 
systems. The event is then passed via an Ethernet connection to a "processor" 
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Figure 3.18: SVT impact parameter resolution. 

node (PR). Each processor node is a separate dual-processor PC. Each of the 
two CPU's on the node process a single event at a time. The processor runs a 
"filter" executable which performs the near-final quality reconstruction. If the 
executable decides to accept an event, it is then passed to the "output" nodes 
of the farm (OU). These nodes send the event onward to the Consumer Server / 
Data Logger (CS/DL) system for storage first on disk, and later on tape. A full 
description of the Level 3 system is found in Gomez-Ceballos [61]. 

For the first two thirds of the data used in this analysis, the COT track 
reconstruction algorithms as described in Section 3.3.3 are performed at Level 
3. The COT tracking returns Pt,Zq,(I)q and cot^ and combines with the do 
measurement from the SVT to create a further improved track. The Level 1 
and Level 2 trigger conditions are confirmed at Level 3 using improved track 
measurements. For the last one third of the data, full SVX-II tracking is available, 
and the trigger conditions are repeated using a combined COT/SVX-II fit of the 
track helices. 

3.5.4 B_CHARM Scenario A Trigger Path 

The B_CHARM Scenario A trigger path requires two tracks from the SVT and cuts 
on the Pt, the minimum and maximum do, the differences between the tracks' 
parameters, such as A0o, ^zo, and L^y. A0o is defined as the opening angle 
between the track pair in the r-0 plane. Azo is the Zq difference of the track pair. 
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Figure 3.19: Event building and Level 3 operating principle: data from the front 
end crates is prepared by Scanner CPU's (SCPU) and fed into the ATM switch. 
On the other side of the switch, converter nodes (CV) assemble events and pass 
them to processor nodes (PR). Accepted events are passed to output nodes (OU) 
which send them to the Consumer Server and Data Logging systems (CS/DL). 



Lxy is the distance of the intersection of two tracks with respect to the SVT beam 
line projected on the direction of the total momentum vector in the r-0 plane. 
The cuts at Level 1-3 triggers of B_CHARM Scenario A trigger path are described 
below. 

The trigger requirements are: 
Level 1 

• a pair of opposite charged XFT tracks 

• each XFT track transverse momentum Pt > 2.04 GeV/c 

• scalar sum + > 5.5 GeV/c 

• A0O < 135 ° 

Level 2 

• a pair of opposite charged SVT tracks 

• each SVT track satisfies: 

- SVT track fit < 25 

- Pt>2 GeV/c 

- SVT impact parameter: 120 /im < do{SVT) < 1000 /im 
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• scalar sum Pj, + P^ > 5.5 GeV/c 

• 2° < A0O < 90° 

• L^y > 200 /im, this cut was added starting with run 150010 [62] 
Level 3 

The following cuts do not change for the whole period of taking data: 

• each Level 3 track: Pt > 2 GeV/c and \r]\ < 1.2 

• scalar sum P^ + P^ > 5.5 GeV/c 

• 2° < A(/)o < 90° 

• L^y > 200 iim 

• Azo < 5 cm 

The following cuts are different for period I and II. All the number of time 
integrated luminosities are after the good run selection. 

Period I: 9*^ February 2002 to 19*^ May 2003, Runs 138809-163113, ~ 120 pb^^. 

• Tracking algorithm at Level 3 uses only the COT hits and Level 3 tracks 
should be matched to Level 2 tracks found by the SVT 

- Track azimuthal angle difference: A(0ol3 — ^osvt) < 0.015 radians 

- Curvature difference: A{Cl3 - Csvt) < 0.00015 cm"^ 

• 120 yum < (io(-^3) < 1000 ^m., is calculated using SVT beamline 

Period II: 19*^^ May 2003 to 6*^^ September 2003, Runs 163117-168889, ~ 50 pb'^ 

• Tracking algorithm at Level 3 uses both the SVX and COT hits 

— require at least 3 hits in different SVX layers 

— no attempt to match Level 3 tracks to SVT tracks 

• 80 //m < do (-^3) < 1000 //m, is calculated using SVT beamline 

Data derived from the above trigger path were written to the tape for subse- 
quent reconstruction and physics analysis. 
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Chapter 4 
Data Samples 



Data used in this analysis are collected with the upgraded CDF detector from 9*'* 
February 2002 to 6*'* September 2003 and cover runs 138809 through 168889. This 
period corresponds to an integrated luminosity of ~237 pb~^. In this chapter, we 
present details of how we arrive at our final data sample, optimize the cuts. We 
describe the data sample used for this analysis in Section 4.1. The cuts that are 
applied to obtain a clean signal with low background are explored in Section 4.2. 

4.1 Data Sample 
4.1.1 Overview 

We wish to reconstruct the decays of the and the A^. The cr of the and 
Aft are about 460 and 370 //m, respectively. The cr of and A+ 

are about 120, 150, 300 and 60 iim. Long hved Kg and A have much longer 
lifetimes, 2-8 cm. Therefore, a minimum requirement on the distance between 
the beamline and the secondary vertex (decay length) removes contamination of 
short-lived charm hadrons. The daughter particles from the B decays also tend 
to have a larger impact parameter {do) than the tracks produced at the primary 
vertex, as illustrated in Figure 4.1, but smaller do than the daughter particles 
of Kg and A. Consequently, a minimum and a maximum cut on the track do 
rejects the background from the primary tracks. Kg and A, or secondary tracks 
from the particle interaction with the detector material. The newly developed 
SVT trigger cuts on the decay length and the track do- It is the best trigger for 
distinguishing B decays from other physics processes. Among all the SVT trigger 
paths, we find B_CHARM Scenario A most suitable for this analysis. Section 3.5.4 
presents the definition of the B_CHARM Scenario A trigger path. The goal of any 
measurement is to have the smallest possible uncertainty. By recording data from 
the same trigger, the systematic uncertainties common to both modes cancel. We 
therefore select all our data from the B_CHARM Scenario A trigger path. 

The data from B_CHARM Sceneirio A are processed with the Production ex- 
ecutable, version 4.8.4, and compressed into the secondary datasets hbotOh and 
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Figure 4.1: Example of a decay: vr"^ from the decay has a larger impact 
parameter, do, than that of the track produced at the primary vertex, tt^. 

hbotli. The total size of hbotOh and hbotli is about 10 Terabytes (150M 
events), which is too big to be analyzed quickly multiple times. We apply loose 
selection cuts and reduce hbotOh, hbotli to smaller, tertiary datasets. Sec- 
tion 4.1.2 discusses the data skimming. Then we optimize the analysis cuts using 
the tertiary datasets in Sections 4.2. 

4.1.2 Data Skimming 

Before skimming the data, we exclude the runs with an incorrect alignment table 
(152595-154012) in hbotOh. The alignment table contains the parameters for 
the positions of the COT and the SVX. The data during runs 152595-154012 
are reprocessed with the correct alignment table and collected into hbotli. We 
further require the following systems declared good by the CDF Data Quality 
Monitoring group: SVX, COT, CMU, Cherenkov Luminosity Counters (CLC) 
and the Level 1-3 triggers. We exclude the runs when SVX is off and when there 
are known high voltage or trigger problems in the COT, CMU or SVT. By making 
these requirements, the Monte Carlo program can better reproduce the response 
of these detectors (see Section 6.1). After making the good run selection, the 
integrated luminosity reduces from 237 pb~^ to 171.5 pb~^. 

The skimming program starts by storing a set of offline reconstructed tracks 
which satisfy the quality requirements on: Pt, the number of COT hits in the 
axial and stereo layers, the number of SVX r-0 hits, and the impact parameter. 
Then, the tracks that are matched to those found by the SVT or to muon stubs 
(CdfMuon) in the CMU, are marked for further use. After saving the SVT 
and muon information, we begin our reconstruction by identifying the charm 
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candidates: D'^ and A+. 

We first cut on the raw mass of tlie cliarm candidates, wlierc tlie raw mass 
is calculated using the track momentum at the point of closest approach to the 
beam line. We determine the charm (tertiary) vertex by performing a vertex fit 
with the CTVMFT package developed by Marriner [63]. CTVMFT determines 
the decay vertex by varying the track parameters of the daughter particles within 
their errors, so that a between the track trajectory and the points is minimized. 
We cut on the fitted charm mass and X^_<^, where "x^-</>" the returned from 
the fit in the r-0 plane. 

The charm candidate is then combined with an additional track to form the 
B candidate. The additional track has a minimum Pt requirement of 1.6 GeV/c. 
Once we have a valid fourth track, we cut on the raw mass of the four tracks. The 
mass window varies depending on whether the fourth track is matched to a muon 
stub. We require that two of the four tracks from the reconstructed 5-hadron 
candidate each matches an SVT track. We confirm the trigger by requiring the 
matched SVT tracks to pass the Scenario A cuts listed in Section 3.5.4. We then 
perform a four-track vertex fit. The four-track vertex fit includes a constraint 
that the tertiary vertex points to the secondary vertex in the r-(j) plane. After 
the vertex fit, we cut on Pt of the charm, the xl-cj>j Pt and fitted mass of the 
four tracks. 

After applying the requirements discussed above for each signal mode, we 

reduce the secondary datasets, hbotOh and hbotli, by a factor of 25 (from 10 
to 0.4 Terabytes). The reduced datasets are then written to tape for further 
use. Detailed information about the skimming is found in the reference by the 
author et al. [64]. In Section 4.2, we present our analysis cut optimization with 
the reduced datasets. 

4.2 Signal Optimization 

From the reduced datasets we reconstruct our signals: 

• S° ^ D*+7r- and B D*+ii-X, where D*+ D°7r+, D° K-7r+ 
m B^ ^ D+Ti- and B D+fx-X, where D+ K-n+Tr+ 

• A+Tr~ and B A+ii'X, where A+ pK-n+ 

The reconstruction procedure is similar to that described in Section 4.1.2 and 
Yu [64] . The following cuts are studied more carefully and optimized : 

• xl-4> of B and charm vertex fit 

• Pt oi B and charm candidates 

• CT ol B and charm candidates: L^.^ x 
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Figure 4.2: Signal optimization: number of signal in the MC and data as a 
function of the ct{B) cut after applying the normalization factor for A;, 



Our semileptonic signals are larger than the hadronic signals, and the statisti- 
cal uncertainty of the relative branching fraction measurement is dominated by 
the uncertainty of the number of events in the hadronic signals. Therefore, we 
optimize the hadronic mode only and apply the optimized cuts to the semilep- 
tonic mode. The optimized quantity is the significance, -r^^, where "S"' is the 
number of signal and "-B" is the number of background events. 

For our optimization, the amount of signal, "5"' comes from a MC as described 
in Section 6.1. The reason for using MC signal is that the data signal is small 
and susceptible to statistical biases. We generate MC with about 20 times more 
events than the data and eliminate this problem. In order to scale the significance 
close to the true value measured from the data, we apply a normalization factor 
fc on the signal MC, 

fc = 1^, (4.1) 



^1 



MC 



where S'data and ^mc are the amount of the signal found in the data and MC 
after applying loose cuts, and 



S = fcXS, 



MC- 



(4.2) 



Note that even though "S"' , after scahng, is the same as S'data now, the uncertainty 
on "S"' is fc X ^/SMC, is much smaller than the uncertainty of S'data, V "S'data- 
Figure 4.2 shows a comparison of the number of signal in the data and in the MC 
after applying the normalization factor. 

We evaluate the background beneath the signal peak from the data. We first 
apply loose cuts on each mode to identify a clear 5° or A^ peak; 

• cr{B) > 50 yUm 

• each track Pt > 0.5 GeV/c 
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Figure 4.3: AcTt invariant mass for data and MC. Left: data, the pink solid lines 
indicate the signal region. Right: MC. 

• 7T from the B hadron is CMU fiducial 

• for 5^ D*+TT-: 

- 1.833 < Mk^ < 1.893 GeV/c^ 

- 0.143 < MKnn - < 0.148 GcV/c^ 

• for # ^ D+7V-: 1.8517 < M^^^ < 1.8837 GeV/c^ 

• for Afc ^ A+TT-: 2.269 < MpK^ < 2.302 GeV/c^ 

We require that both the muon and pion from the B hadron point within CMU 
fiducial volume because we use the CMU only to identify the muons. CMU 
covers the region of pseudo-rapidity (77) less than 0.6. Making the same fiducial 
requirement for the hadronic mode allows the tracking efficiencies from both 
modes to cancel. 

The backgrounds in the signal and in the upper mass regions are mainly 
combinatorial, and may be described by an exponential function, as we will see 
in Section 5.2. Therefore, we fit the upper mass region to an exponential function. 
Finally we extrapolate and integrate the exponential over the mass region of ± 3 a 
around the signal peak to obtain "5" . Figure 4.3 shows the Ab mass distribution 
in the data and MC. The figure also shows the signal region we define and the 
upper mass region we fit to an exponential. 

The optimization follows an iterative procedure which passes through the 
data multiple times. In the first pass, cuts on each variable are scanned and 
optimization points are found. In the second pass, we apply the optimized cuts 
for all but the variable which is being re-optimized. We iterate this process 
several times until the optimization points become stable; usually twice is enough. 
Figures 4.4 shows ^r^^ , § and from the optimization of A^ A+vr" mode, 
as a function of each cut variable, where S'ref is the number of signal events at 
the starting point. Tables 4.1- 4.2 list the final analysis cuts. Note that because 
the MC and the data do not agree well, as shown in Section 6.2, we choose 
to make a loose cut at the plateau region of the significance. The final analysis 
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cuts for the Pt of charm hadrons are tighter than the optimization points. The 
tighter cuts arise from the 4 GeV/c Pt threshold apphcd to the c-quark in the 
MC sample for our semileptonic background study (see Section 7.4.2). This Pt 
threshold makes the reconstruction of charm hadrons below 4 GeV/c inefficient. 
The MC sample is produced by the CDF B group and it would take a prohibitive 
amount of CPU time to generate a new sample more suitable for our analysis. 
Therefore, we increase the Pt cut of our charm hadrons to 5 GeV/c. As the 
significance of the charm is a slowly varying curve, changing the cuts has little 
effect on the signal yield. 

In addition to the cuts which are optimized above, we also require that the 
muon and pion from the B hadron each matches an SVT track. Finally, for the 
semileptonic modes, we make cuts on the four track invariant mass (eg: M(Ac/i)) 
to reduce the backgrounds from the other B decays, see Section 7.2 for more 
details. The signal and sideband distribution of each optimized variable after 
N — 1 cuts can be found in Yu [65]. The optimization yields a, S/B oi 37.6 
and 62.8 for the S° D*+7r- and B D*+n-X modes, 2.6 and 1.3 for the 

— > D+TT" and B — > D'^iJ,''X modes, 1.6 and 0.3 for the A;, — >■ A+tt" and 
B — > A^yU~X modes. Figure 4.5 shows the charm+7r (left) and charm (right) 
mass spectra from the hadronic and inclusive semileptonic signals in the data 
after applying the optimized analysis cuts. 

4.3 Summary 

We have reconstructed our signals in the data collected from the trigger path 
B_CHARM Scenario A. We have optimized our analysis cuts. In the next chapter, 
we will present the fit to the charm and B hadron mass spectra to obtain the 
number of signal events. 
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PT(A,)[GeV/c'] t^^fB) [cm] 




Figure 4.4: Significance (y^j) niarked by circles, signal/background (^) marked 
by squares and signal/reference {-^) marked by triangles for cuts used in — >• 
A^X analysis. 
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Figure 4.5: charm+vr and charm mass spectra from our signals after all cuts. 
From the top left to the bottom right are: M/j*^, M£)0^ — Mj^o, M^)^, Mk-k-k, 
Ma^tt and MpKn- 
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Table 4.1: Final analysis cuts shared by all the modes. 



All 


Pt for all tracks 


> 0.5 GeV/c 


ttb and /ib Pt 


> 2.0 GeV/c 


Pt of 4 tracks 


> 6.0 GeV/c 


Pt of charm hadron 


> 5.0 GeV/c 


fXB CMU xl 


< 9 


every track exits at COT layer 95 


TiB and hb matched to SVT tracks and CMU fiducial 
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Table 4.2: Final analysis cuts for each mode. 

B D*+X 
VertexFit x^-^ < 16 
4 track VertexFit xl-4> < 17 
cr(D" B) > -70 /xm 

ct{B beamspot) > 200 //m 
1.833 < Mk^ < 1.893 GeV/c^ 
3.0< Mk^^^ <5.3 GeV/c2 for B D*+fi-X 
0.143< Am <0.148 GeV/c^ for ^ D*+7r- 

B^ D+X 
D+ VertexFit Xr-^ < 14 
4 track VertexFit x^_^ < 15 
ct(D+ ^ S) > -30 urn 

ct{B — > beamspot) > 200 //m 
3.0< Mx^^/, <5.3 GcV/c2 for B D+fx-X 
1.8517 < Mk^-k < 1-8837 GeV/c^ for ^ D+tt" 

Ab ^ 

Pt of proton > 2 GeV/c 

A+ VertexFit Xr-<^ < 14 

4 track VertexFit Xr-4, < 15 

cr(A+ ^ Ab) > -70 //m 

CT(Ab^ beamspot) > 250 /im 

3.7< MpK^i, <5.64 GeV/c^ for S ^ A+//-X 

2.269< Mpx^ <2.302 GeV/c^ for A^ ^ A+tt" 
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Chapter 5 

Signal Yield in the Data 



In this chapter, we explain how the signal yield in the data is extracted. Ideally, if 
we were capable of fully reconstructing B hadrons in both hadronic and semilep- 
tonic modes, we would use the B hadron mass distribution to obtain the number 
of signal events (yield). However, a neutrino is missing in the semileptonic decay 
and the invariant mass of "charm+/x" is a broad spectrum with a shape which 
is poorly distinguished from the backgrounds. Therefore, a proper variable to 
use for the semileptonic mode is the mass of the charm hadron. We extract the 
yield by fitting the charm+vr (or charm) mass spectra in Figure 4.5 to a function 
which describes both the signal and the background. We integrate the signal 
function to obtain the yield. The signal function for all modes is a Gaussian or 
double-Gaussians. The background function varies with the decay mode. 

All our fits use an unbinned, extended likehhood technique. The general 
extended likelihood function (C) is expressed as: 

log£ = Elog{^^sig • Sim,) + iVbg-E(m,)} - N,,^ - N^^ + logC, (5.1) 

i 

where i represents i*'* event, m represents the reconstructed charm+7r (charm) 
mass. The amounts of signal and background are denoted as A^sig and A^bg, re- 
spectively, while S{m) {B{m)) are the functions which describe the signal (back- 
ground) mass spectrum. The last term in Equation 5.1, C, is a Gaussian con- 
straint on one fit parameter, x: 

C^g{x,^,a)^^e-'^^'^^\ (5.2) 
V27r(7 

where wc constrain the variable x around the mean fi. The difference of x and n 
follows a Gaussian distribution with an imcertainty a. The unbinned likelihood 
fitter calls the MINUIT package developed by James et al. [66]. MINUIT varies the 
fit parameters to minimizes —2 • log C . 

The performance of the fitter was checked on 1000 toy MC samples similar 
to the data distribution. We plot the pull distribution for each parameter, i.e. 
{x — Xq)/(Tx-i where x is the fit value, Xq is the generated (input) value, and (7-^ 
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is the uncertainty from the fit to the toy MC. For a large number of toy MC 
tests, the pull is expected to follow a Gaussian distribution. We examine if the 
fitter returns an output consistent with the input, i.e. if the mean of the pull 
distribution is consistent with zero and if the width is consistent with one. Note 
that the /i and a of the Gaussian constraint in Equation 5.2 are determined from 
a subsidiary measurement using the data and the MG. Therefore, we simulate 
this measurement in the toy MG test, by smearing the mean of the constraint 
with a Gaussian distribution of mean fi and sigma a in Equation 5.2. In order 
to evaluate the quality of the fit, we also superimpose the fit result on the data 
histograms and compute a x^. A complete description about the fitting and the 
pull distributions can be found in Yu [65] . Remark that as the B hadrons are fully 
reconstructed in the hadronic channels, the yields we extract are the true amount 
of signal for this analysis. The yields we extract for the inclusive semileptonic 
channels include the exclusive signals and indistinguishable backgrounds: such as 
muon fakes, decays from bb, cc, or other B hadrons. These backgrounds will be 
estimated in Ghapter 7 and subtracted in the calculation of the relative branching 
ratios. 



5.1 Mass Fit of the Semileptonic Modes 
5.1.1 dstarmu Yield 

As seen in Figure 4.5 (top right), the events with D*iJ in the final state have 
almost no combinatorial background. The combinatorial background is reduced 
largely by requiring M^n be consistent with the world average mass and 
cutting on the variable Mr,*^. We fit the mass difference Mdo^t — Moo instead of 
M^iOtt, because the width of Mpo^ — M^o is significantly narrower than that of 
MiyOj^. The signal to background ratio is thus higher in the signal region. The 
available Q of the D*^ decay is only about 7 MeV/c^, where Q is the momentum 
transferred to the daughters. After the Lorentz boost, carries most of Z}*+'s 
momentum and the bachelor pion from D*^ has a lower momentum. Therefore, 
the width of M^o^ is similar to that of M^o . While in the case of M^o^ — M^o , 
the Mjjo mass resolution is subtracted and only the momentum resolution of the 
soft TT will contribute to its width. 

The M^On — M^o distribution is fitted to a double Gaussian signal and a 
constant background. The extended log likelihood function is expressed as: 

l0g£ = ^loglA^sig- [(I-/2) •^l(mi,//,(7i) + /2-^2(mi,//, (72)] 

i 

+^bg • } - iVsig - A^bg, (5.3) 

-'"max -'"min 

where /2 is the fraction of the second Gaussian, The mass window 0.14 < Mdo^^ — 
Mdo < 0.18 GeV/c^ is specified by Mmin and Mmax- Both Gaussians have the 
same mean but different sigmas. Table 5.1 lists the mean, width of the pull 
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Table 5.1: D*ij, results from the unbinned likelihood fit. 



T 1 

Index 


Parameter 


1000 toy MC 


1000 toy MC 


Data nt value 








1 1 

pull mean 


11 • 1 j_1 

pull width 




1 


-^sig 




-0.023 ± 0.031 


1.006 ± 0.023 


1059 ± 33 





/2 




U.UU^ ± U.Uo4 


1 n7o -I- n no/i 

i.U I ^ ± U.UZ4 


U.OD ± U. iU 


3 


/i 


[GeV/c2] 


0.049 ± 0.033 


1.044 ± 0.024 


0.145410 ± 0.000016 


4 


CTl 


[GeV/c^] 


-0.048 ± 0.033 


1.052 ± 0.024 


0.00031 ± 0.00004 


5 




[GeV/c^] 


0.011 ± 0.032 


1.031 ± 0.023 


0.00071 ± 0.00006 


6 






0.010 ± 0.031 


1.000 ± 0.022 


321 ± 19 



distribution from 1000 toy MC test and the fit value of each parameter from the 
unbinned likelihood fit to the data. Figure 5.1 shows the fit result superimposed 
on the data histogram. We have obtained from the fit: 

^B^D'-^.-x = 1059 ± 33. 

5.1.2 DiJ, Yield 

A first glance of Mx-n-n in Figure 4.5 (middle right) might suggest that we could fit 
Mktttv to a Gaussian signal and a first-order polynomial background. But, since 
we do not apply particle identification (PID) in this analysis, the background 
under the signal contains not only the combinatorial background, but also con- 
tamination from the Ds decays. Not using PID means that a pion mass might 
be assigned to a kaon, and may be reconstructed as D+. Figure 5.2 shows 
the mis-reconstructed D'^ mass spectrum from the Bg — > D^ii~V^jt MC, where 
Ds are forced to decay into the final states listed in Table 5.2. These final states 
are selected after a study to identify the dominant Ds decays reconstructed in 
the D^ mass window. The MC used to assess Ds background is produced as 
described in Section 6.1. 

We need to include the mis-identified Dg mass shape in our hkelihood fit 
so to properly estimate the number of D^ events in the data. Assuming that 
B Dffi~X has a similar mass spectrum as Bs D^ ^~V^, we could use 
the MC for Figure 5.2 to obtain the function which describes the line-shape of 
mis-reconstructed Dg mass spectrum. We find the Dg spectrum {!F) could be 
described by a constant and a triangular function convoluted with a Gaussian 
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Table 5.2: Dominant mis-identified Dg sequential decays in Djj, signal. Branching 
fractions without uncertainties have an upper limit in the PDG. 



Selected final states of Dg decays 



Mode 



B (%) 



relative to B{Ds — > 07r) 



Dt 
Dt 

Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 
Dt 



(f)K+ 



/ovr+ 

p+r] 
p^r( 

k\+ 

K*^K+ 

TT+TT+TT" 

K+K+K- 



3.6 

0.03 

1.7 

3.9 

0.28 

0.04 

0.15 

0.57 

0.35 

10.8 

10.1 
0.4 
0.65 
3.6 
3.3 
0.005 
0.9 
0.02 



± 0.9 
±? 
± 0.5 
± 1.0 
± 0.11 
±? 
±? 

± 0.17 
± 0.12 
± 3.1 
± 2.8 
±? 

± 0.28 
± 1.1 
± 0.9 



± 



h0.022 



0.005 
± 0.4 

±? 



0.008 ± 
0.48 ± 
1.08 ± 
0.077 ± 
0.011 ± 
0.042 ± 
0.16 ± 
0.098 ± 



1 

0.05 
0.09 
0.025 



2.98 
2.78 
0.11 
0.18 
1.01 
0.92 
0.0014± 
0.25 ± 
0.0056± 



± 
± 
± 

± 
± 
± 



0.03 
0.022 
0.44 
0.41 

? 

0.06 

0.16 

0.09 

0.0007 

0.09 



? 
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Figure 5.1: M/jo^ — Mpo from the events fit to a double- Gaussian signal 
and a constant background. The result of the unbinned likelihood fit is projected 
on the histogram and a probability is calculated. 



in 

Tim) = (1 - U,) . ]_ + • Tim), (5.4) 

-'"max -'"min 

where /trg is the fraction of triangular function, Mmax and Mmin specify the mass 
window, 1.767 < Mk^^ < 1.977 GeV/c^, and 

, 2(m-Mo) ^ . X 

Here, ® represents convolution, ^ is the Gaussian and cxtrg is the width of Q. 
The triangular function value starts from zero at Mq and increases as the mass 
increases. When the mass reaches Mog, the function values is at its maximum 
and drops precipitously to zero. A graphical representation of Mos and Mq may 
be found in Figure 5.3. The exact form of Tim) is found in the appendix of 
Yu [65] derived by Heinrich. Figure 5.4 shows the result of the fit to the MC. 

Now with the function form of the Mk-k-k spectrum from the Ds decays, we 
have to normalize the MC yield to the data. The Dg yield may be obtained 
by reconstructing one of the Dg final states in the data: B Dffi^X, where 
07r+, — > K~^K~ , then using MC to determine the ratio of this Dg decay 
to that of all the decays in Table 5.2, R(j)n'- 

_ ^Jt^ 
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M(K7i7i) [GeV/cT 



Figure 5.2: Bg Dfn^Un MC reconstructed as Dfi final state. Here Dg are 
forced to decay into the modes listed in Table 5.2. The arrows indicate the 3 a 
signal region. 



The normalization of Dg is then expressed as: 

''B-^D+ti-X - ^ 



Ns^oU-. = '^"-"•'■'r-*"*^"" - (5.7) 

4>7T 



In order to obtain N-^^n+ii-x Ds^(Ptt (P^kk data, the same analysis cuts for 

Dfi are applied, except that we assign kaon mass to one of the same-sign charged 
tracks and pion mass to the other. We still assign kaon mass to the track which 
has the opposite charge of the other two. In addition, the candidates are required 
to pass the following cuts: 

• 1.767 < MKn^ < 1.977 GeV/c^ 

• \Mkk - 1.0191 < 0.01 GeV/c2 

The cut on Mkk guarantees that there is no mis-identified D~^ in the Df signal 
we reconstruct. We confirm this by reconstructing from the D^ ji~V 



MC and no Dg candidate is found. See Figure 5.5 for the B Dfn X signal in 
the data, we find: 

^B-*D+fj.-X,Ds^<t>7T,(f>^KK = 237 ± 17. 
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Figure 5.3: Graphical representation of the triangular function. 
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Figure 5.4: Mis-reconstructed Dg from the Bg —>■ ^~V^ MC fit to a constant 
and a triangular function convoluted with a Gaussian. The dashed curve indicates 
the result of the unbinned likelihood fit. 
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Figure 5.5: Data: Reconstructed B — > Dfjj,~X, where Df (t>T^^ ■, and 
K^K~ . Mktttt is required to be between 1.767 and 1.977 GeV/c^. There are 237 
±17 events in the peak. 



We then reconstruct the same Dg decay chain in the MC as in the data and 
obtain 

= 0.131 ± 0.007. (5.8) 

Inserting the result of N^^j~,+ ^--^ j-,^^^^ ^^^^^ and R^.„ into Equation 5.7, we 
have: 

NB-.Dt,-x = 1812 ± 160, (5.9) 
The uncertainty in Equation 5.9 comes from the fractional uncertainties on: 

NB^Dt,-x,D.^^.,^^KK (7.2%) and R^^ (5%). 

In the unbinned fit, the extended log likelihood function is expressed by the 
sum of two likelihoods: one describing the data and the other describing the 
Bg — s> D'^fi~T7^ MC since we fit the data and MC simultaneously; 

log/: = iog/:<^^*^ + logc^^, (5.10) 

The likelihood function for the data, logC'^"'^"', is a sum of a signal Gaussian, a 
first-order polynomial for the combinatorial background (7i), and the function 
for the Ds {J-', see Equation 5.4). A Gaussian constraint on the amount of Dg, 
, is employed. 

log C"^'^ = Yl logl^'sig ■ Girrii, /i, a) + iV,o^bg ■ ^(m^) + Nd, ■ ^(m,)} 

i 

- N,,g- N,omhg- Na^+logCD^, (5.11) 
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Table 5.3: D/j, results from the unbinned likelihood fit. 



Index 


Parameter 


1000 toy MC 
pull mean 


1000 toy MC 
pull width 


Data fit value 


1 


N ■ 


-0.012 ± 0.035 


1.004 ± 0.025 


4720 ±100 


2 


H [GeV/c2] 


0.027 ± 0.037 


1.048 ± 0.027 


1.8680 ± 0.0002 


3 


a [GeV/c2] 


0.007 ± 0.035 


0.992 ± 0.025 


0.0084 ± 0.0002 


4 


-^combg 


-0.076 ± 0.038 


1.073 ± 0.027 


15178 ± 197 


5 


Pi 


0.018 ± 0.036 


1.027 ± 0.026 


-5.2 ± 0.7 


6 


Nd. 


0.042 ± 0.037 


1.065 ± 0.027 


1832 ± 155 


7 


/trg 


0.022 ± 0.036 


1.023 ± 0.026 


0.617 ± 0.021 


8 


Mo [GeV/c2] 


0.055 ± 0.035 


1.007 ± 0.025 


1.69 ± 0.02 


9 


Moff [GeV/c2] 


-0.025 ± 0.036 


1.019 ± 0.026 


1.888 ± 0.002 


10 


^trg [GeV/c2] 


-0.035 ± 0.037 


1.056 ± 0.027 


0.010 ± 0.002 



where 

nj I \ 1 I / -^max ~l~ -^min \ 

^(^^) = M -M ■ +^^-(^^ ^ )' 

-"^max -"'-'mm ^ 

From the prediction of Equation 5.9, we have ijP = 1812, and = 160. 

The hkelihood function log>C'^*-' is used to fit Bs — > D'^ii'V^ MC and obtain 
the parameterization of !F{m). Here the normalization does not matter. 

log£^c^^log{^K)}. (5.12) 

i 

Table 5.3 hsts the mean, width of the pulls from 1000 toy MC test and the result 
returned from the unbinned likelihood fit to the data. Figure 5.6 shows the fit 
result superimposed on the data histogram. We have obtained from the fit: 

Nb^d^,-x = 4720 ± 100. 

We also perform a cross-check by removing the constraint on iV^o^ and obtain 
N-B^D+^-x = 4667 ± 139, Nd^ = 2184 ±620, which are consistent with the result 
of the constrained fit. The fit without constraint has a x^/NDF=197.0/199 and 
probability is 52.7%. 

5.1.3 AcM Yield 

When a proton mass is assigned to a kaon or pion, numerous B meson to D meson 
semileptonic decays could be mis-reconstructed as a Ac/x final state. In order to 
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Figure 5.6: Mk-k-k from the Dix events fit to a Gaussian (signal), a first-order 
polynomial (combinatorial background), and a constant plus a triangular function 
convoluted with a Gaussian (marked by the dashed line). The result of the 
unbinned likelihood fit is projected on the histogram and a probability is 
calculated. 



estimate the B meson background shape under our signal, we use generator level 
MC and generate the semileptonic decays (/i channel) of each B meson fiavor 
separately. After applying analysis cuts, we add up the mis-reconstructed mass 
spectrum from each kind of B meson according to the production fractions: 

b^Bd = (39.7 ± 1.3) %, 
Bu = (39.7 ± 1.3) %, 
Bs = (10.7 ± 1.1) %. 

Figure 5.7 shows a smooth mass spectrum from the generator MC. The shape is 
best described by a second-order polynomial, with x^/NDF = 36.6/42, prob = 
70%. A first-order polynomial fit yields x^/NDF = 56.6/43, prob = 8%. Because 
the combinatorial background may be parameterized by a first-order polynomial, 
and adding a first- to a second-order polynomial gives a second-order polynomial, 
we fit the combinatorial and the B meson background together to a second-order 
polynomial (H). The extended log likelihood function could be expressed as: 

log£ = ^ log{Ar,ig ■ g{m^, /i, a) + A^bg ■ ^(m^)} - A^^ig - A^bg, (5.13) 

i 

where 

= Xf M + ^1 ■ ^""^ ~ ^'^''^'> + ^2 ■ (12 ■ (mi - M„,id)' - M^g). 

-'^max -'"min 
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Here, Mmax and Mmin specify the A+ mass window: 2.19 < Mp^-n- < 2.37 GeV/c^. 
The average of Mmax and Afmin, or the mid point in the mass window is Mmid- The 
difference of M^ax and Mmin is Mjig. Table 5.4 hsts the mean, width of the pulls 
from the toy MC test and the parameter value from the fit to the data. Figure 5.8 
shows fit result superimposed on the data histogram. We have obtained from the 
fit: 

NB-.At,-x = 1237 ± 97. 



Table 5.4: Ac/x results from the unbinned hkelihood fit. 



Index 


Parameter 


1000 toy MC 
pull mean 


1000 toy MC 
pull width 


Data fit value 


1 




0.018 ± 0.030 


0.997 ± 0.022 


1237 ±97 


2 


/i [GeV/c2] 


0.017 ± 0.033 


1.070 ± 0.024 


2.2850 ± 0.0005 


3 


a [GeV/c2] 


-0.069 ± 0.032 


1.036 ± 0.023 


0.0074 ± 0.0006 


4 




0.004 ± 0.031 


1.021 ± 0.022 


16576 ± 157 


5 


Pi 


0.010 ± 0.031 


1.007 ± 0.022 


-4.3 ± 0.8 


6 


P2 


0.020 ± 0.031 


1.012 ± 0.022 


-3.7 ± 1.8 
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Figure 5.7: B ^ X^v^ MC mis-reconstructed as Ac/i. The mass spectrum fit 
to a second-order polynomial. x^/NDF = 36.6/42, prob= 70.7%. Two arrows 
indicate the signal region of A+. 
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Figure 5.8: MpXn from the Ac/i events fit to a Gaussian (signal), a second-order 
polynomial (combinatorial + B meson). The result of the unbinned likelihood fit 
is projected on the histogram and a probability is calculated. 
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5.2 Mass Fit of the Hadronic Modes 



Figure 4.5 shows that to the left of the hadronic signal peak, the "charm+7r" mass 
spectrum exhibits an interesting structure. In order to extract a correct number 
of events observed in the hadronic channels, we have to take into account the 
background structure when fitting the charm+vr mass spectrum. For the B 
D~^7r~ and A;, — > A^tt" modes, we import the and A5 mass functions derived in 
the analyses of Furic [67], and Maksimovic [68], respectively. Several parameters 
that describe the background shapes or normalizations are fixed. We find small 
modifications are needed for the numerical values of the fixed parameters in the 
B — > D~^7i~ mode, as a few variables we apply cut on are different from those 
in Furic's analysis. Wc apply our cuts on the MC used in Furic's analysis and 
refit the MC to extract the numbers for our analysis. For the B — * D*^tt~ 
mode, we produce an inclusive B — > D*^X MC sample to study the background 
composition. The decay modes with distinguished mass shape are separated from 
the other modes. The decays with similar mass spectra are lumped together and 
fit to the same background function. Figure 5.9 shows the B^ and A;, mass 
spectra from the contributions of different decays. 

Our hadronic mass spectra share several common features: It is clear that the 
background from the B hadron decays only contribute to the mass region below 
the signal, while in the data, the upper mass region is composed of combinato- 
rial background, which may be described by an exponential or a constant. The 
combinatorial background extends down to the lower B mass region as well. In 
the region 40 to 70 MeV/c^ below the signal peak, Cabibbo suppressed decays, 
S° D*+K-, B° D+K^, Afe A+is:-, with a branching ratio about 8% 
of our Cabibbo favored signals, produce a small contamination. Going further 
down in the charm+vr mass, we have partially reconstructed B decays from the 
semileptonic modes, and other mis-identified B hadronic decays. 

Note that since both B D*~^7r~ and Af, — > A+tt" have low statistics, we 
constrain the widths of their signal Gaussians in the following way: We first fit the 
width of Mdtt {o-'^^^) from the high statistics B — > D+tt" sample (~600 events) 
in the data. Then we multiply cr^^^^ with the MC width ratio: crf^^^.^/a^^ and 
predict ai^^^jj.^. 

5.2.1 D*+7r- Yield 

The study from the B — > D*'^X MC shows that the background in the lower 
mass region is dominated by the following decays: Cabibbo suppressed decay 
B° D*+K-, 5° D*+p-, and the remaining B D*+X. See the texts 
below for the detailed descriptions. 

1. 5° — > D*^K~: fully reconstructed Cabibbo suppressed decays. The mass 
spectrum is a peak about 40 MeV/c^ below the B D*^tt~ signal, with 
small tails on the lower mass side. The shape is modeled by a lifetime 
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Other B decays 
Other decays 
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Figure 5.9: Background composition for B D*^-n^ (top), B D^n^ (mid- 
dle) from Furic [67] and A;, A^tt^ (bottom) from Maksimovic [68]. Note that 
the mis-reconstructed Bg -D+tt^ and Af, — A+tt" in the Dn mass is not shown. 
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function; 

S{m) = exp(m, rD*K) <S> G{m, ij,d*k, (^d*k), (5-14) 
where td*k is the hfetime, fiD*K is the zero point of the hfetime function also 
the mean of the Gaussian. The width of the Gaussian also the resolution 
of the lifetime function is (T£,*k- The exact form of S{m) is found in the 
appendix of Yu [65]. See Figure 5.10 (top) for the fit to E° ^ D*+K- MC. 

2. ^ D*~^p~ , where p~ — > tt'^tt": modeled by a triangular function convo- 
luted with a Gaussian; 

See Figure 5.10 (middle) for the fit. 

3. Continuum: remaining B — > D*'^X decays partially reconstructed. These 
backgrounds have similar mass spectrum and are group together. The shape 
is modeled by a first-order polynomial with a negative slope and a turn-off 
at M°^^^'^; when m < M°^^^'^: 

'Him) = ,,,,, • (ikri'^"'"^ - m), (5.16) 

V ) (M°ti^-B _ ^^.^2 V off \ ) 

and when m > M°^^^'^: 

n{m) = 0. (5.17) 

The lowest boundary of the L>*7r mass window, M-^An-, is 4.6 GeV/c^. See 
Figure 5.10 (bottom) for the fit to these MC samples. 

In the unbinned fit, the extended log likelihood function is expressed by the 
sum of five likelihoods: one describing the data, and the other four describing 
the MC samples from each type of background and the signal: 

\ogC = log£^^*^ + log£^?^ + log£^F^ + log£^S, + logC^B, (5-18) 

The likelihood function log£'^'^*'^ is a sum of a signal Gaussian, a constant com- 
binatorial background, the functions for D*K (S), D*p (T), and the continuum 
(H). In addition, there is a constraint on each of the following parameters: the 
signal width, relative amount of D*K to the signal {fD*K), and the fraction of 
D*p in D*p -\- remaining B D*'^X {fD*p) The reason for the last constraint 
is because B D*~^p~ and the remaining B D*~^X decays occupy the same 
mass region. Therefore, the likelihood fit converges faster if we constrain fo'p- 



+iVbg ■ [/c 



combg 



+ (1 - /combg) • [fo'p • T'irrii) + {1- fo'p) ■ ^(m^)]} 
- A^sig - iVD*K - iVbg 

+ logCi + logC2+logC (5.19) 
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Figure 5.10: Various MC samples reconstructed as 5 ^ D*~^n . From the top 
to the bottom are W D*+K', B D*+p-, and the remaining; B D*~^X. 
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Table 5.5: Branching ratios and relative efficiencies for B — > D*~^7r background. 







> D*+K- 




B(%) 


0.276 ± 0.021 


0.020 ± 0.005 


e ratio 




1 


1.02 ± 0.02 


fD*K 


0.071 ± 0.019 






D*+p- 


remaining B D*^X 


Nmc 




758 


2371 


fD*p 


0.242 ± 0.008 



where £{mi), T{mi) and ^.{mi) are expressed in Equations 5.14-5.16. The Mmax 
and Mmin specify the mass window: 4.6 < Afo*^ < 5.6 GeV/c^. The parameters 
fD*K, -/Vbg, /combg and /dv are defined as follow: 



D*K 



K 



N ■ 



A^bg = Ncombg + -^otherB + Nd*p-i 
-^combg 



f combg — 



bg 



-^otherB + 

The constraints are expressed as: 

Cl — G{fD*K, 

Ca = S{a, iJ,p,ap), 

where = 0.071, ai = 0.019, /^a = 0.242, (T2 = 0.008, Hp = 0.0259 GeV/c^, and 
ap = 0.0012 GcV/c2. 

Here, /Xp and are determined using the 5° — > D'^Tr~ signal in the data, 
5° — > D'^TT" and 5° — > D*+7r" MC as described earlier. The /ii and (7i are 
determined using the world average branching ratios, and the efficiencies from 
the MC listed in Table 5.5: 

The /i2 and ca are determined by counting the number of reconstructed D*p and 
the remaining B — > D*^X events in the MC after all the analysis cuts. 
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The three hkeUhoods for the background MC are used to obtain the parame- 
terization of S{'m), T{m), and 7i(m). The normahzations do not matter here. 

Elog^K), (5.21) 

i 

Elogr(m,), (5.22) 

i 

^logUimi). (5.23) 

i 

In addition, log jC^^^ is used to obtain the reconstructed mass difference between 
MC and data, maifr- In the log£^^*^, all the parameters except the normahzation 
and the resolution parameters (a) for the signal Gaussian and the background 
functions, differ by mdifr from those in the logC^^'^. The resolutions for all the 
backgrounds arc kept the same between MC and data, while the resolution of the 
signal Gaussian in the data is a separate free parameter from that in the MC. 

We use the total likehhood to fit the data and MC simultaneously. Table 5.6 
lists the pull means and widths of toy MC test and the unbinned likelihood fit 
result to the data. Figure 5.11 shows the fit result superimposed on the data 
histogram. We have obtained from the fit: 

If we remove the constraint on the signal width, we find N-^ = 110 ± 11 

and (Jdata ~ 0.0295 ± 0.0033GeV/c . Removing the constraint on /d^k gives us 
^B°^D*+i:- ~ ±11 and fD*K — 0.053 ± 0.053. Removing the constraint on 
fD*p gives us Ns°-^D*+Tr- ~ ± 11 and fo^p — 0.38 ± 0.07. In conclusion, the 
un-constrained fits return a value consistent with the constrained fit, but with 
larger uncertainties. The fit xV^DF are 20.0/12, 20.8/12, 16.5/12 and the fit 
probabilities are 6.7%, 5.4 %, 16.9% for the three different unconstrained fits. 



loff f.^^ 

^'-'6 ''-otherB 



5.2.2 W D+n- Yield 

As noted earlier, we make use of the mass function derived in Furic's analysis 
for the B — > D^t:~ mode. The parameters which are kept constant in Furic's 
mass function remain constant in our analysis. The following backgrounds con- 
tribute to the mass spectrum of from Furic's study: Cabibbo suppressed decay 
B^ D~^K~ , D*+7r~, S° D~^P~, remaining B D^X and the combi- 

natorial background. Recent study by Bclloni, Martin and Piedra et al. [68] [69] 
shows that the mis-reconstructed Bg -D^tt" and A;, — > A+tt" produce small 
contamination in the — > D^t:~ signal. See the text below for the detailed 
descriptions. 

1. B^ ^ D^K~: fully reconstructed Cabibbo suppressed decays. The mass 
spectrum is a peak about 60 MeV/c^ below the B — > D^n" signal. The 
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Table 5.6: B — > D*'^7r results from the unbinned likelihood fit. 



Index 


Parameter 


1000 toy MC 


1000 toy MC 




Data fit value 










pull 


mean 




pull width 






1 


N ■ 




-0 


019 ± 


0.031 





964 


± 


0.022 




106 ± 11 


2 




[GeV/c2] 





013 ± 


0.033 


1 


032 


± 


0.024 


5 


2772 ± 0.0002 


3 


0"MC 


[GeV/c^] 


-0 


043 ± 


0.033 


1 


036 


± 


0.024 





0262 ± 0.0002 


4 


fD*P 




-0 


009 ± 


0.032 


1 


006 


± 


0.023 




0.244 ± 0.008 


5 




[GeV/c2] 





027 ± 


0.033 





943 


± 


0.020 




4.43 ± 0.01 


6 




[GeV/c^] 





045 ± 


0.032 





993 


± 


0.023 




5.134 ± 0.001 


7 




[GeV/c2] 


-0 


077 ± 


0.033 


1 


001 


± 


0.024 




0.026 ± 0.001 


8 


/combg 




-0 


045 ± 


0.032 





940 


± 


0.023 




0.09 ± 0.03 


9 






-0 


010 ± 


0.033 


1 


019 


± 


0.023 




428 ± 21 


10 




[GeV/c2] 


-0 


220 ± 


0.031 





972 


± 


0.022 




5.174 ± 0.004 


11 






-0 


047 ± 


0.033 


1 


016 


± 


0.023 




0.069 ± 0.018 


12 


IJ'D*K 


[GeV/c2] 


-0 


032 ± 


0.033 


1 


024 


± 


0.024 


5 


2345 ± 0.0009 


13 


td*k 


[GeV/c2]-i 


-0 


017 ± 


0.032 





986 


± 


0.023 





0287 ± 0.0009 


14 


aD*K 


[GeV/c2] 





029 ± 


0.033 


1 


029 


± 


0.024 





0254 ± 0.0006 


15 




[GeV/c2] 


-0 


034 ± 


0.032 





992 


± 


0.023 




0.005 ± 0.003 


16 




[GeV/e^] 


-0 


050 ± 


0.031 





971 


± 


0.022 




0.026 ± 0.001 



63 




Figure 5.11: M/j.^r from the B D*~^7r~ events is fit to a Gaussian (signal), 
a constant (combinatorial), and the background functions for the lower mass 
spectrum as described in the text. The result of the unbinned likelihood fit is 
projected on the histogram and a probability is calculated. Note that the bins 
with less than 20 entries are combined. 



shape is modeled by a single Gaussian; 

P/C(m) = g{m, /i - AMdk, ctdx), (5.24) 

where the shift of Gaussian mean from the b'^ D^n^ signal, AMjjx, 
and the width, ajjK, are extracted from the MC. 

2. Bg ^ DfTT~, where —>■ (f)TT~^ and (p ^ K^K~: this decay produces a 
peak at around 5.31 GeV/c^ when the pion mass is assigned to one of the 
kaons. The spectrum is modeled by double Gaussians with the same mean; 

Bs{m) = fi ■ g{m,ij,B,,(Ti) + (1 - /i) ■ g{m,ij,B,,a2), (5.25) 

where the fraction /i, fis^, cri and (72 of each Gaussian are obtained from 
the fit to the MC as shown in Figure 5.12 (top left). 

3. Afe A+7r~, where — > pK^n^: this background produces a broad peak 
around 5.4 GeV/c^, the region where the pion mass is mis-assigned to the 
proton. The spectrum is modeled by a lifetime function; 

Ci^{m) = exp(m, taJ (g) g{m, /ia^, (TaJ, (5.26) 



64 



where //a,, and are the zero point and the resolution of the hfetime 
function. See Figure 5.12 (top right) for the fit to the A+tt" MC 

when reconstructed as D+tt". 

5° D~^P~, where p~ 7r°7r~ and D*'^'k~ where D*+ D+7r°: 

These two backgrounds are combined. The spectrum of B ^ D~^p~ looks 
like B D*~^p~ in Figure 5.10 (middle) and is modeled by a lifetime 
function. The spectrum of 5^ D*^ti~ is composed of two horns and is 
modeled by two Gaussians with different means. 

The structure of double horns arises for the following reasons: When B^ — > 
D*+7r'~, — > D+7r°, is reconstructed as D'^7r~, the mass is lower than 
the world average S° mass due to the missing 7r°. The amount of the 
negative mass shift, AM, is determined by the angle between the ir^ and 
the D*~^ flight direction, d(f). Because both B^ and vr" are scalars (spin=0), 
to conserve the total angular momentum in the decay, the vector particle 
(spin=l), D*^, is transversely polarized. The angle dc/) from a transversely 
polarized D*'^ is cos^ 9 distributed and the most probable dcf) is either or 
180 degrees. Therefore, AM is quantized and this forms a double-horns 
spectrum. 

After combing 5° — > D'^p~ and B^ D*~^'k~, we have: 

7^(m) = (1 - fn) ■ exp(m, r^f) Q{m, //ref, CTref) 

+ fn • (0.5 • G{m, jl^ei - i^ref - (^ref , CTr) 

+0.5 • Q{m, Prei - J^ref + (^ref, CTh))- (5.27) 

The exact form of the lifetime function is found in the appendix of Yu [65]. 
The zero point of the lifetime function is //ref and z/ref is the offset of the 
mid point between two horns from the lifetime function. The //ref and i^ref 
are left free in the likelihood fit to the data. The values of the following 
parameters are extracted from the fit to the MC, as shown in Figure 5.12 
(bottom left), and kept constant in the fit to the data: the lifetime (rref), 
the fraction of horns {fn) , the half distance between the peak of two horns 
((^ref)) the resolution of the lifetime function (a"ref) and the width of both 
horns {(Th)- 

Continuum: remaining B — > D^X decays and partially reconstructed. 
These backgrounds have similar mass spectrum and are group together. 
The shape is modeled by a first-order polynomial with a negative slope and 
a turn-off at Mog; when m < Mog: 

= (M ^2 ■ (^off - (5.28) 

(Moflf - M^i^y 

and when m > M^s: 

n{m) = 0. (5.29) 
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The lowest boundary of the Dtt mass window, Mmin, is 4.6 GeV/c^. See 
Figure 5.12 (bottom right) for the fit to these MC samples from Puric's 
analysis [67]. 

6. combinatorial: modeled by an exponential function. When the slope of the 
exponential, po, is not zero, 

Sxpim) = Po Tj T:r- ■ e-^'^-^^-^-''^), (5.30) 



and when po is zero. 



^-^("^^ = M -M. ' ^'-^'^ 



where Mmax and Mmin specify the mass window: 4.6 < Afojr < 5.6 GeV/c^ 
and Mmid is the average of Mmax and Mmin- 

In the unbinned fit, the extended log likelihood function is expressed as a sum 
of a signal Gaussian, the functions for the DK mode ("P/C), Dgii {Bs), Acvr (£/?), 
D*7i plus Dp {TV), the remaining B D^X decays (7i), and the combinatorial 
background {£xv)'- 

log £ = ^ log{A^sig • [Q{mi, //, a) + foK ■ T>JC{mi) 

i 

+fBs ■ Bsijrii) + /a, • CBirrii)] 

+7Vbg ■ [(1 — /combg) ■ [(1 — /otherB) " TlilTli) + /otherB " 7Y(mj)] 
+/combg • ^Xvi^h)]} 

- A^sig-(1 + /DK + /B.+/Aj-A^bg, (5.32) 

where VKirrii), Bs{mi), Cts^rrii), IZirrii), 7{{mi) and Exvij^i) are expressed in 
Equations 5.24-5.31. The fractions foKi fss and /a^ are the ratios of N^k, 
and to the signal, iVgig. The total amount of combinatorial background, the 
backgrounds from the Dp, D*7r, and the remaining B decays is denoted as A^bg- 
The parameters /combg and /otherB are defined as follows: 

r A^combg 

/combg = ^ ) 
-'^bg 

r _ -^otherB 

/otherB = 



A^otherB + A^D*,r + Nop 



All the fractions and ratios here except /combg are kept constant in the likelihood 
fit. The — > D'^K~ fraction, fox, is determined from the world average 
branching ratios; 
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Table 5.8 lists the values of the branching ratios in Equation 5.33. We have 
fDK = 0.073 ±0.023. 

The Bg fraction, Jb^, is obtained using the formula: 

fs B(B,^D-n-) B{D- -ron+)B{c:>-rK+K-) 



-r(Dt - HK^K-)^-) • ' 

where the branching ratios are from the 2004 PDG and the CDF II measurement 
jr^ ■ ^|^u^^-^^_| by Furic [67]. The efficiency ratio is obtained by applying our 

Dtt analysis cuts on the Bg MC. Inserting the numbers listed in Table 5.8 into 
Equation 5.34, we obtain /bs — 0.006 ± 0.001. Note that the uncertainties from 
the branching ratios of 0, Dg, and D decays vanish after multiplying Furic's result 

with the ratio: V^-:^!^";)""^ 

The Ab fraction, /a^, is obtained using a similar formula; 

_ aA,(PT>6.0)g(A,^A+7r-) g(A+^pA-7r+) <^a+.- 

anoiPr > 6.0)B(b' - D^n') B{D^ - X-7r+7r+) eM,^_^^^^_ ' ^'■"''^ 

where the product of the first and the second terms come from 2004 PDG and 
CDF II measurements by Le, et al. [68]. The uncertainties from the branching 
ratios of Ac and D decays vanish in Equation 5.35. The efficiency ratio is obtained 
using the A^ A+tt" MC. The value of fA, is then 0.031 ± 0.005. Table 5.8 lists 
the numerical values of Le's result and the MC efficiency. Finally, /otherB is 
obtained using Furic's B D~^X MC. We apply our analysis cuts and count the 
number of D*7t + Dp and the remaining B —>■ D^X events. We find /others = 
0.569 ±0.011. 

Table 5.7 lists the constant parameters with their values and uncertainties 
obtained from the fit to the MC. Table 5.9 lists the mean, width of the pulls from 

the toy MC test and the value of each fit parameter from the fit to the data. 
Figure 5.13 shows the fit result superimposed on the data histogram. We have 
obtained from the fit: 

ALo ^ , = 579 ± 30. 



67 



Table 5.7: Fixed parameters in the B — > D+tt unbinned likelihood fit. 



Parameter 


Meaning 


Value 


foK 


^B°^D+ K- 1 

mass shift of 5° ^ D+iT" [GeV/c^] 


0.073 


±0.023 


AMdk 


0.067 


±0.006 




width of 5° ^ D+i^-[GeV/c2] 


0.032 


±0.009 


fBs 




0.006 


±0.001 




mean of Bs background [GeV/c^] 


5.307 


±0.001 


/l 


fraction of the narrow Bs Gaussian 


0.773 


±0.002 


(Tl 


width of the narrow Bg Gaussian [GeV/c^] 


0.021 


±0.002 


<y 2/(^1 


width ratio of the Bg Gaussians 


1.8 


±0.3 






0.031 


±0.005 




mean of A;, [GcV/c^] 


5.416 


±0.002 




width of A5 background [GeV/c^] 


0.024 


±0.002 


^A, 


Ufetime of background [GeV/c^~^] 


0.052 


±0.002 


Tref 


lifetime of Dp background [GeV/c^~^] 


0.36 


±0.06 


O'rcf 


width of Dp background [GeV/c^] 


0.039 


± 0.008 


/// 


fraction of D*7r horns 


0.20 


±0.06 


(^ref 


distance between two horns [GeV/c^] 


0.039 


± 0.003 


(7h 


width of the horns [GeV/c^] 


0.019 


± 0.003 


/otherB 


fraction of the remaining B — > D^X 


0.569 


±0.011 


Moff 


cut off for B D+X mass [GeV/c^] 


5.112 


± 0.007 



68 



Table 5.8: Parameter values used to determine /dk, fs^ ^^^^ /aj- 



B(B^ D+K-) 


(2.0 ± 0.6) xlO-^ 




(2.76 ± 0.25) X 10-^ 


foK 


0.073 ± 0.023 


f, B(Bs^Dt-n-) 


0.35 ± 0.05(sto^) ± 0.02 {syst) ± 0.09 {BR) 


B{Dt ^ 0vr+) 


(3.6 ± 0.9)% 


^ K+K-) 


(49.1 ± 0.6)% 


Y<( /-)+ V T(^ — ^—\ 

L li-Vg — rl\ ' J\ tT ) 

r{D+-*<t>(K+K-)7T-) 
MC UMC 
^B,-^D.+7r-/ B°-»D+7r- 


(9.2 ± 0.6)% 
0.81 ± 0.08 
0.071 ± 0.004 




0.006 ± 0.001 


aAj^(PT>6.0)B(Af,^A+7r-) 


0.82 ± 0.08 {stat) ± 0.11 {syst) ± 0.22 {BR) 

(5.0 ± 1.3)% 


(^j30(^'T>6.0)i3("B"^D+7r-) 

B(A+ ^ pX-7r+) 


^A.^A+TT-/ V^n+^- 


0.069 ± 0.002 


/a6 


0.031 ± 0.005 



Table 5.9: B — > D+tt results from the unbinned likelihood fit. 



Index 


Parameter 


1000 toy MC 


1000 toy MC 


Data fit value 








pull mean 


pull width 




1 


N ■ 




0.012 dz 0.035 


1.021 ± 0.026 


579 ±30 


2 


/X 


[GeV/c2] 


-0.026 ± 0.034 


0.989 ± 0.025 


5.278 ± 0.001 


3 


a 


[GeV/c2] 


-0.040 dz 0.035 


1.015 ± 0.026 


0.0235 ± 0.0012 


4 






0.017 ± 0.034 


0.990 ± 0.025 


4049 ± 67 


5 


A*ref 


[GeV/c2] 


-0.036 ± 0.036 


1.037 ± 0.026 


5.145 ± 0.015 


6 




[GeV/c2] 


-0.037 ± 0.038 


1.085 ± 0.028 


0.068 ± 0.020 


7 


/combg 




-0.145 ± 0.034 


0.988 ± 0.025 


0.583 ± 0.044 


8 


Po 




-0.051 ± 0.034 


0.976 ± 0.024 


1.75 ± 0.15 
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Figure 5.12: Various MC samples reconstructed as 5 D^n . From the 
top left to the bottom right are 

B^ D^p^, and the remaining B D^X. The fit probabilities are 32.5%, 
66.8%, 16.5% and 32.9 %. 
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M(D7r) [GeV/c^] 

Figure 5.13: iW^o^ from the D^n~ events is fit to a Gaussian (signal), an 

exponential (combinatorial), and the background functions for the lower mass 
spectrum as described in the text. The result of the unbinned likelihood fit is 
projected on the histogram and a probability is calculated. 



5.2.3 Ab A+TT- Yield 

We use the mass function derived in the analysis of Maksimovic for the fit to A;, 
A^vr^ data [68]. We cross-check the values of the background shape parameters 
by applying our analysis cuts on the MC used in the Maksimovic analysis. We 
find the same numbers can be used for this analysis. The following backgrounds 
contribute to the mass spectrum of Af, from their study: Cabibbo decay A^ 
A+JC^, four-prong mis-identified B meson, the remaining B meson decays, the 
remaining Af, decays and the combinatorial background. 

1. Afe ^ A+i^": fully reconstructed Cabibbo suppressed decays. The shape 
is modeled by two Gaussians of different mean and width; 

CcfCim) = fi ■ g{m, (Ti) + (1 - /i) ■ Gim, /i^^^, aa), (5.36) 

where /i, fi\^K^ "^i; f^A^K '^2 are from the fit to the MC. 

2. mis-identified four-prong B mesons: all the B mesons with four tracks in 
the final states and fully reconstructed. B —>■ D^n^ contributes about 
50% of this type of background. Since these decays have similar final state 
as our Ab — > A+vr" signal, they produce a distinguished peak to the left of 
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the signal Gaussian. This background (Bi-pnoMg) is modeled by the sum 
of a Landau (>C^x>) and a Gaussian function: 

B4rnoMg{m) = /l • CM/v{m, hbpl, ctl) + (1 - /l) • Q{m, /^bpg, ctg), (5.37) 

where /l, /iBPL and (Tl are the fraction, mean and the width of Landau 
distribution. The mean and the width of the Gaussian are denoted as /xbpg 
and (Tg- These parameters are extracted from fit to the MC as shown in 
Figure 5.15 (bottom). 

3. remaining B meson decays: this background {Ob) spectrum is modeled by 
the sum of an exponential function and a product of a bifurcated Gaussian 
(Bj^) with a step-down function: 

Oeim) = £xv{m)+fuig-Br{m, /^ob, ^ob^ ^^b)-(l- ^ ^ g(M„i-m)/agb ) ^^-^^^ 

where 8xv{i^) is expressed in Equations 5.30-5.31. The parameters /bifg, 
A*ob) (^ohi '^ob fraction, mean, left sigma, right sigma of bifurcated 

Gaussian. The step-down function parameters, Hohst and ag^, together with 
the parameters for the bifurcated Gaussian, are extracted from the MC as 
shown in Figure 5.15 (top left). The exact form of the bifurcated Gaussian 
is found in the appendix of Yu [65] . 

4. remaining Afc decays: this background {Oc) spectrum is modeled by the 
sum of two Gaussians and the product of a bifurcated Gaussian and a 
step-down function 

+ -B^(m,;.oi,a,^.0-(l- ^^^(,i_^VaJ (5.39) 

where the parameters in the function are from the fit to the MC as shown 
in Figure 5.15 (top right). 

5. combinatorial background: described by an exponential function 

In the unbinned fit, the extended log likelihood function is expressed as a 
sum of a signal Gaussian and the functions for K^K (CqIC), four-prong B meson 
{B4pTioj\fg) , remaining B meson decays (Cb), remaining decays (Oc) and the 
combinatorial background {Sxv)- In addition, there is a constraint on the width 
of the signal Gaussian determined using the — > D'^Tr~ data, — > D^t:~ and 
Afe A+7r~ MC as described earlier. 

\ogC = log{^sig ■ [G{mi, fj,, a) + /a,x ■ Cc}C{mi)] + A^B4prong ■ B4vnojvg{mi) 

i 

+NoB ■ OB{mi) + A^oL • Oc{mi) + N^^^^^ ■ £xv{mi)} 

- AT^ig • (1 + J^^k) - A^B4prong " A^OB - A^OL - A^combg + log C (5.40) 
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where CclCirrii), ^ivnoMgij^i)) Osij^i)) Ocijrii) and Exvij^i) are expressed in 
Equations 5.36-5.39 and Equations 5.30 -5.31. The fraction fx^K is defined as: 



and is fixed to 0.08; the number is suggested by the branching ratio of the 
Cabibbo suppressed relative to the Cabibbo favored decay in the B meson sys- 
tem. The mean and the sigma of the Gaussian constraint for the signal width 
are 0.0231 GeV/c^, and 0.0012 GeV/c^, respectively. 

Table 5.10 lists the mean, width of the pulls from the toy MC test and the 
fit result to the data. Table 5.11 lists the values of the constant parameters. 
Figure 5.14 shows the fit result superimposed on the data histogram. We have 
obtained from the fit: 

iVA,^A+.- = 179 ± 19. 

We also cross-check by removing the constraint on the signal width and obtain 
AT + = 177 ± 2: 

result in Table 5.10. T 
fit probability of 20.2 %. 



^Ab^A+TT- = 177 ± 22, and a = 0.022 ± 0.004, which are consistent with the 
result in Table 5.10. The fit without constraint has a x^/NDF of 123.2/111 and 




Figure 5.14: Ma^tt from the — * A+vr^ events is fit to a Gaussian (signal), an 
exponential (combinatorial), and the background functions for the lower mass 
spectrum as described in the text. The result of the unbinned likelihood fit is 
projected on the histogram and a probability is calculated. 
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Figure 5.15: Various MC samples reconstructed as —>■ A^tt^. From the top 
left to the bottom are the remaining B meson decays, the remaining A^ decays 
and four-prong B meson decays. 
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Table 5.10: — > A+tt results from the unbinned likelihood fit. 



Index 


Parameter 


1000 toy MC 


1000 toy MC 


Data fit value 








pull mean 


pull width 




1 






0.007 ± 0.032 


0.995 ± 0.023 


179 ±19 


2 


/X 


[GcV/c2] 


0.021 ± 0.033 


1.031 ± 0.024 


5.621 ± 0.003 


3 


a 


[GeV/c2] 


0.026 dz 0.031 


0.976 ± 0.022 


0.023 ± 0.001 


4 


-^B4prong 




0.002 ± 0.032 


1.018 ± 0.023 


150 ± 32 


5 


A^OB 




0.038 ± 0.033 


1.046 ± 0.024 


3170 ± 291 


6 


iVoL 




-0.048 ± 0.033 


1.030 ± 0.023 


962 ± 324 


7 


A^combg 




-0.023 ± 0.032 


1.013 ± 0.023 


1971 ± 171 


8 


Po 




-0.027 ± 0.032 


1.010 ± 0.023 


0.63 ± 0.10 
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Table 5.11: Fixed parameters in the A;, — > A+tt unbinned likelihood fit. 



Parameter 


Meaning 


Value 


/AcK 


NA,-.A+K-/NA,^A+n- 


0.080 


fi 


fraction of the narrow AqK Gaussian 


0.902 




mean of the narrow Gaussian [GeV/c^] 


5.573 


CTl 


width of the narrow AqK Gaussian [GeV/c^] 


0.029 


I^AcK 


mean of the wide AqK Gaussian [GeV/c^] 


5.529 




width of the wide AqK Gaussian [GeV/c^] 


0.075 


h 


fraction of the Landau, 4-prong 


0.413 


A*BPL 


mean of the Landau, 4-prong [GeV/c^] 


5.486 


Ol 


width of the Landau, 4-prong [GcV/c^] 


0.025 


/^BPG 


mean of the Gaussian, 4-prong [GeV/c^] 


5.526 


og 


width of the Gaussian, 4-prong [GeV/c^] 


0.078 


So 


slope of the exponential, other B 


2.180 


/bifg 


fraction of the bifurcated Gaus, other B 


0.106 




mean of the bifurcated Gaus, other B [GeV/c^] 


5.598 


t^ob 


left a of the bifurcated Gaus, other B [GeV/c^] 


10.0 


<yl 


right a of the bifurcated Gaus, other B [GeV/c^] 


4.800 


l^obst 


mean of "step-down", other B [GeV/c^] 


5.436 


„ob 
"0 


slope of the "step-down" , other B 


0.079 


A*ol 


mean of the bifurcated Gaus, other Af, [GeV/c^] 


3.469 




left a of the bifurcated Gaus, other A^ [GeV/c^] 


10.0 




right a of the bifurcated Gaus, other A^ [GeV/c^] 


1.236 


A*olst 


mean of "step-down", other A;, [GeV/c^] 


5.451 


at 


slope of "step-down", other A^ [GeV/c^] 


0.091 


ff 


fraction of first Gaus, other A^ 


0.0005 




mean of first Gaus, other A;, 


5.644 




width of first Gaus, other A^ 


0.019 


fdl 
J2 


fraction of second Gaus, other A;, 


0.0034 




mean of second Gaus, other A;, 


5.459 


^2 


width of second Gaus, other A^ 


0.030 
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5.3 Summary 



Using the unbinned, extended log likelihood technique, we fit the charm and B 
hadron mass spectra to obtain the number of events. The yield for each mode 
is hsted below. The performance of the fitter is vahdated using 1000 toy MC 
test for each mode. In general, the mean of each pull distribution from the toy 
MC test is consistent with zero and the width is consistent with one. For the 
fit parameter with a pull mean deviated from zero and a width deviated from 
unity, the fitter only indicates a less than 1% bias on the central value. Besides, 
these fit parameters are not correlated with the number of signal events and do 
not affect the yield we obtain. The fit result to the data is also superimposed on 
the data histograms and a is computed. We have obtained good ior each 
mode. 



Mode 


Yield 


5" ^ D*+7r- 
B D*+fjL-X 
B° D+TT- 
B^D+ii-X 
Ab A+TT- 
B ^ Atpi-X 


106 ± 11 
1059 ± 33 

579 ± 30 
4720 ± 100 

179 ± 19 
1237 ± 97 
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Chapter 6 



Monte Carlo Samples, 
Acceptance and Efficiencies 

With the raw yield in hand, we now turn to the correction which must be apphed 
to obtain the value of the ratio of branching fractions, that is the acceptance, 
trigger and reconstruction efficiency which may only be calculated using a Monte 
Carlo program. The Monte Carlo (MC) simulation plays a crucial role in this 
analysis. In addition to the acceptance and efficiencies for our signals and back- 
grounds, as described in Section 6.3 and Chapter 7, the MC is used for the 
optimization of signals in Section 4.2. MC is also used to find out the function 
form that describes the mass spectrum of the background due to partial- or mis- 
reconstruction in Section 5.2. In this chapter, we first explain the components 
of Monte Carlo samples and show that, in general, the MC reproduces the data. 
Then we present the acceptance, trigger and reconstruction efficiencies obtained 
from the MC. 

6.1 Monte Carlo Simulation Components 

There are several components in the MC simulation: 

• production and decay of the B hadrons 

• detector simulation 

• trigger simulation 

Production and Decay of B Hadrons 

We use two types of event generators: Bgenerator [70] and PYTHIA [71]. The 
Bgenerator is the primary generator used in this analysis for calculating the 
acceptance and efficiencies of our signals and most backgrounds. Bgenerator 
generates a single &-quark according to the Prib) spectrum which follows the 
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NLO calculation by Nason, Dawson, and Ellis (NDE) [72]. The rapidity of the b 
quark, y{b), is generated according to 1 + qiPr), where 



g(Py) = -0.0456 - 0.00289Pt. (6.1) 

The mass of the b quark is set to 4.75 GeV/c^. For the B meson MC sample, 
the b-quarks are generated with a Pt threshold of 4.0 GcV/c over the range in 
rapidity \y\ < 2.5, and then fragmented into B mesons with the CDF default 
Peterson fragmentation parameter [73], e^, set to 0.006. Figure 6.1 shows a small 
discrepancy in the reconstructed Pt{B^) between data and MC. The slope of the 
data to MC ratio is about 2 a away from zero. The MC events which survive 
the trigger simulation, reconstruction and the analysis cuts, will be rc-wcightcd 
according to the ratio numerically, i.e. we multiply each event with the ratio, w. 
We then calculate the efficiencies using the re-weighted MC events; 

-Rpass = (6.2) 
i 



i-pass 
'gen 



(6.3) 



Figure 6.1 also shows a Aa discrepancy in the reconstructed Pri^b) between 
data and MC from the Bgenerator. As the discrepancy is significant, in order 
to correctly assess the acceptance and efficiency of the A5, the fragmentation 
process inside Bgenerator has to be turned off. The needs to be gener- 
ated directly with a Pt spectrum which reproduces the data. This spectrum is 
obtained in the following way: We first obtain the default generated A;, Pt spec- 
trum from the Bgenerator. Then, the default generated Pt{Ai,) is rc-wcighted 
with the exponential slope of the ratio data/MC shown in Figure 6.1, using the 
"acceptance-rejection (Von Neumann)" method [74]. See Figure 6.2 for the A^ Pt 
spectra before and after our reweighting. We also confirm that the reconstructed 
PT{^b) from the MC using the re-weighted spectrum reproduces the data, see 
Figure 6.10. 

The other event generator, PYTHIA, is a program for the generation of collisions 
at high energies. PYTHIA simulates physics processes using leading-order matrix 
elements, supplemented by the initial and final state radiation. The program 
also includes fragmentation and hadronization of the quarks and gluons in the 
final state. Unlike Bgenerator, PYTHIA includes the beam remnants that are 
left when a parton from the beam particle is removed to participate in the hard 
QCD interaction. PYTHIA provides more realistic simulation of an event than 
Bgenerator, and produces multi-particle final states similar to the hadron collider 
data. However the generation using the PYTHIA is also more time consuming than 
the Bgenerator. This makes PYTHIA inefficient to understand the acceptance and 
efficiency of a single decay mode. Therefore, PYTHIA has been used in this analysis 
only to study the background from bb and cc decays. 
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After the event generation, the hadrons are allowed to decay using the EvtGen 
software package developed by Lange and Ryd [75]. This package is maintained 
by BABAR and mainly tuned by the results from the experiments at the T(4S') 
resonance. The decay model and branching ratios for and are well de- 
scribed but not necessarily those of the Bs and the B baryons. As a proper decay 
model for the A;, semileptonic decays is not implemented in the EvtGen yet, we 
use a flat phase space to decay the A5 first. Then, we will apply a scaling factor 
on the acceptance after taking into account the effect of the semileptonic form 
factors (see Section 6.3). 

Detector, Trigger Simulation, Production and Reconstruction 

The particles from the output of Bgenerator and EvtGen are then run through a 
full ( "realistic" ) simulation of the CDF detector and trigger. The software version 
for the simulation is 4.11.2 with patches which implement the most up-to-date 
configuration of SVT. The geometry and response of the detectors active and 
passive components arc simulated using the GEANT software developed by Brun, 
Hagelberg, Hansroul, and Lassalle [76]. Most of the detector subsystems, like 
COT and CMU, are assumed to be in a time-independent and perfect condition, 
which means there are no dead channels and the high voltages are constantly at 
full value. Selecting the data when these systems are in good condition helps to 
ensure that MC reproduces the data, see Section 4.1. Because the SVX active 
coverage and the configuration for the XFT and SVT systems change on various 
occasions (see Section 3.5), we divide the data taking period into eight sub- 
periods, where the detector and trigger performance is constant. We generate 
our MC samples for these eight sub-periods by choosing the runs with maximum 
number of L3 triggered events as the representative runs. Each run has its own 
parameters for the performance of the detector and triggers. For the sub-periods 
with large integrated luminosity, we choose more representative runs so that 
each run corresponds to a period with integrated luminosity around 3-6 pb~^. 
See Table 6.1 for the representative runs in the MC. The number of generated 
events is proportional to the integrated luminosity of the sub-period each run 
represents. The positions of the beamline for each run is taken directly from the 
database and simulated in the MC. 

After the detector and trigger simulation, the MC events are run through a 
trigger decision program, svtf ilter. svtf ilter takes the information from the 
simulated SVT data and makes the B_CHARM Scenario A requirements described 
in Section 3.5.4. The events which pass svtf ilter are processed with the same 
Production executable (version 4.9.1hpt3) as that which is run on the data. 
The Production executable reconstructs higher level objects, such as electrons, 
muons, tracks and missing energy, from the simulated detector and trigger data. 
The resulting MC events have the same structure and format as the data and are 
then run through the same analysis program described in Section 4.1. 
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Figure 6.1: Comparison of reconstructed -B° (left) and (right) Pt spectrum 
between data and Bgenerator (MC). The top figures show the Pt distribution 
while the bottom figures show the ratio data/MC. The curves in the bottom 
figures are the result of an exponential fit to the ratio. It is evident that the MC 
Pt spectrum is harder than that of the data. 



xlO' 




Generator level P ^{A^^) [GeV/c] 



Figure 6.2: Generator level Af, Pt spectra before and after the reweighting. 
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Table 6.1: Simulated Runs in the MC sample. 



Range 


Run 


JCdt (pb-^) 


Comment 


138809-143000 


140129 


3.4 


Scenario A implemented 


143001-146000 


145005 


4.0 


Tevatron incident 


146001-149659 


148824 
149387 


4.2 
2.9 


SVX coverage improved 


149660-150009 


149663 


0.6 


SVT optimization (covcrage+patterns) 


150010-152668 


150820 
151844 

152520 


4.1 
3.7 

3.5 


-^xy ^ 200 fim cut added 


152669-156487 


152967 
153327 
153447 
153694 

154452 
154654 
155364 
155795 
155895 
156116 
156484 


3.6 
3.7 
3.7 
2.4 
4.2 
4.9 
4.3 
2.5 
3.6 
3.7 
2.6 


XFT from 2-miss to 1-miss 




lUUZOU 

1 fil ^7Q 
lUlU ( o 

1 fi9i "^n 

±u^ ±ou 

163431 


O. 1 

7 
o.o 
o.o 
o.y 
o.u 
o.u 

(J.U 

7 
o . ( 

o. / 

4.3 


data taken after the shutdown 


164303-167715 


164451 


4.6 


SVT change from 4/4 to 4/5 




164844 


3.5 




165121 


2.9 






165271 


3.9 






165412 


3.6 






166008 


6.0 






166063 


2.9 






166567 


5.2 






166662 


5.3 






167053 


5.9 






167186 


2.2 






167506 


4.0 






167551 


2.7 




Total 


170.9 
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M(D^)[GeV/c^] 



Figure 6.3: Md^ from the B D^fifake data. The distribution is sideband 
subtracted using Mk-ktt- 



6.2 Monte Carlo and Data Comparison 

To confirm that the simulation accurately reproduces the data, we compare var- 
ious reconstructed distributions from the MC with the same distribution from 
the data. To ensure a fair comparison, the combinatorial background present in 
the signal region of data has to be removed . We perform a sideband subtraction 
for the 5° ^ D*+Ti~, B D*+fi-X, 5° D+n- and B D+fi-X decays. 
For the A;, A^tt^ and B —>■ A^fi~X decays, a sideband subtraction can not 
remove all the backgrounds in the signal region as explained later in the text 
and in Section 5.2. Instead, a signal distribution of variable "X" is obtained by 
fitting Ma^tt and MpXn to get the number of signal events in bins of variable 
"X". For all the semileptonic modes, we include the MC samples of the physics 
backgrounds described in Section 7.2. The distribution from each physics back- 
ground is scaled according to the assumed or measured branching ratio for that 
background. In addition, the distribution of each compared variable from the 
fake muons is subtracted from the data. The distribution from the fake muons 
is obtained by reconstructing the "fake muon-charm" final state as described in 
Section 7.3. The combinatorial background in the "fake muon-charm" is removed 
using the same method as described above for the real muon. See Figure 6.3 for 
the Md^ from the muon fakes. 

For the B meson semileptonic channels, the mass difference between D*~^ 
and Z)° in the B D*~^n~X mode {Mk-k-k - Mk-k), and mass of in the 
B — s> D^fi^X mode {Mkittt), are used as the variables to perform the sideband 
subtraction. The signal region for both these modes is defined as: 

|M-Mpdg| < 2(T, (6.4) 
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Figure 6.4: Invariant mass difference Mk-ktt — Mk-k (left) and invariant mass Mkittt 
(riglit) sfiowing our reconstructed B D*^fi^X and B D^^^X signals. The 
vertical solid (dashed) lines indicate the signal (sideband) regions. 



and the sideband region is defined as: 

4cr < |M-Mpdg| < 6cT. (6.5) 

The background function is assumed to be a straight line. Therefore, the amount 
of background in our signal region is the same as that in our sideband regions. We 
obtain a clean signal distribution by subtracting the histogram in the sideband 
region from the histogram in the signal region. Figure 6.4 displays the signal and 
sideband regions of Mk-k-k - Mxn and Mktttt- 

For the B meson hadronic modes, we use the upper mass sideband above the 
signal peak to perform the sideband subtraction. The lower mass region below the 
signal peak consists of both combinatorial background and partially reconstructed 
B decays. However, the background in the signal region and in the upper mass 
region above the peak is mainly combinatorial as shown in Figures 5.11 - 5.13. 
We have learned in Section 5.2 that the combinatorial background is adequately 
described by an exponential function. Therefore, we fit the upper mass region 
to an exponential function. We further extrapolate the exponential to the signal 
region and obtain the ratio of the background in our signal region to that in our 
upper mass sideband, i?bg- The histogram of the compared variable extracted 
from the upper mass sideband is scaled by i?bg and subtracted from the histogram 
in the signal region. See Figure 6.5 for the B^ mass signal region we define and 
the upper mass region we fit to an exponential. 

For the B A'^fx'X and A;, A^tt" modes, there are non-negligible back- 
grounds under the signal peak from the refiections due to a mis-assignment of 
the mass for one of the particles, see Chapter 5 for more details. This type of 
background has a different behavior from the combinatorial background in the 
sideband region. Since a background-free sideband subtraction is difficult to per- 
form, we choose to fit the number of signal events in each bin of the variables 
which we want to compare. For the number of B ^ A+/i~X candidates, the 
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Figure 6.5: Invariant mass Mk-k-rtt for B D*^-k^ (top) and B D^n^ 
(bottom) signals. The vertical solid lines indicate the signal region. The upper 
mass regions on the right are fitted to exponential functions. 



MpKn distribution is fitted to a signal Gaussian and a second-order polynomial 
background as shown in Figure 6.6. For the number of Af, —>■ A^Tr^candidates, 
the MpK-K-K distribution is fitted to a simplified model: a Gaussian signal and 
an exponential background, as shown in Figure 6.7. Note that although the A;, 
fit model is simplified, the systematic uncertainty due to the naive model is no 
more than 3% of the number of signal events in each bin compared with the 15% 
statistical uncertainty. The widths of Mp^n and M^n are fixed to the values 
obtained from the full statistics when doing the fit. Figure 6.8 shows the data 
and MC comparison using the fit values obtained from Figures 6.6- 6.7. 

When comparing the MC and data distributions, if the number of data signal 
events in one bin is less than 20, that bin is combined with the next bin until the 
sum of the events is over 20. Then a is computed. 



where i stands for i*'' bin and total number of bins in a histogram is n. The 
number of degree freedom is n — 1. For the B — > A^n~X and Af, — > A^vr" 
modes, a is also calculated except that the bin width of each variable is fixed 
in this case. Besides the test, we also plot the ratio data/MC. We fit the ratio 
to a first-order polynomial and check if the slope, Ai, is consistent with zero. 

In the first pass, we find discrepancies in the Pt spectra of B^ and A;, between 
MC and data (see Figure 6.1). As the semileptonic modes are three-body decays 
and the hadronic modes are two-body decays, the efficiency of the trigger and 
analysis Pt cut depends strongly on the Pt of B hadron (see Figure 6.9). We 
decide to reweight the Pt spectra of -B° and Af, as described in Section 6.1. 
Figures 6.10- 6.11 show the comparison between MC and data for the analysis 
cut variables of the A;, — > A^ modes. Figure 6.12 shows the comparison for the 
Md*^, Md^, and the Ma^^ from the phase space MC before and after multiplying 
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each bin entry with a scahng factor. The scahng factor is obtained by dividing the 
Ma^^ distribution from the form factor weighted (see Section 6.3) by that from 
the phase space generator-level MC. The agreement of the MC Ma^^ distribution 
with that from the data has significantly improved after applying the scaling 
factor. In addition, we compare the efficiency of the M^^^ cut in the MC and 
data, given that the other (N-1) analysis cuts are apphed. The contribution of the 
fake muons and physics backgrounds are included. Since there are uncertainties 
from the fit to both data and MC, the efficiency is defined as: 



where n is the number of events after making all the analysis cuts and N is the 
number of events after making the A'^ — 1 cuts. The uncertainty on the efficiency 
is derived by Heinrich [77]: 

= y(6-^)2+(l-26)-(^)2, (6.8) 

where and ct^ are the uncertainties from the fit. We find the data give an 
efficiency of 0.77±0.04, while the form factor weighted MC gives an efficiency of 
0.72±0.05, which is in good agreement with the data efficiency. Other distribu- 
tions important for this analysis may be found in Yu [65] . 

In general, the MC describes the data well except for the pseudo cr of Ac/i, 
and the Xr-,p of the B and charm vertex fits. For the disagreement in the vertex 
fit 

beyond the scope of this analysis to scale the measurement errors 
in the MC, we choose to make a loose cuts on the data Xr-</>- Section 8.1.3, 
we perform a cross-check of the relative branching ratio variation by dividing the 
data into two subsets, according to the cuts on the pseudo cr of Ac/i, xl-(p 
other variables. We do not sec significant inconsistency. Therefore, we do not 
assign any systematic uncertainties for these variables. 
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Figure 6.6: Example of mass fit for the MC and data comparison. The 
variable to compare is the Pt of proton, from 2 to 8 GeV/c, in bins of 0.5 GeV/c. 
MpKn is fitted to a signal Gaussian and a second-order polynomial. 
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Figure 6.7: Example of Ab mass fit for the MC and data comparison. The variable 
to compare is the Pt of proton, from 2 to 10 GeV/c, in bins of 1 GeV/c. Mpx-Kn 
is fitted to a signal Gaussian and an exponential background. 
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Figure 6.8: MC and data comparison of PT(proton). The data points come from 
the fit to MpK-n in Figure 6.6 (left) and the fit to MpKmr in Figure 6.7 (right). 
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Figure 6.9: MC Efficiency ratio of A;, — > A+tt to A;, A+/x z/^ as a function of 
the transverse momentum of A^. 



89 




Figure 6.10: A+tt^ MC and data comparison: from the top left to the 

bottom right are: PT{Ah), cr(Afe), vertex fit vertex, Pt{A^), 

cr(A^), and vertex fit xl^s the A^ vertex. 
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Figure 6.11: B A^/i~X MC and data comparison: from the top left to the 
bottom right are: PriAcfi), cr(Ac/i), vertex fit Xr-if, ^^i the Ac/i vertex, Pt{A'^), 
cr(A+), and vertex fit Xr~6 for the A+ vertex. 
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Figure 6.12: four track invariant mass MC and data comparison: from the top left 
to the bottom right are: M^*^, M^/x, ^Acm (phase space MC without scahng), 
Ma^^ (phase space MC after scahng). 
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6.3 Acceptance, Trigger and Reconstruction Ef- 
ficiencies of Signal 

We obtain the product of the acceptance, trigger and reconstruction efficiencies 
using the MC described in Section 6.1. The total efficiency is defined as: the 
number of reconstructed events surviving the trigger simulation and analysis cuts, 
divided by the total number of events generated. Efficiencies for the backgrounds 
are found in Chapter 7. Our data could be divided into eight sub-periods under 
different trigger and hardware configurations, see Table 6.1. In this analysis, 
because the final states are nearly alike, we expect that the ratio of the efficiencies 
to be independent of the detector, trigger and cahbration effects. To confirm 
this, we divide our signal MC samples into eight sub-periods and calculate their 
efficiencies and the ratio of hadronic to semileptonic modes. Tables 6.2- 6.4 list 
the efficiencies for our signals. Figure 6.13-6.14 show that the absolute efficiency 
may vary dramatically in each period, however the efficiency ratio of the hadronic 
to the semileptonic mode is quite stable. 

Note that because Bloom and Dagenhart [78] find a difference in the CMU 
muon reconstruction efficiency between MC and data, wc apply a scaling factor on 
the efficiencies of the semileptonic signals and backgrounds: R = 0.986±0.003. In 
addition, Giagu, Herndon and Rescigno [79] [80] notice that there are differences 
in the XFT efficiencies for the charged kaons, pions, and protons, when the XFT 
configuration is switched to the "1-miss" mode, i.e. when the tracking algorithm 
in the XFT requires at least 11 hits of 12 wires from each COT axial superlayer. 
The COT frond-end electronics requires a minimum input charge from the ion- 
ization of the incident particle. At a fixed momentum, protons and kaons deposit 
less charge than the pions, have more hits below the electronics threshold, and 
fail the stringent XFT "1-miss" requirement. Therefore, in general, the proton 
and kaon XFT efficiencies are lower than that of the pion. Figure 6. 15 (top) shows 
that kaon and pion XFT efficiencies are identical in the MC and need to be cor- 
rected. Giagu, Herndon and Rescigno measure the ratio, data/MC for the XFT 
efficiencies of pions, kaons and protons, as shown in Figure 6.15. We rewcight 
the MC events numerically according to the ratio: 

C, = 1.002 -"-"^^ 



Ck = 0.969- 



Pt ' 
0.094 



T 



1.3 3.2 2.2 
C, = 1.06- — -f-^--3, 

IJ' Irp JTrp 

where Pt is the transverse momentum of the track that passes the trigger cuts 
in our reconstruction program. 

Finally, one additional scaling factor has to be apphed on all the decays 
with Ac// in the final state. We have mentioned in Section 6.1 that a phase space 
decay model was used for these decays. In a phase space, the event density in the 
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u^-cos 9 plane is a constant within the kinematic boundary. The w is the scalar 
product of the A5 and four-velocities, and 9 is the angle between the muon 
and the neutrino momentum vectors in the A5 rest frame. The form factors that 
describe the strong interaction in the A;, semileptonic decay modify the event 
distribution in the phase space and change the fraction of events accepted. We 
obtain the scaling factor in the following way: Using the "acceptance-rejection" 
method, we reweight the generator-level A^ — > A+//~I7^ MC according to: 



where the differential semileptonic decay rate, — ^^dw ^ '^^ ' obtained from 
Huang [36]. The T{cos9,w) inchides the W spin effect and describes the corre- 
lation between the /x and V^, and 



Here, cos^max and cos^min specify the kinematic range and are functions of w. 
Figure 6.16 shows the phase space and the w distribution from the phase space 
and the form factor reweighted MC. Then, we apply generator-level analysis cuts 
to obtain the acceptance. We further divide this acceptance by that from the 
phase space MC and obtain a scaling factor of 0.994±0.025, where the uncertainty 
is dominated by the size of the MC sample. 

6.4 Summary 

We have described the procedure of generating MC samples and compared the 
MC distributions with those in the data. In general, the MC and data are in good 
agreement. We also obtain the signal efficiencies from the MC. It is confirmed 
that the efficiency ratios from both the A;, and B° modes are insensitive to the 
time variation of beam lines and SVT trigger configurations. 



dr{Ab A+/x-I7^) T(cos 9, w) 
dw P{w) 



(6.9) 




(6.10) 
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Table 6.2: Total efficiency and ratios for B — > D*^t: and B — > D*^jji 



Run Range 


J UjV 






C IVCttliJ 






(pb ^) 


(10 ^) 


(10 ^) 




138809- 


-143000 


3.4 


1.72 ± 0.15 


0.77 ± 0.07 


2.22 ± 0.28 


143001- 


-146000 


4.0 


1.42 ± 0.12 


0.73 ± 0.06 


1.94 ± 0.24 


146001- 


-149659 


7.1 


1.57 ± 0.10 


0.82 ± 0.05 


1.92 ± 0.17 


149660- 


-150009 


0.6 


2.35 ± 0.40 


1.13 ± 0.19 


2.08 ± 0.50 


150010- 


-152668 


11.3 


2.65 ± 0.10 


1.13 ± 0.05 


2.34 ± 0.13 


152669- 


-156487 


39.2 


2.74 ± 0.06 


1.23 ± 0.03 


2.23 ± 0.07 


159603- 


-164302 


52.7 


3.10 ± 0.05 


1.38 ± 0.02 


2.24 ± 0.05 


164303- 


-167715 


52.7 


4.27 ± 0.06 


1.90 ± 0.03 


2.25 ± 0.05 


Total average 


171.0 


3.22 ± 0.03 


1.44 ± 0.01 


2.24 ± 0.03 



Table 6.3: Total efficiency and ratios for B —>■ D~^tt and B jj, 



Run Range 


JCdt 






e Ratio 






(pb-i) 


(10-^) 


(10-^) 




138809- 


-143000 


3.4 


2.99 ± 0.19 


1.35 ± 0.09 


2.21 ± 0.21 


143001- 


-146000 


4.0 


2.39 ± 0.16 


1.36 ± 0.09 


1.76 ± 0.16 


146001- 


-149659 


7.1 


2.93 ± 0.13 


1.28 ± 0.06 


2.30 ± 0.15 


149660- 


-150009 


0.6 


4.82 ± 0.57 


2.23 ± 0.27 


2.16 ± 0.37 


150010- 


-152668 


11.3 


4.12 ± 0.14 


2.01 ± 0.06 


2.05 ± 0.09 


152669- 


-156487 


39.2 


4.79 ± 0.08 


2.24 ± 0.04 


2.14 ± 0.05 


159603- 


-164302 


52.7 


5.43 ± 0.07 


2.48 ± 0.03 


2.19 ± 0.04 


164303- 


-167715 


52.7 


7.49 ± 0.08 


3.37 ± 0.04 


2.22 ± 0.03 


Total average 


171.0 


5.67 ± 0.04 


2.58 ± 0.02 


2.20 ± 0.02 



Table 6.4: Total efficiency and ratios for — > A+tt and A^ — > A+// u 



Run Range 


JCdt 
(pb-i) 


(10-^) 


(10-^) 


e Ratio 


138809- 


-143000 


3.4 


1.84 ± 0.15 


0.50 ± 0.06 


3.69 ± 0.51 


143001- 


-146000 


4.0 


1.23 ± 0.12 


0.36 ± 0.04 


3.43 ± 0.53 


146001- 


-149659 


7.1 


1.42 ± 0.09 


0.48 ± 0.04 


2.98 ± 0.30 


149660- 


-150009 


0.6 


2.23 ± 0.39 


1.13 ± 0.19 


1.98 ± 0.48 


150010- 


-152668 


11.3 


2.13 ± 0.12 


0.63 ± 0.04 


3.37 ± 0.27 


152669- 


-156487 


39.2 


2.37 ± 0.05 


0.75 ± 0.02 


3.15 ± 0.11 


159603- 


-164302 


52.7 


2.67 ± 0.05 


0.82 ± 0.02 


3.26 ± 0.09 


164303- 


-167715 


52.7 


3.76 ± 0.06 


1.14 ± 0.02 


3.30 ± 0.08 


Total average 


171.0 


2.86 ± 0.03 


0.86 ± 0.01 


3.31 ± 0.05 
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Figure 6.13: B D^n MC total efficiency (left) and the ratio of that to the 
— > D+yU~z7^ MC efficiency (right) in eight different hardware configurations. 
The shaded area represents the average efficiency ratio including the uncertainty. 
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Figure 6.14: Total efficiency ratio of B D*+n- to B D*+ fi (left) and 
Af, A+7r~ to A;, — i> A+/i~z7^ (right) MC in eight different hardware configura- 
tions. The shaded area represents the total average efficiency ratio including the 
uncertainty. 
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Figure 6.15: Giagu and Rescigno [79] find that K and vr XFT efficiencies in the 
MC (top left) are identical. The relative proton (top right) and vr (bottom) 
XFT efficiencies between MC and data, in bins of 1/Pt, are fit to a third or- 
der polynomial by Herndon [80], and to a ffist order polynomial by Giagu and 
Rescigno [79], respectively. 
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Figure 6.16: The phase space (top) and the w (bottom) distribution obtained 
from the A;, — >■ A+/i~I7^ phase space MC, before (left) and after (right) reweight- 
ing the events according to the semileptonic form factors from Huang [36]. 
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Chapter 7 

Backgrounds of the Semileptonic 
Modes 



Unlike the e^e~ collider experiments, BELLE [24] and BABAR [25], the initial en- 
ergies of the B hadrons are unknown in CDF. For the hadronic channels, such 
as B D+TI-, the B meson is fully reconstructed using the momenta of the 
daughter particles {D^ and tt^) at the point of the parent decay. However, for the 
semileptonic channels, such as S ^ £)+//~I7^, without the information of initial 
energies, it is difficult to constrain the momenta of the missing neutrinos and fully 
reconstruct the B candidates. Therefore, the D* and Ac// combinations we 
observe in the data consist of the exclusive semileptonic signals, B Z}*+/i^77^, 
B — > jjL^V^, Afo — i> Iv^ji^V^, in the presence of other backgrounds. These 
backgrounds arise from three sources: 

• physics backgrounds: B hadron decays into similar final states, a charm 
hadron, a real muon and other tracks. 

• muon fakes: a charm hadron and a track which fakes a muon. 

• &6 and cc: two B or charm hadrons from the bb and cc pairs decay into a 
£)*+ A+) and a muon, respectively. 

The goal is to measure the relative partial decay widths of the exclusive semilep- 
tonic decays to hadronic decays. The backgrounds listed above should be sub- 
tracted from the observed inclusive semileptonic signal in the data. The ratio of 
branching fractions is then calculated as follows: 

^semi _ / -^inclusive semi ~ -^physics ~ ^faken ~ ^ cc, bb\ ^had ^'7 1 ^ 

~R ~^ N ' ~e ' \'-^) 

'-'had ^ ' had tgemi 

where B stands for the branching ratio, e is the efficiency from the MC. We 
estimate the contribution from the physics and 66, cc backgrounds, using the 
efficiencies from the MC and the branching ratios from the Particle Data Group 
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(PDG) [8]. We normalize the backgrounds to the observed number of events in 
the fully reconstructed hadronic signal in the data, 



-^physics {bb,cc) — -^had " XI » Z " 

*-'had ■ Chad 

Substituting Equation 7.2 into Equation 7.1, A^had cancels. The estimate of bb, cc 
and physics background contributions is then free from the uncertainties in the 
hadronic yields. In the case of fake muons, we subtract the fake muon candidates 

measured in the data directly. 

The branching ratios of hadronic modes in Equation 7.2 come from the 
world average in the PDG. For the A;, mode, we extract B{Ah A+tt^) from the 
recent CDF result, '^H^^T>G-o)B{Ab-*AcTv ) , -^^ rgg^ gg^g^^gg Lg'g measure- 

ment requires the ratio of the A;, to production cross sections as an input, we 
correct for the ratio in Section 7.1 and obtain i3(A(, —>■ A+tt"). Sections 7.2- 7.4 
estimate the amount of backgrounds in the semileptonic signal. We will show that 
the dominant signal contamination is from the physics background. The second 
largest background arises from muon fakes. The smallest background source is 
from bb and cc. 



7.1 Correction of ^ and B(Ab A+tt") 

CDF has made a number of measurements of the relating information from the 
decay of 5-hadron to that of another with similar topology, such as 

aAjPT>6.0)g(Ab^A+7r-) r aA,(PT>4.0)g(A6^J/V-A) r^^. 

a^o{PT>6.0)B(B"^D+n-) " (TBo{PT>4.0)S(BO-.J/VifO) L^-'J" 

In order to extract the individual branching fraction, one must have the knowledge 
of the production cross section. Since no production cross section measurements 
exist for Aj,, one must infer the cross-section using other means. Using the total 
6-quark cross-section [41], the fragmentation fraction (/„, fd, fs and /baryon), and 
assuming the composition of the 6-baryons is dominated by Af,, i.e. /a^ = /baryoiu 
we may infer the production cross-section of any 5-hadron species. However, the 
fragmentation fractions assume the entire Pt spectrum is measured for both the 
particles and A;,) in the above ratios. If the Pt spectrum is incorrect, the 
kinematic acceptance used in the above ratios will be over- or under-estimated. In 
addition, since most analyses at CDF require a Pt threshold to improve the sig- 
nal to background ratio, any difference in the Pt spectrum between the particles 
participating in the ratio will modify the effective fragmentation fraction. Fig- 
ure 7.1 illustrates this effect for a very small difference in the spectrum (top) and 
a large difference in the spectrum (bottom). In this section, we derive corrections 
for the production cross sections ^— ^ to account for their different Pt spectra and 
overestimated kinematic acceptance from the previous measurements. 
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Figure 7.1: Pt spectrum of and A;, (left) and the dependence of the production 
ratio on the Pt threshold, (right). The top figures show the case where 

both hadrons have an almost identical spectrum. The ratio of the integrated 
areas underneath the spectrum, from and above, does not depend much on 
the value of P^™. The bottom figure shows a softer A;, Pj- spectrum. The ratio 
of the integrated areas will depend strongly on the value of P|f^™. 



101 



To simplify the notation, we define the relative production cross sections of 
Aft to as p, 

= ^Ati^T >xGcY/c) 
^^^^-^^..(Pt >xGeV/c)' ^^■'^> 

where x is the Pt threshold of the B hadrons in GeV/c. When estimating the 
backgrounds for our semileptonic signals, we choose to use the hadronic signals 
as the normalization because they are fully reconstructed, and free from indistin- 
guishable physics backgrounds. Efficiencies and branching ratios of the hadronic 
modes are required when performing the normalization, see Equation 7.2. We use 
Monte Carlo to calculate the efficiencies and external measurements of hadronic 
decays for the branching ratios. While precise measurements of 13{B — > D*+7r~) 
and B{B^ — > D+tt") exist [8], there is not yet a direct measurement of B(Ab — > 
A+TT-). The CDF Run II measurement of p(6) x ^'^^Zm'^J-] [^^^ 
provides an input but requires a good understanding of p(6). In addition, we 
need to know p(4) to normalize the B meson background, 5°'"*" — > A+/xX in the 
inclusive Ac/i signal, to A;, A+tt", 

A lower minimum Pt requirement of 4 GeV/c is applied here to cover the entire 
kinematic range of the B hadron semileptonic decays, see Section 7.2.2 for more 
details. 

The fragmentation fractions fd, fu, fs and /baryon in the PDG are defined as 
the probability for a b quark to hadronize into B^, B~^, B^ or baryons, 

(7.5) 
(7.6) 
(7.7) 
(7.8) 

with the assumption that 

/d + /« + /. + /baryon = l. (7.9) 

The PDG [82] determines /baryon by combining the measurements of mean life- 
times, S° mixing parameters and branching ratios by LEP, SLD and Taylor in 
CDF-I [83]. Nevertheless, it is still questionable whether the /baryon should be the 
same for the LEP and the Tevatron experiments. As the collider and the detector 
environments for Taylor's analysis are most similar to those for this analysis and 
we have access to the details of the measurement, we discuss the corrections for 
Taylor's /baryon//d rcsult, 0.236 ± 0.084. 



fd 


= B{b- 




fu 


= B{b- 




fs 


= B(b- 




/baryon 


= B(b- 


-> b — baryon) , 
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Correcting the Kinematic Acceptance 

Taylor's analysis uses electron-charm final states, such as S — > £)*+e~X, B — > 
D'^e~X, and B — > A^e~X to measure /baryon- Accurate and B^ Pt spectra 
from the fully reconstructed decays were not available at the time. The MC 
samples [84] in Taylor's analysis for calculating the acceptance and efficiencies 
were generated using the default settings in the package Bgenerator as described 
in Section 6.1. We have shown in Section 6.2 that the and B^ Pt spectra from 
the Bgenerator are stiffer than those measured in the data. This leads to an 
over-estimate of the acceptance and efficiencies, particularly for the A^ decays. 
One must correct for this effect first. 

For technical reasons, it is impossible to repeat a full CDF Run I detector 
simulation. Consequently, we obtain the correction using generator-level MC. 
We generate B° and A^ with the default Bgenerator Pt spectra as described 
in Section 6.1. Then we decay B^ and Af, using the QQ software package [85] 
as in Taylor's analysis. We apply the cuts listed in Table 7.1 on the generator 
level quantities. These cuts mimic those applied in Taylor's analysis as much as 
possible. We divide the number of events which pass the cuts by the number 
of events with B^ (A^) Pt > 10 GeV/c and rapidity \y\ < 2.0, to obtain the 
acceptance and efficiencies for the exclusive semileptonic decays: B —>■ D~^e~Ve, 
and Ab A+e~z7e. We repeat the same process using the corrected B^ and A;, 
Pt spectra as described in Section 6.1, except that we use QQ to decay A;, and 
B^ to be consistent with Taylor's analysis. The production cross-section ratio 
derived from the Taylor's analysis {j^^'^^^"^) could be expressed as a ratio of 
the number of signal events divided by the ratio of the product of the branching 
ratio and efficiency for the A^ — > A^e~Ve mode, to the same expression for the 
5° — > D^e~Ve mode: 

O^^Taylor ^ . . ^Taylor^ 

a BO 

where we use the following shorthand notation: 

Nr = ^^^-^te-V. ^ 

= (7.12, 

B{Ab A+e u^) 

ER = ^ ^^^^ . (7.13) 
Using the same notation, we could express the corrected ratio as: 

(7 A corrected , , /ja Taylor f, 

1^ ^Nr-Br-eT^^^^ ^— -Ce^-^-C, (7.14) 

CfiO a BO /baryon 

where 

C, = e^°'"'"""*"Ve^^^^°\ (7.15) 
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Table 7.1: Generator- level cuts as Taylor's analysis. 



Parameter 


Cut Value 




Prie) 


> 7 GeV/c 


1.206± 0.006 


Transverse energy -E'r(e) 


> 8 GeV/c^ 


1.332± 0.009 


\v(e)\ 


< 1.1 


1.343± 0.012 


all tracks Pt 


> 0.4 GeV/c 


1.420± 0.015 


daughters of charm do /cdo 


> 1.5 


1.428± 0.015 




> 0.5 GeV/c 


1.428± 0.015 


Pt{K) 


> 1.2 GeV/c 


1.588± 0.020 


Pt{p) 


> 2.0 GeV/c 


1.808± 0.026 




> 1 


1.777± 0.037 


M{De) 


< 5.0 GeV/c^ 


1.777± 0.037 


M(Aee) 


< 5.3 GeV/c^ 


1.805± 0.039 



This gives us the first correction factor, (7^=1.81 ± 0.04 (stat) ^^'22 {Pt)- The last 
uncertainty comes from the uncertainty of the exponential slope used to obtain 
the corrected Pt spectrum as described in Section 6.1. Table 7.1 also gives the 
value of Ce after each cut is apphed. Note that using the generator level MC, we 
do not obtain the reconstruction efficiencies for Taylor's analysis. However, the 
reconstruction efficiency is the same for both f^^^^^""^ an^j corrected cancels 
in Equation 7.15. 



Correction due to the Difference in the Pt Threshold 

The second correction is due to a difference in the Pt threshold of Taylor's and 
our analysis. Data used in Taylor's analysis come from an electron trigger which 
cuts on the transverse energy, Et, of electron below 8 GeV/c^ and probe the 
B hadrons with Pt greater than 10 GeV/c. Our data come from a two track 
trigger with a looser Pt requirement and extend the minimum Pt of B hadrons 
down to 4 GeV/c. If the Pt spectra of B^ and Af, are almost identical, the 
value of p will be independent of the Pt threshold of the B hadrons, as shown 
in Figure 7.1 (top). If the Pt spectra of B^ and A;, are quite different, the value 
of p will strongly depend on the Pt threshold, as shown in Figure 7.1 (bottom). 
Figure 7.2 shows that Af, Pt spectrum is softer than that of the B^ . Therefore, we 
need to apply another correction factor, Cp^, on the previous efficiency corrected 
Taylor's result, p(10) = ^-^^^ • C„ to obtain p(4) and p(6) for this analysis. 
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The p{x) is expressed as follows: 

fd 



p{x)^^^-a-Cp,{x), (7.16) 



where 

^^-^^^ - N^o{Pr > xy N^^iP^ > P^^y^-y ^^-^^^ 

Here, a; and P^^^^"^ stand for the Pt thresholds of the B hadrons in our and 
Taylor's analysis. Wc obtain Cp^{x) using the generator level MC produced with 
the corrected Pt spectra of A;, and . About 20M events of A^ and B^ are 
generated. No additional cuts are apphed except that all the B hadrons have 
rapidity less than 2.0. We count the number of B hadrons above 4 and 6 GeV/c, 
and divide that by the number of B hadrons above certain P^'^^^"'^. A scan of 
pTayior ^^^^ g QqY^q ^o 20 GcV/c is performed but the value at 10 GeV/c is 
used in the analysis. We have Cp^(4) = 1.480 ± 0.002 {stat) +J;J^° (Pt), and 
Cpj,(6) = 1.309 ± 0.002 (stat) i'io'^Qg (-fr), where the last uncertainty comes from 
the uncertainty of the A{, Pt spectrum. Figure 7.3 presents Cpj,(4) and Cpj,{6) 
as a function of P^^^^°\ 

After applying corrections Cg and Cp^,, we calculate p(4) and p(6) to be: 

+0 24 

p(4) = 0.63±0.23{stat®syst) _^_^^{Pt), (7.18) 

+0 19 

p(6) = 0.56±0.20{stat®syst) q'-^-^(^t). (7.19) 

Then, we could go back to CDF Run II measurements of p(6) x ^C^fc^^c-^ ) g^^^^^ 
extract B{Ab ^ A^tt"): 

B(Ab ^ A+TT-) ^G-^-B(^^ D+n-), (7.20) 

p(6) 

G = p(6)x^i^^^^. (7.21) 

The values for each of the parameters are hsted in Table 7.2. We find 

B{Ab ^ A+TT") = (o.41 ± 0.19 [stat © syst) ^Jj'Jjg (Pt)) %, 

which is in good agreement with the prediction by Leibovich, et al. [33], 0.45%. 

As noted earlier, B{Ai, A+vr^) will be used later for the normalization 
of the amount of physics, bb and cc backgrounds to the observed number of 
events in our hadronic signals. Several variables, such as hadronic to semileptonic 
efficiency ratios, Cg, Cpr^ and G, depend on the A^ Pt spectrum and are correlated. 
We would like to study the systematics on the relative branching fractions from 
these variables simultaneously. Consequently, the uncertainty due to the A^ Pt 
spectrum is separated from the other systematic and statistical uncertainties. 



where 
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m 



3.5 r 
3|- 
2.5 

2H 
1.5 
1 

0.5 





f I ndf 
Prob 
pO 
p1 



4.052 / 7 
0.7737 
0.328710.118 
-0.0892210.02933 
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Pj of b-hadron [GeV/c] 



Figure 7.2: Comparison of and A;, Pt spectra measured in the data. The 
positive slope (2 a away from zero) of the ratio of to A;, histograms indicates 
that Pt{B^) is harder than Pxi^^b). 




16 18 20 
pTaylor [QeV/c] 



Figure 7.3: correction factor, Cpj,, as a function of P^'^-^'™. The two curves 
represent the correction factor for two Pt thresholds of B hadrons: 4 GeV/c 
and 6 GeV/c, respectively, pj'^-^'™ is the Pt threshold of B hadrons in Taylor's 
analysis. 
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Table 7.2: Parameters for calculating i3(A(, — > A+tt ). 



Paraiiieter Vi 



aiuo 



G 0.82 ± 0.25 (stat syst) ± 0.06 (Pt) 
4.2 ± 1.5 

/baryon 

Cp^(6) 1.309 ± 0.002 (stat) tlfol (^t) 

Ce 1.81 ± 0.04 (stat) toil (Pt) 

B(B^ D+TT-) (2.76 ± 0.25)- 10"^ 

BjAb A+TT-) (%) 0.41 ± 0.19 (stat gyg^) (Pr) 



7.2 Physics Backgrounds 

Physics backgrounds come from the decays of B hadrons into similar final state 
as our semileptonic signal: a D*'^ [D^, A^), a fi~ and missing particles. Branch- 
ing ratios and efficiencies of these physics decays are needed to normalize the 
background contribution to the observed number of hadronic signal events in the 
data; 

-^physics _ J2 Bj ■ €{ 
^had 

For the backgrounds to the 5^ — > D*^ ii~V^ and 5^ — > D^ii~V^ decays, we find 
the modes which give similar final states as our semileptonic signals in the decays 
listed in the Particle Data Group (PDG) summary [8] and the default decay table 
inside EvtGen package [75]. Many decays of B and D mesons have been measured 
by other experiments, such as CLEO, BELLE and BABAR. These measurements serve 
as inputs to the EvtGen decay package. Since BELLE and BABAR also use the 
EvtGen package, they have included decay modes into EvtGen which have not yet 
been measured and estimate the branching ratios. For the backgrounds to the 
A5 A'tfi~Vfj_ decay, none of the B hadrons decays to Ac/U final states have been 
measured in the CDF-1 and other experiments, or estimated inside EvtGen. We 
use the results from the preliminary measurements by Litvintsev, et al. [86] and 
the prediction of PYTHIA to obtain the background branching ratios. 

After we have a hst of decays which share similar final states as our signals, 
we use the generator level simulation to estimate the composition of the inclu- 
sive semileptonic signal from each physics background. Details of this procedure 
maybe be found in the study by Tesarek, et al. [87] [88]. The decays which 
contribute > 1% to the semileptonic signal after trigger-hke and the four track 
invariant mass Md*(/j a^)^ cuts are selected for further consideration. We generate 
each selected decay separately and run through the full CDF detector simulation 
as described in Section 6.1. Then we run the same signal reconstruction program 
used for the data on the MC and divide the number of reconstructed events by 
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the number of generated events to obtain the efficiency. 



7.2.1 Physics backgrounds of — ^ z^^ and — ^ 

Physics backgrounds of — > D*~^n~T'^ and — > D^ii~Vn fall into two cate- 
gories: 

1. Scmileptonic decays of B'^, B^, Bg, which include either additional particles 
(eg: — > D+7r°/i.~I7^) or a higher mass charm meson with subsequent 
decay into our charm signal (eg: B D*~^iJ,~V^, D*'^ D'^n^) 

2. Hadronic 5° decays into two charm mesons, one charm meson decays as 
our charm final state, the other charm meson decays semileptonically (eg: 

Tables 7.3- 7.4 summarize the physics background in 5^ D*^ ^~V^ and 
B D~^fi~T7^ which contribute > 1%. The second column in the table lists the 
measured or estimated branching ratios. All the numbers in parentheses are esti- 
mated uncertainties (100% for B and 5% for charm) for the unmeasured branching 
fractions. The third column fists their efficiencies relative to the hadronic signals 
with statistical errors. The fourth column lists the normalization of each back- 
ground relative to the exclusive scmileptonic signal. The last column lists the 
number of events from each background after multiplying the relative branching 
ratio and efficiencies with the number of hadronic signal in the data, as expressed 
in Equation 7.22. The uncertainty in the last column only includes the statistical 
uncertainty of the hadronic yield. When normalizing backgrounds from B^ and 
Bs to B^ signals, the following fragmentation fractions quoted in the 2004 PDG 
are used: 

b^Bd = (39.7 ±1.3)%, 
b^Bu = (39.7 ±1.3)%, 
b^Bs = (10.7 ±1.1)%. 

Note that most of the backgrounds which contribute > 1% share the same Feyn- 
man diagram as our scmileptonic signal: 

• resonant mode: a spectator quark and a 6 to c quark transition via a virtual 
W boson exchange (eg: background of 5° ^ £)+/i~z7^, ^ — > D*^ix~Vf^) 

• non-resonant mode: with additional uu or dd quark pair created from the 
vacuum (eg: B — > L'"'"7r°/x~I7^) 

Note that the non-resonant modes tend to have smaller branching ratios and 
smaller efficiencies than the resonant mode. The dominant background of S — > 
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D*~^IJ, v^is B — > z/^ where D\ — > D*+7r . The total physics background 
in Tabic 7.3 is about 15% of i? — D*^ ix~X events in the data after all the 
cuts. The dominant background of -B ^~T'^ is B ^ D*~^fi~Vn where 

D*^ — > D~^7T^. The total physics background in Table 7.4 contributes about 40% 
of S — > D^jjrX events in the data after all the cuts. As shown in Figure 7.4, a 
cut on the invariant mass of can reduce or eliminate the background 

from B'^, B^ decaying semileptonically to more particles or a higher mass charm 
state. 

7.2.2 Physics backgrounds of A;, — > A+//~F^ 
Physics backgrounds of — > A^fi'V^ fall into two categories: 

1. Other semileptonic decays of A{,, which either include additional parti- 
cles (eg: Ah A+7r"'"7r~/i~I/^) or include a higher mass charm baryon 
with subsequent decay into our charm signal (eg: A^ Ac{2593)~^ ii'V^, 
A,(2593)+ ^ A+7) 

2. Baryonic semileptonic decays of Bu,d,s^ which decay into A+ or higher mass 
charm baryon, a proton or a neutron and leptons (eg: B~ —>■ A+p/x~I7^) 

Other Aft Semileptonic Decays 

We have learned from Section 7.2.1 that these physics backgrounds should have 
the same or similar Feynman diagrams as A;, A+/x~I7^. The observation of 
Ac(2593)+ with spin | and Ac(2625)+ with spin | [89] [90] indicates the existence 

of decays A;, —>■ Ac(2593)"'"/i~I7^ and A;, — Ac(2625)"''/i~!7^. Leibovich and Stew- 
art [91] predict a range of the relative decay widths of A5 to these two excited Ac 
decays to the inclusive semileptonic decay. A;, —>■ A+/i anything: 

0.083 < r(A> - A.(2593) V-l. ) ^ ^^^^ 

T{Ab ^ Ajfi anything) 

0.079 <r(f;-W)V-T^^)< 0.166. 
r(Aft — i> A+/i anything) 

However, the wide range gives a systematic uncertainty as large as 100%. In addi- 
tion, the following decays have similar Feynman diagrams as A^ — > Ac(2593)"'"/x~I7^ 
and Ab Ac(2625)+/i-I7^: 

Ab Atffi-u^, Ab SOtt+^u-F^, 
Ab A+n^n^fi-u^, Ab S+ttV"'^/., 
Ab A+7r+7r~^~l7^, A^ E++7r~^~l7^. 

Corresponding decays in the r channel can also produce contamination in our 
semileptonic signal, as seen in Section 7.2.1. 

In order to reduce the systematic uncertainties due to the branching ratios of 
these backgrounds, Litvintsev reconstructed the following decays [86]: 
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Table 7.3: Physics backgrounds in B ^ D*^ jj, 



Mode 


BR (%) 


e ratio 


Norm 






0.276± 0.021 


1 


- 


106 ±11 


D , n*+i/— V 

JD — ^ U jl A 








1 nc;o 4- QQ 




5.44 ± 0.23 


0.447 ± 0.006 


1.000 




JD— , 7~>0,, — T7 

— > U TT 


U.OD ± U.iD 
00. 0( ± (^C).C)C)J 


U.ooD ± U.UUo 


U.Uoo 


Oi ± 




U.OO ± (^U.ODJ 

33.33 ± (1.67) 


n Q/i o -I- n nns 
u.o4y ± U.UUo 


n no 7 
U.UZ / 


OK -L- Q 


B- D'^ji-V^ 

^ D*+TT- 


0.37 ± (0.37) 
66.67 ± (3.33) 


0.351 ± 0.008 


0.036 


33 ± 4 




0.37 ± (0.37) 
33.33 ± (1.67) 


0.336 ± 0.008 


0.017 


16 ± 2 


B- ^ D*+7r-^-V^ 


0.20 ± (0.20) 


0.242 ± 0.007 


0.020 


19 ± 2 


B^ D*+Tr^l2-u^ 


0.10 ± (0.10) 


0.239 ± 0.006 


0.010 


9± 1 


^ D*+T-Vr 


1.60 ± (1.60) 
17.36 ± 0.06 


0.136 ± 0.005 


0.016 


15 ± 2 



Table 7.4: Physics backgrounds in 5 — > D'^fj, i/^ 



Mode 


BR (^X) 


f ratio 


Norm 


^rr( lit 


^ D+n- 


0.276± 0.025 


1.000 




579 ± 30 


B D+ii-X 








4720 ±100 


B" ^ D+ii-V^ 


2.14 ± 0.20 


0.455 ± 0.004 


1.000 




^ D*+^i-V^ 


5.44 ± 0.23 
32.30 ± 0.64 


0.372 ± 0.005 


0.671 


1373 ± 71 


B" D+7rV-17^ 


0.30 ± (0.30) 


0.165 ± 0.004 


0.051 


104 ± 5 


B^ D^n^ fi'Vij, 


0.60 ± (0.60) 


0.165 ± 0.004 


0.102 


208 ± 11 


B- ^ Z^?/i-z7^ 

^ D*+7i- 


0.56 ± 0.16 
66.67 ± (3.33) 
32.30 ± 0.64 


0.278 ± 0.005 


0.034 


70 ± 4 


B- ^ L>?/i-I7^ 

D*+7r- 


0.37 ± (0.37) 
66.67 ± (3.33) 
32.30 ± 0.64 


0.273 ± 0.005 


0.022 


46± 3 


B" D+T-Vr 


0.70 ± (0.70) 
17.36 ± 0.06 


0.100 ± 0.004 


0.013 


26 ± 1 


Bs ^ D+K^fi-V^ 


0.30 ± (0.30) 


0.137 ± 0.005 


0.011 


23 ± 1 
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Figure 7.4: Invariant mass of D~^jj,~ for the signal and physics backgrounds from 
semileptonic decays [87]. The top and bottom histograms are the same plot on 
a linear and log scale, respectively. Note that the backgrounds are concentrated 
in the low mass region. The signal to background ratio is larger at the higher 
mass region. 
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Ab A,(2593)+Ai-I7^, Ab Ae(2625)+/i"i^M> 
Ab E°7r+//-i7^, Ab E++7r-/x-I7^, 

using the data processed with the Production executable, version 5.3.1, and 

compressed into the secondary datasets xbhdOd and xbhdOe. The cuts apphed in 
Litvintsev's analysis are similar to the ones for this analysis. After taking into 
account the efficiency difference between the reconstructed backgrounds and the 
Ab — > A+/x~I7^, we extract the relative branching ratios of these backgrounds to 
the Ab A+ji-V^: 



B{Ab ^ Ae(2593)+/i-l7^) 

B{Ab A+fi-V^) 
B{Ab Ae(2625)V-77^) 
B{Ab A+f^-V^) 
B{Ab ^ Sg7r+/x-I7^) 
B{Ab ^ A+/.-l7^) 
B{Ab ^ S++7r-//-l7^) 
B{Ab ^ Atii-v,) 

Assuming the isospin symmetry, we can infer 

B{Ab ^ S+ttV-^^m 



(4.7± 1.6) X 10-^ (7.23) 

(7.9± 1.5) X 10"^ (7.24) 

(4.2± 1.6) X 10-^ (7.25) 

(4.2 ± 1.6) X 10^1 (7.26) 



B{Ab ^ A+ii-V^) 



= (4.2 ± 1.6) X 10-^ 



However, we need the input of B{Ab A+/i~I7^) to obtain the absolute branch- 
ing fractions of the decays listed above. Assuming heavy quark symmetry, we 
expect the semileptonic decay width for all B hadrons are the same, i.e. Fy^™' = 
■p^mi ^ r^mi = r^™'. Then, the semileptonic branching ratios of the B hadrons, 
psemyptotai^ Only Vary due to the lifetime difference which result in a difference in 
ptotai gjn(.g ^jig decays to a spin-| Ac, we expect contributions from both S 

and P wave amplitudes. Taking a weighted average of B{B D+fiu + D*+iiu) 

and B{B~ D^fiu-l > D*^/ju), where the branching ratios of the and 

(D*+ and D*°) final states correspond to the S (P) wave amplitudes, we obtain 
7.83 ± 0.26 %. We further scale the number by the ratio of hfetimes, J^^^ = 0.80, 
and estimate B{Ab A+n'V^)^ (6.26 ± 0.21)%. For the B{Ab A+t-V^), we 
multiply B{Ab ^ AtiJi-V^) by the ratio of phase space p';,'sltl^%^i^-l^}) = 

Adding the branching fractions of A^ to A^, Ac(2596), Ae(2625''), E0'+'++ 
semileptonic decays in the and r channels after correcting for the B{t — > 
IiV^Ut^) = 17.36 ± 0.06%, we get 8.2%. The branching ratio of the inclusive A^ 
semileptonic decays in the 2004 PDG is: 

B{Ab A+fx anything) = 9.2 ± 2.1%. 

We fill the difference, 1.0%, with the following decays: 
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where the branching fraction of A^ — > A+tt+tt"//"!/^ is estimated to be twice of 
A5 ^ A+7r'^7r"/ir77^ based on the isospin invariance. A more detailed description 
about the estimate of these A5 semileptonic decays can be foimd in Tesarek [88]. 
Table 7.5 summarizes the physics background from the A^ semileptonic decays 
discussed above and their relative efficiencies to the hadronic signal. The dom- 
inant backgrounds are A;, — > Ac(2593)+//""I7^ and A;, — > Ac(2625)+//~i7^. Total 
physics background in Table 7.5 is about 9.2% of the B — > K'^jjL'X events in the 
data. 



B Meson Baryonic Semileptonic Decays 

While there are branching ratio measurements of the B baryonic hadronic decay, 
eg: B —>■ A+p7r+7r~, there is only an upper limit for the semileptonic decay of 
Bu/ Bfi mixture: 

B{ByB+ Kpeue) < 0.15%. 

In order to obtain a list of B meson baryonic semileptonic decays that could 
contribute to our signal, we make use of the predictions from the PYTHIA. We 
generate 5+, and Bg mesons using PYTHIA. We force the mesons to decay 
semileptonically and let PYTHIA handle the fragmentation. We count the num- 
ber of events for each specific B meson to Aj!" decay, and find the maximum 
contributions are from the decays B^ A+pyU~z7^ and 5^ — * A^n^~z7^. We 
assume the upper limit for the muon-ncutron or muon-proton final state should 
be the same as those decays with a proton and electron in the final state (see 
above) . The value of this upper limit is then used for each branching fraction of 
B" — > A^p/i~I7^ and B — > A^ri//""!/^. 

We obtain the efficiencies for these two decays from the MC. Since we find 
the Pt spectra of B mesons and A;, are quite different, it is least ambiguous to 
calculate the quantity: 

where the production cross section ratio is the p(4) in Equation 7.18. We use a low 
Pt threshold (Pt > 4 GeV/c) because we wish to accurately assess the acceptance 
of the B hadron after applying the reconstruction requirements. Specifically 
we are concerned about the case where the neutrino is emitted in the direction 
opposite to the direction that B hadron is traveling. This case increases the P^ 
of the remaining daughters and may make their total P^ greater than the P^ 
of the B hadron. Therefore, the denominator of the efficiency is the number of 
events with B mesons or A;, Pp > 4 GeV/c and rapidity < 2.0. The numerator 
is the number of events which pass all the trigger and analysis cuts. Table 7.6 
summarizes the B meson to A+/i~ backgrounds. The contribution of S° and B'^ 
in the B A+/x~X events is about 0.4% each. 
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Table 7.5: Physics backgrounds in — > A+/x from other A^ semileptonic 
decays. 



Mode 


BR (%) 


e ratio 


Norm 


Af event 


Afc ^ A+TT- 


0.41 ± 0.21 


1 




179 ± 19 


S ^ A+^-X 








1237 ± 97 


Afo A+^"z7^ 


6.26 ± 0.21 


0.300 ± 0.004 


1 




Ab ^ Ac(2593) V'^/. 
E++7r- 
^ A+7r+ 

^A+vr- 

^A+vrO 

■-^ A+TT+TT" 

^ A+7r°7r° 
^A+7 


0.30 ± 0.13 
24 ± 7 

100 
24 ± 7 

100 

24 ± (1.2) 
100 
18 ± 10 
9 ± (0.45) 
1 ± (0.05) 


0.196 ± 0.003 


0.031 


26 ± 3 


Ab ^ Ac{2625)+fi~V^ 
^ A+7r+7r~ 
A+7r°7r° 

--A+7 


0.49 ± 0.17 
66 ± (3.3) 
33 ± (1.7) 

1 ± (0.05) 


0.191 ± 0.003 


0.050 


42 ± 4 


A AH- -f , —rr-: 

Ab A+/> i/^ 


0.25 ± (0.z5j 


0.023 ± 0.002 


n nno 

O.OOo 


2.6 ± 0.0 


A5 A+7r+7r-^-77^ 


0.50 ± (0.50) 


0.032 ± 0.002 


0.009 


7± 1 


Ab ^ A+7r°n^fi-V^ 


0.25 ± (0.25) 


0.033 ± 0.002 


0.004 


3.6 ±0.4 


Ab E°7r+/i"z7^ 
A+vr- 


0.26 ± 0.13 
100 


0.081 ± 0.004 


0.011 


10 ± 1 


Ab S+7r°/i"I7^ 
^ A+ttO 


0.26± 0.13 
100 


0.072 ± 0.004 


0.010 


8± 1 


Ab E++7r~/x~I7^ 
^ A+7r+ 


0.26± 0.13 
100 


0.077 ± 0.004 


0.011 


9± 1 


Ab ^ A+T-I7^ 


1.74 ± (1.74) 
17.36 ± 0.06 


0.040 ± 0.003 


0.007 


5± 1 
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Table 7.6: Physics backgrounds in — > A+yu from B mesons. 



Mode 


BR (%) 


Relative e 


Norm 


■^event 


Aft ^ A+TT- 


0.41 ± 0.21 


1 




179 ± 19 


B ^ A+/i-X 








1237 ± 97 


Aft A+//-I7^ 


6.26 ± 0.21 


0.265 ± 0.004 


1 




-B~ ^ A+p;U-z7^ 


0.15 ± (0.15) 


0.045 ± 0.002 


0.006 


4.7 ± 0.5 




0.15 ± (0.15) 


0.044 ± 0.002 


0.006 


4.6 ± 0.5 



7.3 Fake Muons 

Another source of background originates from a charm hadron together with a 

hadron track (vr, proton) misidentified as a muon. A hadron is misidentified as 
a muon when it has higher energy and punches through the hadron calorimeter, 
or when it decays into a muon before being stopped in the hadron calorimeter via 
decays like tt"*" — > /I'^i'/x, f^~^^ij. or — > tt^/x+i/^^. Physics processes that 

generate these hadrons are direct production, inelastic coUisions with the detector 
material, fragmentation or the decays of charm and B hadrons. Fragmentation is 
the process by which a 6 or c quark combines with additional quarks and gluons to 
form a qq or qqq bound state. Fake muons from the first three categories tend to 
have a softer Pt spectrum than the real muons from B decays. A tighter Pt cut 
on the muon candidate largely removes these backgrounds. Fake muons from the 
charm hadrons which are produced at the primary vertex are also suppressed. 
For the reason that we require the muon candidate should be matched to an 
SVT track with a do greater than 120 yum, while fake muons from the promptly 
produced charm tend to have smaller impact parameter. Also for the reason that 
we require the charm hadron and the muon candidate to form a vertex displaced 
from the beam line and make a strict requirement on the pseudo cr. 

Mb 

= P.(charm + ;.) ' ^^'^^^ 

Here (charm + /x) is the total transverse momentum of charm hadron and the 
muon. 

Therefore, our principle source of fake muons comes from two types of B 
hadron decays: 



B — > D*+ (D"*", A+)Xhad anything: hadronic decays of any B hadrons, 
where Xhad is tt, K or proton which fakes the muon. 

B D*^ {D'^: A+)Xhad I T^i anything: semileptonic B decay into a charm, 
a hadron track Xhad) and any leptons (e, ji, r). The muon is not recon- 
structed but Xhad fakes the muon. 
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In this section, we estimate the fake muon contamination for our three signals: 
— > D*~^iJ,~V^, ^ — > D^ii~V^ and A;, — > h.'^ijrVfj^. Note the charge conjugates 
of the modes hsted are also included. 

7.3.1 Background Estimate 

We use two methods to estimate the amount of contamination from fake muons in 
our semileptonic signal. Each method uses a different way to obtain the number 
of (hadron track, charm hadron) candidates in our data. Both methods apply the 
previous CDF measurements of the probabilities for a real pion, kaon and proton 
to be misidentified as a muon. These measurements are performed by Ashman- 
skas and Harr [92] using a pion and kaon sample from the D*^ D^'tt"'" decays, 
where K~7r~^, and by Litvintsev [93] using a proton sample from the 

A — > p+TT" decays. The fake probability {Vk, Vp) is defined as the number of 
pions (kaons or protons) that pass the following muon identification cuts divided 
by the total number of pions (kaons or protons) inside the fiducial volume of CDF 
Central Muon Detector (CMU) and matched to an SVT track. 

• The track is fiducial to the CMU and matched to an SVT track 

• The track is associated with hits in the CMU 

• The matching between the track and the hits in the CMU is less than 9. 

Figure 7.5 shows the P^, Vk measured in sixteen and Vp measured in twelve 
transverse momentum bins for positive and negative charged tracks, separately. 

Method I 

The first method uses data to obtain the number of (hadron track,charm hadron) 
candidates, then Monte Carlo to determine the ratio of pions, kaons and protons 
in the hadron tracks. We run the same signal reconstruction program on the 
secondary datasets hbotOh and hbotli. We do not use the skimmed tertiary 
datasets (see Yu [64] ) as the samples are biased by requiring at least one track in 
the event matched to a muon stub in the muon detector. We look for a charged 
track which fails the muon identification cuts (TRK-^"''). TRK-^"^' and a charm 
hadron should form a displaced vertex and pass the same analysis cuts we apply 
to the signal. Each event is weighted with the fake probability (Pavg) according 
to the momentum and the charge of TRK-^"*'. We then fit the weighted charm 
hadron mass distribution, i.e. Mk-ktt — Mxn, Mk-ktt and MpXn, using the same 
functions as described in Section 5.1, to obtain the signal contamination from the 
fake muons. Since an event- weighted likelihood fit will not give a proper error for 
the yield, a binned fit is performed. "Pavg is a weighted average of pion, kaon 
and proton fake probabihty {V-jr, Vk, Vp). The weight Ri is determined by the 
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fraction of pions, kaons and protons in the B — > D*~^ {D^^ ■^^)-'^had ^ anything 
and B — > !?*+ A+)Xhad anything MC after analysis cuts: 

= R^V^ + RkVk + RpVp, (7.29) 

where 

and i is TT, or proton. The Monte Carlo is generated as described in Section 6.1. 
Decays of B^, B^, Bg and Af, are generated separately and decay tables include all 
the possible decays. Each kind of B hadron gives different Ri and is weighted with 
the product of the production fractions, total branching ratios and the number of 
generated events. Table 7.7 summarizes the pion to kaon ratio and the number of 
fake muon candidates before and after weighting the events in our three different 
signals. See Figure 7.6 for the weighted mass distribution of the A;, mode. 

The uncertainties of the number of fake muon candidates come from three 
sources: 1. the uncertainty from the binned fit, 2. the uncertainties on 
the pion, kaon and proton fractions due to the finite Monte Carlo sample size, 
the uncertainties of the branching ratios and production fractions, and 3. the 
uncertainty on the measured fake probability. For the last source, we vary the 
fake rate in each momentum bin ± 1 sigma, independently. We then add the 
systematic shifts in quadrature to get the accumulative uncertainty. The number 
of fake muons using this method is about 3-5% of the inclusive semileptonic 
signals in the data. 



Method II 

The second method relies on the MC and the external input of the branch- 
ing ratios from the Particle Data Group (PDG) [8]. MC is run through the 
same reconstruction program for the data. Then we apply the same cuts as 
signal reconstruction and obtain the trigger and reconstruction efficiencies of 
B D*+ {D+, A+)Xhad anything and B D*+ {D+, A+)Xhad I t^i anything de- 
cays. We weight the MC events which pass the analysis cuts with the fake proba- 
bility according to the momentum, the charge and the particle type of the track, 
Xhad- The particle identification of the track, Xhad, is obtained by matching the 
hits on the reconstructed track with those on the input simulated track. Together 
with the efficiency of hadronic mode, branching ratios of our hadronic signals, 
B ^ D*+ (D+,A+)Xhad anything and B D*+ (D+,A+)Xhad / Vi anything 
from the PDG, we normalize the background to the observed number of hadronic 
signals in the data, 

ATfa,,^ _ B(B ^ D*+ A+)X) • efajce^ 

Nhad fd ,baj:yon • ^had^had 

Equations 7.31- 7.32 use as an example to show how we derive the B — >■ 
D*+ A+)Xhad anything and B D*+ (£>+, A+)Xhad I i^i anything branch- 
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ing ratios from the existing information in the PDG. 



— 7 — 

B{B L>+Xhad li^i anything) = - • B{B L>+7r+// anything) 

7 

■ -■ ^ L>+7r-/x anything), (7.31) 

B(B L'+Xhad anything) = B(B D+anything) 

7 _ 

- - • B{B D^ji anything), (7.32) 

where the factor, |, comes from the fact that the branching ratios of muon and 
electron channels are equal and the branching ratio of the tau channel is scaled 
down by the ratio of the phase space, ~ \. Therefore, we have to scale up the 
branching ratio of the muon channel byl + l + | = |to get the total branching 
ratio of all the lepton channels. 

Table 7.8 summarizes the parameters used to calculate the number of fake 
muon events, where the decay B — > D*^ (L>+,A+)Xhad anything is denoted as 
mode "1" and B £)*+ (£)+, A+)Xhad I T^i anything is denoted as mode "2" in 
the tabic. The uncertainties on the number of fake muons originate from: the 
uncertainty on the hadronic yield, the relative efficiency, the uncertainty on the 
fake rate and the relative branching ratios. The dominant uncertainty is from the 
relative branching ratios. The number of fake muon backgrounds from method I 
is consistent with the result using method II. We use the results of method I in 
the calculation of our final result of the relative branching ratios. In general, the 
fraction of fake muons is about 5% of the total semileptonic yield in the data. 



Like-sign Combination 

Note that wc do not use the like sign combination (i.e. the charm hadron and the 
muon have the same sign of charges) to estimate the fake muon background for 
two reasons: First, two different B hadrons from the hb in the event can produce 
a real muon and a real charm of the same charge sign when the B hadrons in the 
event have opposite flavors and one B hadron decays semileptonically. Second, 
the two track trigger, used for this analysis, requires a pair of tracks with opposite 
charges. The trigger requirement greatly reduces the number of like-sign (wrong- 
sign) candidates and introduces large statistical errors for the number of fake 
muons. 



Fcike muons from hh and cc 

One type of fake muons is not included in the previous subsections. These fake 
muons stem from hb, cc to two B or charm hadrons then decay into a charm 
signal, a hadron track misidcntified as a muon and other missing particles. A 
study at the generator level for the B D^/j^V^ mode is done using the bb 
and cc PYTHIA [71] Monte Carlo datasets as described in Section 7.4.2. We apply 
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Table 7.7: Parameters for the number of fake muons: Method I. 



N before weighting 

Rk 
Rp 

-^fake 






Aft ^ A+ij.-Vi^ 


2953 ± 57 
0.937 ± 0.009 
0.063 it 0.009 


15343 ± 303 
0.909 ± 0.005 
0.091 ± 0.005 


3560 ± 198 
0.71 ± 0.16 
0.05 ± 0.08 
0.24 ± 0.16 


44 ± 3 


230 ± 19 


40 ± 9 



Table 7.8: Parameters for the number of fake muons: Method II. 









Aft A+/i 


^had % 


0.276 ± 0.021 


0.276 ± 0.025 


0.41 ±0.19 


Bi % 


10.9 ± 2.1 


17.7 ± 2.4 


4.8 ± 3.0 


Chad 


0.0038 ± 0.0004 


0.0022 ± 0.0002 


0.0029 ± 0.0003 


B2 % 


1.3 ± 0.3 


1.8 ± 0.6 


< 1.23 


^had 


0.0005 ± 0.0002 


0.0010 ± 0.0002 


0.0002 ± 0.0001 




106 ±11 


579 ±30 


179 ±19 


A'fake /i 


45 ± 11 


220 ± 41 


28 ± 34 
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Figure 7.5: The probability for a pion, kaon or proton being misidentified as a 
muon in bins of transverse momentum (Pt) from the measurements by Ashman- 
skas, Harr [92] and Litvintsev [93]. From the top left to the bottom right are 
7r+, TT^, K~^, K~, p and p fake probabilities. 
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Figure 7.6: Fit of B ^ A+/ifake yield after weighting the charged track which fails 
the muon ID cut with an average muon fake probability. There 40 ± 6 events in 
the peak. Fit x^/NDF = 55.6/39, probability = 4.1%. A sideband subtraction 
yields 44 ± 25 events in the signal peak. 



Table 7.9: Fake muons from bb and cc. 





bb 


cc 


N 

gen 

Real muon l^pass 
Fake muon Np^ss 


43454949 
15 
1.8 


89718181 
35 
0.4 



analysis-like cuts on the Monte Carlo. We weight the events that pass the cuts 
with the muon fake probability according to the "muon" candidate momentum, 
charge and the true particle identification: a kaon, a pion or a proton. Then we 
compare the number of weighted events with the number of charm hadron and 
real muon combinations, i.e, a bb and cc background as described in Section 7.4. 
We find that fake muons from bb and cc is about 10% of the bb and cc background 
with real muons. See Table 7.9. From Section 7.4, we show that the bb and cc 
background with real muons is at the 1% level. Therefore, we conclude that 
bb and cc background with fake muons is about or less than 0.1% and can be 
ignored. 
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Figure 7.7: Charm hadrons from cc with small (left) and big (right) A0. In 
the left figure, the muon from the semileptonic decay of and the forms 
a secondary vertex and fake our B ^'V^ signal. In the right figure, A0 

between two charm hadrons is too big and the daughters can not form a secondary 
vertex. 



7.4 hh and cc Backgrounds 

When the azimuthal angle (A0) between hb or cc quark pair is small, daughters of 
two heavy flavor hadrons from the fragmentation of bb or cc appear to come from 
the same decay vertex, see Figure 7.7 . Here A0 is defined as the opening angle in 
the plane perpendicular to the proton and antiproton beam axis. If one hadron 
decays semileptonically, and the other hadron decays into a charm final state, 
such as D*^ — >■ D^TT'^, — >■ K~7T~^7r~^, and — ^ pK~7r^, the muon from the 
semileptonic decay, together with the charm, may fake our semileptonic signal. 
Production mechanisms and an estimate of the amount of bb and cc backgrounds 
are discussed below. 

7.4.1 bb and cc Production Mechanism 

In pp collisions, the b or c quarks may be single or pair produced by the elec- 
troweak and the strong (QCD) processes. The b or c quark production cross- 
section for the electroweak process a ■ B{pp ^ W ^ be) is around 0.01 fib and 
is derived from the CDF measurement of the inclusive W cross-section by Halki- 
adakis, et al. [94]. The bb and cc production cross-sections for the QCD process 
are around 50 and 200fib respectively from the PYTHIA [71] Monte Carlo, when the 
total transverse momenta of the hard scattering, i.e. the part of the interaction 
with the largest momentum scale, is greater than 5 GeV/c and at least one b or 
c quark has Pt > 4.0 GeV/c, pseudo-rapidity r] < 1.5. The bb and cc production 
rates from the electroweak process are about five thousand times smaller than 
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FigTirc 7.8: Representative lowest order Feynman diagrams (without loops or 
radiative corrections) of quark annihilation, gluon fusion, flavor excitation and 
gluon splitting. Details of these processes may be found in Lannon [95] . 



the QCD processes. Therefore, only pair production by the QCD processes are 
discussed here. 

Figure 7.8 shows the leading and next-to-leading order Feynman diagrams for 
bb{cc) production by the QCD processes from Lannon [95]. The QCD process 
that contributes the production at leading order is flavor creation; which includes 
quark anti-quark annihilation {qq — > bb or cc) and gluon fusion (gg — > bb or cc). 
The distribution of the azimuthal angle (A0) between two produced b(c) quarks 
peaks at 180 degrees. The reason is that the q and q (or g and g) come from the 
proton and the anti-proton separately. The initial total momenta of the gluon qq 
pair is zero. The bb (cc) pair are produced back-to-back to balance the momentum 
in the final state. 

The next-to-leading order (NLO) processes, flavor excitation and gluon split- 
ting, contribute at the same level as the flavor creation [95]. Flavor excitation 
refers to the following process: The gluons within one of the beam particles in 
the initial state split into a bb {cc) pair. One of b(c) quarks is scattered out of the 
initial state into the flnal state by a gluon or a light quark from the other beam 
particle. The other b(c) quark is not involved in the hard scattering process. 
The A0 of bb {cc) from flavor excitation is more evenly distributed than flavor 
creation. 

Gluon splitting refers to the process when no b(c) quarks are involved in the 
hard scattering. One gluon in the flnal state splits into bb {cc) pair. If the gluon 
is soft, A(f) will be a flat distribution. If the gluon is hard, the daughters of the 
gluon, bb {cc) tend to move co-linearly and have small A0. Figure 7.9 shows the 
azimuthal angle distribution between two b quarks from the study of Field [96] 
using PYTHIA CTEQ4L prediction. 
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Figure 7.9: A0 between b and b from the study of Field [96] using PYTHIA 
CTEQ4L. Both b quarks have \y\ < 1. In the top plot, one b quark has Pt > 5 
GeV/c and the other b quark does not have any Pt cut. In the bottom plot, 
both b quarks have Pt > 5 GeV/c. 
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7.4.2 Background Estimate 



The amount of bb and cc background is normalized to the number of events 
observed in the hadronic modes in the the data, 

^bb,cc _ ^bb,cc ' f^^j^j gg.^ 

Here, i represents the species of b(c) hadrons and j represents the decay modes 
which could contribute to bb and cc backgrounds. /* stands for the production 
fraction ratio for species i. Bj and are the branching ratio and the efficiency 
of j*^ decay mode. The following subsections detail the methods to estimate 
(^bb cc ' f^^j^j '^B",Ab ' -^hadChad in Equation 7.33. We do not use detector 

and trigger simulations to obtain the efficiencies for the following reasons: First, 
detector and trigger simulations are time and CPU intensive. Second, we will 
find the contribution of this background is quite small compared with the other 
backgrounds. Third, we care about the efficiency ratio of the background to 
the signal, not the absolute efficiency. Our studies show that generator level 
Monte Carlo gives a good approximation. For instance, the relative efficiency 
f.^^''^ +° , is 3.31 ± 0.05 from the full detector simulation and 3.23 dz 0.01 
from the generator level simulation. The difference is only about 2.5%. Similar 
results are obtained from the relative efficiencies of our other signals. 



Background: (j;,^ ^- • Ei Ej P^j^j 

Our estimate of cr^^^c ' Si Y^j P^j^j relies heavily on the Monte Carlo. We use 
PYTHIA version 6.2 [71] and to generate bb and cc events, we include the QCD 
processes mentioned in Section 7.4.1: flavor creation, flavor excitation and gluon 
splitting (MSEL=1). We further require the Pt of the hard scattering be greater 
than 5 GeV/c. Events with b quarks Pt greater than 4.0 GeV/c and pseudo- 
rapidity less than 1.5 are collected into the nbot90 sample. Events with c quarks 
which satisfy the same kinematic cuts are collected into the nbotaO sample. Note 
that nbot90 and nbotaO have small overlap when both b and c quarks are pro- 
duced and arc above the Pt and pseudo-rapidity thresholds. Details of nbot90 
and nbotaO datasets could be found in [97]. PYTHIA bb and cc cross-sections are 
used for (j;,^^^^. 

The product of the efficiency, branching ratio and production fraction from all 
the modes, J2i J2j P^j^j, is obtained using the following steps: First, we identify 
the hb (cc) production mechanism to which we are most sensitive. In this pass, 
we only study the background that forms a D~^fi~ signature since has longer 
hfetime than D*+, A+ and W D+ 11 suffers larger cc background contami- 
nation compared to the other two modes. We need to achieve higher accuracy for 
the estimate of P^j^j- Consequently, for the bb background, we re-decay 

nbotQO sample ten times, i.e. we rc-usc the same kinematic distribution of the 
parent hadrons from nbot90 ten times but decay the hadrons with independent 



125 



Table 7.10: Summary of bb and cc production mechanisms and our relative sen- 
sitivity for reconstructing the event in our semileptonic sample. 





7 7 11 1 

bb background 


cc background 




219093011 


21996889 


N 


75 


62 


^gluon 


70 


57 


^excitation 


5 


5 


^creation 








igluon{%) 


93 ± 3 


92 ± 3 




7 ± 3 


8 ± 4 


^creation{%) 









random numbers and force the decay K~7r~^7r~^. For the cc background, 

we force the decay of — >• K~7i~^n'^ and require that all the negative charged 
charm hadron decay semileptonically. 

Then, a generator level two track trigger filter (SvtFilter) is applied. We 
further identify any combination of and a muon which passes the generator- 
level, analysis-like cuts found in Table 7.11. In order to avoid double counting 
(count cc as bb background in nbot90 and bb as cc background in nbotaO) due 
to the overlap of nbot90 and nbotaO samples, the ancestors of the muon and 
charm hadrons are retrieved by tracing the true information from the generator. 
If both muon and charm hadron come from the same B hadron, the combination 
is rejected. If the muon and charm hadron come from different B hadrons, the 
combination is categorized into bb background, otherwise, the combination is 
categorized as a cc background. We find that for both bb and cc, more than 
90% of the events that pass the cuts are from gluon splitting. Therefore, wc arc 
most sensitive to the "gluon splitting" mechanism. Table 7.10 summarizes the 
background contributions from different production processes. 

Second, we filter the gluon splitting events and re-decay the b(c) hadrons 
in nbot90 and nbotaO ten times with the procedure described above. For the 
bb background estimate, we let all the b hadrons and negative charged charm 
hadrons decay freely, but force the decays of the positive charged charm hadrons 
in two ways: 

, D*+ D\+, K-7r+ for the background oiB°^ D^+ii'V^ 

• K~'K'^'K'^ and A+ — > pK~7r'^ for the background of — > D~^ijrV^ 

and — > A+//~i7^ 

Then the SvtFilter and the cuts hsted in Table 7.11 are applied. We divide 

the number of reconstructed events by the number of generated events and get 
HiHj r^j^j- Table 7.12 lists the parameters for the bb background. 
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For the cc background estimate, we force the decays of both positive and 
negative charged charm hadrons. The positive charged charm hadrons are forced 
to decay into the modes hsted above. The negative charged charm hadrons are 
forced to decay into semileptonic modes individually for D~ , , D~ and A~. 
As the semileptonic decay modes of these four charm hadrons are all different, 
we separate the events into four classes denoted by the parent charm particles. 
After applying SvtFilter and the cuts listed in Table 7.11, we obtain P^j^j 
for each class. Then we muhiply the semileptonic branching ratios for each kind 
of charm hadron with its f^^j^j ^^^^ sum them up to get the total amount of 
Y,j r^j^j for the cc background. 



i j i j 

+ E E rBjej{A-)B{A- ^ piX). (7.34) 



Table 7.13 lists the parameters for the cc background. 

In both bb and cc background estimates, since we force the positive charged 
charm hadron to decay into the same final state as our charm signals, we have 
to multiply the final result by two to include the contribution from both charge 
states. Table 7.15 lists N^^ and Ncc in our three different signals after multiplying 
the ratio in Equation 7.33 with the observed number of events in the hadronic 
signals. 

Hadronic signal: ctbo a, ■ -BhadChad 

In order to normalize the background to the observed number of events in the 
hadronic mode, the or A;, production cross-section, the efficiency and the 
branching ratio of the hadronic signal, have to come from external input or must 
be calculated using MC (see Equation 7.33). We obtain ctbo by multiplying the 
previous CDF aB+ measurement by Keaffaber, et al. [98] with the production 
fraction ratios, fd/fu, from the 2004 PDG [8]. decay branching ratios are 
also obtained from the PDG. The product of a\^^ and B{Ah — > A+tt^) is obtained 
by multiplying the CDF measurement of "'H^^'r'^^-^^'^^^~*^'' ) rggi 

with the cr^o we derive and the PDG B{B^ D^n^). 

Since we reconstruct both b and anti-b hadrons in the data, we should multiply 
the measured cross-section by two. For the efficiencies, we use the MC to generate 
and decay B hadrons into our signals as described in Section 6.1. The CDF aB+ 
measurement is restricted to the B~^ with greater than 6 GeV/c and rapidity 
{y) less than 1.0. Therefore, the denominator of the efficiency is the number of 
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Table 7.11: Generator- level analysis- like cuts for bb and cc background study. 



Parameter 


Cut Value 


Pt of all tracks 


> 0.5 GeV/c 


Pt of ijl {-Kb) 


> 2.0 GeV/c 


rj of all tracks 


< 1.2 


77 of ^ (ttb) 


< 0.6 


Pt of four tracks 


> 6.0 GeV/c 


Pj of cliariii liadroii 


> 5.0 GeV/c 


CT of four tracks 


> 200 //m (B), > 250 //m (At) 


cr of charm hadron 


> -70 Atm {D*, AJ, > -30 //m {D+) 


3.0 < M^(,)^ < 5.5 GeV/c^ 


3.7 < Ma,^ < 5.7 GeV/c^ 


H (ttb) match to a SVT track 


charm hadron and fi 


(ttb) have opposite charge signs 


2 out of 4 tracks of B candidate pass two track trigger cuts 



events in which the B hadrons have Pt > Q GeV/c and \y\ < 1.0. The numerator 
of the efficiency is the number of events which pass the cuts in Table 7.11 except 
the cut on four track invariant mass. Table 7.14 lists the parameters that are 
used to calculate (Jbo,A(, ■ -Shadehad- 

Semileptonic signal 

Table 7.14 also lists the efficiency and branching ratio of the semileptonic signal 
mode for a comparison with the background N^j^^-. In the case of the A5, we 
lack the external input for the branching ratio of A5 — > A+yU~I7^. Therefore, 
instead of the exclusive mode, we list the branching ratio of the inclusive mode: 
B — > h.'^jjrX from the 2004 PDG as an upper bound. We multiply the ratio p(6) 
from Section 7.1 with the KcafTabcr 0"^+ result to get ctaj, for Pt greater than 
6.0 GcV/c. The efficiency of Aj K'^ji^v^ is listed for a comparison with e^^^^. 
Note that the amount of hb and cc background relative to the signal is around 
1%. The numbers from three different modes should not be compared directly 
without multiplying the branching ratios of the charm decays. 

7.4.3 Comparison of Data and MC Cross Section 

While there are precise measurements of the single charm hadron, B^ and in- 
clusive b hadron cross-sections, there are no accurate measurements of the total 
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Table 7.12: Parameters used for bb background estimate. 



PYTHIA a,-^ ii^h) 


49.6 






bb Djj 


bb KcjJi 


N 

'gen 


221606748 


221619610 


221619610 


N 

^ 'pass 


43 


80 


9 


E^E.-M-e. (10-^) 


1.9 ± 0.3 


3.6 ± 0.4 


0.41 ±0.14 


2-<7,5-E.E, M-6. (10-5 ;,b) 


1.9 ± 0.3 


3.6 ± 0.4 


0.41 ± 0.14 



Table 7.13: Parameters used for cc background estimate. 



PYTHIA Ucc (Afb) 


198.4 


TV 

^ 'gen 


cc D*ii 


cc Dfj, 


cc Ac/i 


720741510 


698988700 


698988700 


TV 

^ 'pass 


214 


396 


7 


D- : B{D- liX) = 14.22 (%) 


117 


205 


4 


:D° : 80" ^ ^xX) = 6.15 (%) 


76 


157 


2 


D- : B{D- nX) = 13.32 (%) 


19 


34 


1 


A,- : BiA: ^ i^X) = 4.5 (%) 


2 








EiEj f%ej {10-') 


3.3 ± 0.6 


6.2 ± 1.0 


0.12 ± 0.05 


2 ■ ■ EiEj f%ej (10-5 //b) 


1.3 ± 0.2 


2.5 ± 0.4 


0.047 ± 0.020 
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Table 7.14: Parameters used to calculate (7^0 ^^, Bhad^had and Bsemi^i 



(7b+ (A*b) 


3.6 ± 0.6 




B D*+X 


B D+X 


Afc K+X 


f If 


1 00 ± 04 


1 00 ± 04 


25 zb 04 


(Jro a iuXy] 


3.6 ± 0.6 


3.6 ± 0.6 


2.2 ± 0.5 






B° D+TT- 


Ai, — > A+TT" 


(TA(,(PT>6.0)B(Ai,^A+7r-) 
(TgO (PT>6.0)B(B°^D+ir-) 






0.82 zb 0.26 




n 07R -I- n noi 


n 97fi -4- r\ r\OK 

U.Z / D zt U.UzO 










n nns + n nn'? 

U.UUO zn U.UUO 


N 

^ 'gen 


/1 9/191 nn 


4Z4Z1UU 


'iQQQQQQfi 
oyyyyyyo 


N 

^ 'pass 




10U400 




Chad {^^ ) 


i.D04 zt U.UUD 


Q n7i^ -I- n nna 

o.U / zt U.UUo 


9 1 HQ -I- n nn9 
z.iuy zt u.uuz 




oo zn o 


fil -1- 1 9 
oi zn iz 


o4 zn io 








Af, -> A+/i"z7^ 


^semi (%) 


5.44 ± 0.23 


2.14 ± 0.20 


9.2 ± 2.1 (%) 


N 

^ 'gen 


4242100 


4242100 


39999996 


TV 
'pass 


32620 


66854 


264484 


^semi (10 ) 


0.769 ± 0.004 


1.576 ± 0.006 


0.661 ± 0.001 


2o-B0^A(,^semiesemi (10"^ /ib) 


300 ± 50 


240 ± 50 


270 ± 90 



Table 7.15: The amount of hh and cc background. 





B- 


D*+X 


B- 


D+X 




-.AtX 


A^had 


106 


± 


11 


579 


± 30 


179 


± 19 


N ■ 

' semi 


1059 


± 


33 


4720 


± 100 


1237 


± 97 




6 


± 


0.6 


34 


± 2 


2.1 


± 0.2 


Ncc 


4 


± 


0.4 


23 


± 1 


0.2 


± 0.03 
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bb and cc cross-section (a^^, acc) at the Tevatron, yet. To understand how well 
PYTHIA predicts o"^^ and a^c, we cross-check indirectly by comparing the "differ- 
ential cross-section" of D^, and inclusive b hadrons in PYTHIA with CDF Run 
I and 11 measurements by Chen [99], Keaffaber [98], and Bishai [41] et ai. We 
count the number of D'^, or b hadrons from nbot90 and nbotaO in bins of 
Pt{D^), Pt{B'^) and Pt{J/iP). The bin width and the Pt ranges are the same 
as Chen, Keaffaber and Bishai analyses. We divide the number of hadrons in 
each Pt bin by the total number of generated events. Then we multiply Pythia 
assumed cx;,^ and (see Tables 7.12- 7.13) to get the cross section of hadrons 
in each Pt bin. We further divide the number by the bin width to obtain the 
"differential cross-section" . The agreements between Monte Carlo and data cross- 
sections are generally within 10% for charm hadrons and 40% for B hadrons (see 
Figure 7.10). 

Besides the total cross-section of bb and cc, the ratio of gluon splitting relative 
to the other two processes, flavor creation and flavor excitation, also affects the 
amount of bb and cc backgrounds. Previous CDF Run I measurement of bb 
azimuthal production correlations by Lannon [100] concludes that Pythia gives 
reasonable prediction of the relative bb production rates from the three processes. 
However, due to the lack of measurements of the cc relative production rates, we 
do not yet have a comparison of the fraction of cc gluon splitting between Monte 
Carlo and data. Therefore, we assign 100% uncertainty when calculating the 
systematic errors for the estimate of bb and cc backgroTinds. As the contribution 
of bb and cc background is at the 1% level, the systematic errors from 100% 
uncertainty is also about 1%. 

7.4.4 Comparison of Data and MC Impact Parameter 

We compare the distribution of the impact parameter of charm hadrons with 
respect to the beam spot in MC and data. Figure 7.11 shows a good agreement 
of the MC with the data. No excess of charm hadrons with small do is found in 
the data. This indicates that the promptly produced charm from cc has negligible 
contribution to the background in the semileptonic B decays, which is consistent 
with our estimate using PYTHIA. 

7.5 Background Summary 

The fraction of each type of background in our semileptonic signal is summarized 
below. The dominant signal contamination is from the physics background. The 
second largest background arises from muon fakes. The smallest background 
source is from bb and cc. 
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- Pyttiia D° cross-section (stat. err only) 

flHDF^e.SIB, prob=23.2% 



2000- 


1.8 

1.6 

1.4 

1.2 

1 

0.8 
0.6 



ii 



15 20 
<Pt(D )> [GeV/c] 



„ 1500 
o 
> 

0) 

O 



^ 1000 - 



500- 




1.2 



1- 

0.8 - 
0.6 - 
0.4 1 



0.2 



- Run I cross-section measurement 

- Pytlnia B* cross-section (stat. err only) 

X^/NDF=30.3/3, prob=0.0001% 



20 25 

<Pt(B*)> [GeV/o] 



> 

CD 

o 



DC 
CO 



1.5 



1.8 
1.6 

1.4 
1.2 

1 

0.8 
0.6 



- Run II inclusive b cross-section measurement 

- Pythia inclusive b cross-section {stat. err only) 

X=^/NDF=40.3/10, prob=0.0015% 



15 20 

<Pj{J/^)> [GeV/c] 



Figure 7.10: D'^ (top left), (top right), and inclusive b (bottom) differential 
cross-sections. The upper plot in each figure shows the differential cross-section 
for data (closed circles) by Chen [99], Keaffaber [98], and Bishai [41] and MC 
(open squares). The lower plot in each figure shows the data to MC ratio. 
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-^Data 

prob:0.1008 
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Figure 7.11: Impact parameters of the charm hadrons, from the top left to the 
bottom are D^, and A^: MC and data comparison. The good agreement of 
the MC with the data indicates that background from the promptly produced 
charm (cc) is negligible. 
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Table 7.16: Summary of the backg: 



;rounds to the semileptonic modes. 





-^bg / -^inc sem 


A%) 


Background Type 


D*fi 






Physics 


15 


40 


9.8 


Muon fakes 


4.3 


4.9 


3.2 


bb and cc 


0.9 


1.2 


0.2 


Total 


20.2 


46.1 


13.8 
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Chapter 8 



Relative Branching Fraction 
Results and Systematics 

Recall that the equation for the ratio of branching fractions is: 

^semi , -^inclusive semi ~ -^physics ~ -^fake/i ~ ^ cc, bb\ ^had /o i\ 

» = ( )X- , 1«-1J 

'-'had -''had tgeuij 

where the yields, iVhad and iVinciusivesemi, are obtained in Chapter 5. In Section 6.3, 
the efficiencies, ehad and egemi are determined and in Chapter 7, the backgrounds to 
our semileptonic signals are estimated. We are now in a position to calculate the 
relative branching fractions. The relative branching fractions with the statistical 
uncertainties only are: 

17.7 ±2.3, 

9.8 ±1.0, 
20.0 ±3.0. 

In the remainder of this chapter, we first discuss and estimate the systematic 
uncertainties. Then, we show the result of the relative branching fractions for each 
mode. Finally, implications of our measurements are discussed and a conclusion 
is given at the end of the chapter. 

8.1 Systematic Uncertainties 
8.1.1 Sources of Systematics 

Systematic uncertainties in our measurements may arise from the difference in 
the semileptonic and hadronic decays, from the lack of knowledge of certain 





D*+7r-) 








D+n-) 


B{K- 





B(Af, ^ A+TT-) 
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backgrounds, and from the uncertainties on the external information. Most of 
the sources are common to all the decay modes. Systematic uncertainties which 
affect only one mode are discussed separately. To simplify the notation, we define 
our measurement of the relative branching fractions as R: 

B ■ 

*-^semi 

and (7r is the systematic uncertainty. We also denote the term, "branching ratio" 
in the text, as BTZ, while the branching ratio of one specific mode is denoted as 
S(mode), eg: B(B° D+tt"). 



Mass Fitting 

• B ^ D*^ X and B A+/i~X: The mass functions are general and 
cover all the possible backgrounds. Because the functions do not involve 
any external BTZ or MC efficiencies, we do not assign systematic uncertainty 
for the mass fitting of these two modes. 

• S — > D'^ji'X: The uncertainties on the Dg decay BTZ can modify the mis- 
reconstructed Dg mass spectrum in S — > D'^fj,~X. In addition, the mean of 
the Gaussian constraint for the amount of Dg background in Equation 5.11, 
Hp, also changes accordingly. We study this effect by varying the following 
numbers in Table 5.2 ± 1 cr independently: the BTZ of — > (j)TT~^, and the 
BTZ of each selected Dg decay relative to B{Df — > ^tt"*"), because these Dg 
BTZ were measured relative to the 07r mode [8] . The corresponding Dg back- 
ground shape and the Hp are re-evaluated for each change of Dg BTZ. The 
changes in the yield are added in quadrature to get the accumulative varia- 
tion. Table 8.1 summarizes the yield variations. Modes which are not listed 
give identical results to the central value. We find total A(iV-Q^^+^_^) — 
33 events and ctr = 0.13. Note that a few Dg decays in Table 5.2 only have 
an upper limit in the PDG and the estimated values in the EvtGen decay 
table are used. We assign 100% uncertainty for these modes. 

• 5° — > D*'^7r~: the composition of the remaining B — > D*'^X background 
can affect the shape of its mass spectrum, and its ratio to the D*p back- 
ground. The latter changes the mean of the Gaussian constraint, ^2, in 
Equation 5.19. We study the systematics by varying the BTZ of B^ — > 
D*'^p~ and the dominant modes in the remaining B — > D*'^X background. 
The change of signal yield from each variation of BTZ is listed in Table 8.2. 
The accumulative yield change is only ^'^'^ events, which is insignificant. 
Therefore, we do not assign systematic uncertainty for the mass fitting of 
this decay mode. 

• —>■ D~^Tr~: the systematic uncertainties come from three sources: the 
normalizations of the Cabibbo suppressed decay, Bg and backgrounds. 
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Table 8.1: 5 — > D^jj, X yield change due to the variation of Dg BTZ. 



D+ K+K-7r+ ± 23 

Df^(f)K+ ± 2 

D+^ri7r+ ± 2 

D+^ri'7r+ ± 3 

D+^a;7r+ ± 1 

D+ ^ p^K+ ± 1 

Dt^f27r+ ± 1 

L>,+ ^p+77 ± 3 

Dt^pW ± 2 

1:)+ ^ ± 2 

± 22 

D+^F*V ± 2 

Total ± 33 



the uncertainties of the background function fit to the MC, and the BTZ of 
the backgrounds. We study the effect of the first two sources in the following 
way: We vary each constant parameter including the normalizations and 
the shape parameters in Table 5.7 ± 1 cr, independently. The changes of 
yield {AN) are listed in Table 8.3. The normalizations of the backgrounds 
from the Cabibbo suppressed decay, Bg and A5 decays are independent 
from the shape parameters. But, several shape parameters for the same 
background are correlated, as shown in Tables 8.6-8.8. In order to take 
into account the correlation properly, the correlated shape parameters are 
grouped together. We calculate the product of the correlation coefficient 
matrix (M), with the row and column vectors of AN, to obtain a total 
systematic uncertainty. For instance, the systematic uncertainty from the 
Bs background shape parameters is: 



a% = ( AN^ AN,^ ANf, ATV^ ) M 

\ <'i / 

where M is a 4 x 4 correlation coefficient matrix returned from the fit to 
the Bg MC (see Table 8.6). The value of AN for each parameter is hsted 
in Table 8.3. 



AN^, 
ANf, 
ANc2 



1.2) 
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Table 8.2: B — > D*~^7r yield change due to the variation of the background BTZ. 



Mode 


Bn (%) 






D*+e- 




5.44 


± 0.23 


< 0.1 




D*^ 






5.44 


± 0.23 


< 0.1 




D*^ 






0.7 


± ? 


< 0.1 




D*^ 






1.30 


± 0.27 


< 0.1 


b"^ 


D*^ 


'P~ 




0.68 


± 0.09 


+0.1 
-0.2 


Total 


+0.1 
-0.2 



For the systematics associated with the BTZ, we vary the BTZ oi B — > 
D*~^n~ , B — > D'^p~ and the dominant modes in the remaining B — > D'^X 
backgrounds, ± 1 cr independently. We re-fit the background shapes using 
the MC, fix the shape parameters and re-fit the data. Table 8.4 lists the 
signal yield change due to the variation of the BTZ. Table 8.5 summarizes 
the signal yield change from the variation of the shape parameters and 
the BTZ. These changes arc added in quadrature to get the accumulative 
difference. The total change in the yield is 13 events, which modifies R by 
0.38. 

• Afo — > A+7r~: we follow the same scheme applied by Martin [101]. Using a 
generic S-decay MC, we first extract the top twenty largest contributing 
modes in the mass region 5.3 < Ma^^ < 6.0 GeV/c^, from each type of 
background: four-prong B meson, the remaining B meson decays, and the 
remaining A;, decays. Each dominant decay contributes iVbase events. Then, 
we generate a new distribution for each dominant mode, according to the 
shape determined from a large single-decay MC. The normalization of the 
new distribution is first Gaussian fiuctuatcd with a mean A^^ase' ^ sigma 
of A{BTZ)/{BTZ) and then Poisson fluctuated. For the measured decays, 
A{BTZ) is the uncertainty reported in the PDG. For the unmeasured B 
meson decays, A{BTZ) is assumed to be three times the uncertainty of the 
closest equivalent mode in the measured B meson decays. For the unmea- 
sured A5 decays, A{BTZ) is hypothesized to be of the BTZ. These 
Gaussian and Poisson fluctuated distributions are then re-combined with 
the other non-dominant modes. The combined background mass spectrum 
is refitted and the newly derived shape parameters are fixed in the fit to the 
data. The whole procedure is repeated 1000 times with different random 
seeds for the Gaussian and Poisson fluctuations. We plot the distribution 
of the A5 A+7r~ yield and record the RMS as the change in the yield due 
to the variation of the BTZ. Figure 8.1 shows an example of the A^ — > A+tt" 
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Table 8.3: B — > _D+7r yield change due to an independent variation of the fixed 
parameter value. 



Parameter 


AN 


f 

JDK 

A A/T 
I\MdK 


mass shift of S° ^ D+K' 


±4.97 
±3.82 


ctdk 


width of S° ^ 


±0.82 


f 

IBs 
l^Bs 


mean of Bg background 


±0.51 
-I- n no 


J 1 


fraction of the narrow Bg Gaussian 


±0.00 




width of the narrow Bg Gaussian 


±0.10 


(72/ (7l 


width ratio of the Bg Gaussians 


±0.04 


hb 

^^^b 


mean of A;, 


±0.22 
±0.22 


f^Ab 


width of Aft background 


±0.02 


^Ab 


lifetime of A5 background 


±0.21 


Tref 


lifetime of Dp background 


±2.62 


Cref 


width oi Dp background 


±0.27 


fH 


fraction of D*ti horns 


±6.38 




distance between two horns 


±1.04 


cth 


width of the horns 


±2.70 


/otherB 


fraction of the remaining B — > D^X 


±1.90 


Moff 


cut off for B — > D^X mass 


±1.02 
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Table SA: B — > D+tt yield change due to the variation of the background BTZ. 



Mode 



BTZ (%) 



b"^ 

b"^ 

B^ 
Total 







0.77 ± 


D*+7r 




0.276± 


D+e- 




2.14 ± 


D+p- 




2.14 ± 


D*+e- 




5.44 ± 


D*+p 


"^/^ 


5.44 ± 


D+n- 


7r° 


0.1 ± 


D*+7r 




0.7 ± 


D+aT 




0.60 ± 


D*+ai 


1.30 ± 


D*+p 




0.68 ± 



0.13 

0.021 

0.20 

0.20 

0.23 

0.23 

? 

? 

0.33 

0.27 
0.09 



+0.5 
-3.4 
+2.2 
-1.3 
+0.6 
-0.2 
+1.0 
-0.4 
+0.1 
-0.2 

± 0.2 
± 0.5 
± 1.0 
± 3.0 

± 0.1 
± 0.6 
± 4.5 



Tabic 8.5: Systematic uncertainty on the B D'^n yield from each indepen- 
dent parameter group. 



AA^ 




± 


5.0 


DK shape 


± 


3.9 




± 


0.5 


Bg shape 


± 


0.1 


hb 


± 


0.2 


Ab shape 


± 


0.4 


Dp + D*7r 


shape ± 


9.9 


/otherB 


± 


1.9 


Moff 


± 


1.0 


BTZ 


± 


4.5 


Total 


± 12.8 
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Table 8.6: Correlation coefficients returned from the fit to Bg MC. 





N 






fl 


(72/(Ji 


N 


1.000 












0.024 


1.000 








CTl 


-0.003 


-0.084 


1.000 






fl 


0.007 


-0.004 


0.847 


1.000 




(^2/(^1 


0.023 


0.133 


0.268 


0.647 


1.000 



Table 8.7: Correlation coefficients returned from the fit to MC. 





N 






^A, 


N 


1.000 










0.000 


1.000 








0.000 


-0.624 


1.000 






0.000 


0.699 


-0.508 


1.000 



Table 8.8: Correlation coefficients returned from the fit to D*7r and Dp MC. 





N 


Tref 






fH 


Sref 




^ref 


N 


1.000 
















T"rcf 


0.000 


1.000 
















0.000 


0.013 


1.000 












O'rcf 


0.000 


0.169 


-0.841 


1.000 










fH 


0.000 


0.549 


-0.688 


0.720 


1.000 








^rcf 


0.000 


0.294 


0.284 


-0.429 


-0.435 


1.000 






O-R 


0.000 


0.407 


-0.507 


0.507 


0.699 


-0.381 


1.000 






0.000 


0.029 


0.971 


-0.786 


-0.624 


0.261 


-0.455 


1.000 
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Table 8.9: A(, — >■ A^vr yield change due to the variation of the background BTZ. 



AN 


A, ^ A+X- 


+1.1 

-2.8 


four-prong B meson decays 


± 2.9 


remaining B meson decays 


± 0.9 


all the other A^ decays 


± 2.8 


Total 


+4.3 
-5.0 



yield distribution from the BTZ variation of the four-prong B meson back- 
ground. In addition, we vary the fraction of the Cabibbo suppressed mode, 
IacK, ^50°^ ^^'^ record the yield change. Table 8.9 summarizes the change 
of Afe 



A+TT yield. The accumulative an is 0.63. 



350 
300 
250 
200 
150 
100 
50 



MC (1000 iterations) 
4-prong B meson background 
RMS = 2.9 events 



^60^ iiO lis 166 1^5 1^0 1^5 200 



N. 



Figure 8.1: A5 A+vr yield from 1000 variations of the 4-prong B meson 
background BTZ. The RMS is recorded as the yield change. 



Measured Branching Fractions 

We use the BTZ from the world average in the PDG to estimate the physics 
backgrounds in our semileptonic signals as described in Section 7.2. We vary the 
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BTZ of these measured physics backgrounds by ± 1 a. We then calculate a^. Note 
that here the variation of the B{Aiy A+tt") does not include the uncertainty 
due to the measured A5 Px spectrum. The due to the measured BTZ is 0.43, 
0.75 and for ^ D*+,B^ ^ and ^ A+ modes, respectively. For 
the remainder of this section, we quote the systematic uncertainties in the same 
order. 

Unmeasured Branching Fractions 

The BTZ of several physics backgrounds in Section 7.2 have not yet been mea- 
sured, e.g.: A5 — > A+/°/x~z7^, or have just been measured by us for this analysis, 
e.g.:h.}y — > Ac(2593)+/x~z7^. For the first case, we use the estimated BTZ from the 
decay file of EvtGen, and our own derivation based on HQET. As we have no 
uncertainty input from the estimated BTZ, we assign a 5% uncertainty to the BTZ 
of the excited charm meson decays and a 100% uncertainty to the BTZ of the B 
hadron decays. Because the excited charm hadrons decay via strong interaction 
and conserve isospin symmetry, their BTZ could be inferred from Clebsch-Gordan 
Coefficients. While for the weak decays of B hadrons, allowable decay spectrum 
is wider. For the second case, we add (20 ® 20)% uncertainty in quadrature with 
the uncertainty from the preliminary measurement (see Table 7.5). The first 
20% arises from the unresolved disagreement of measured with that from the 
HQET prediction. The second 20% is due to the difference of the soft pion recon- 
struction efficiency between MC and data. We vary the BTZ by the uncertainties 
we assigned and calculate the shift of our measurement. The shift due to the 
unmeasured BTZ is 1.09, 0.91 and 0.50. 

Fcike ji estimate 

As noted in Section 7.3.1, the systematic uncertainties from the fake estimate 
originate from: 

1. The uncertainty from the fit to the charm mass spectra. 

2. The uncertainty on the probabilities for the pions, kaons and protons to 
fake muons. 

3. The uncertainty on the fraction of pion, kaon and proton in the hadron 
tracks. 

For each category, we vary the central value ± 1 o", independently. More de- 
tailed description can be found in Section 7.3.1. The resulting uncertainty on 
the amount of fake muons together with the central value are summarized in 
Table 8.10. We then vary the number of fake muons ± 1 cr and insert the new 
number into Equation 8.1. The total variation on R due to the fake estimate 
is 0.07, 0.07, 0.17. 
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Table 8.10: Summary of fake muon contamination. 



B D*^^fake 


B D^Hfake 




15 ± 3 


230 ± 19 


10 ± 9 



bb and cc background 

In Section 7.4, we notice a 10-40% discrepancy of the and inclusive 

b cross-section from PYTHIA with those from the data. In addition, we do not 
possess information about the relative bb and cc production rates between flavor 
creation, flavor excitation, and gluon splitting. Therefore, we assign a 100% 
uncertainty to the amount of bb and cc backgrounds. This changes R by 0.22, 
0.22, 0.04. 

MC sample size 

We have generated large MC samples for calculating the efficiencies of our signals 
and backgrounds, but there is a small statistical uncertainty due to the finite MC 
sample size. We use the uncertainties on the efficiencies to calculate aR. (Tr is 
0.28, 0.18, and 0.32. 

MC Pt{B) Spectrum 

We find discrepancies between data and MC in the Pt spectrum of and Af,, 
as described in Section 6.1. After reweighting the Pt spectrum of B^ and A^, 
we have observed good agreement of MC with the data as seen in Section 6.2. 
However, there is an uncertainty on the exponential slope of data/MC, pi in 
Figure 6.1, which is limited by the amount of data used for comparison with 
the MC. We vary pi ± 1 cr and re-weight the MC events after the analysis cuts 
numerically to calculate the efficiency change. In addition, we vary the variables 
which depend on the A^ Pt spectrum accordingly, eg: cross-section correction 
factors and the result of '^H^^T>^-'i)B(Ki,^K-^ ) , rggi ^ described in 

agol^PT>&.0)B{B^D+-K-) ' L J 

Section 7.1. The total variation on R due to the MC Pt spectrum of B hadrons 
is 0.38, 0.32, and 

' ' —0.50 

Pion Interaction with the Material 

One difference between our semileptonic and hadronic final states is the muon 
and the pion. The muon does not interact with the material via the hadronic 
(strong) interaction while the pion does. In order to model the track reconstruc- 
tion efficiency correctly, two things have to be right: 

1. The type and the amount of material in the detector. 
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2. The model that describes the hadronic interaction cross-section, the final 
state multiplicities and kinematics. 

We generate MC for the signals as described in Section 6.1 except that we switch 
off the hadronic interaction in the detector simulation. Wc compare the difference 
in the hadronic to semileptonic signal efficiency ratio, between the normal MC 
and the MC with the hadronic interaction off. This difference gives us an idea 
for the extreme situation, when the material is 100% wrong. For both the 
and Af, modes, the efficiency ratio changes by 4%. From the study of Korn [102], 
we know that the available CDF detector simulation underestimate the amount 
of material by 15%. In addition, a comparison between two programs which 
model the hadronic interactions, GHEISHA and FLUKA [103], has been done by 
Michael [104]. The FLUKA package is known to better reproduce the experimental 
data but currently it is not available in the CDF detector simulation. The effect 
of the hadronic interaction model estimated by Michael is 20%. Adding 15% and 
20% in quadrature, we get 25%. We multiply the 100% efficiency ratio difference 
described earlier, with 0.25, and get 1%. We apply a scaling factor, 1.01, to all the 
relative efficiencies, including the semileptonic background to hadronic efficiency 
ratios. We then re-calculate it! and find an is 0.22, 0.17 and 0.22. 

CMU reconstruction efficiency scaling factor 

The scaling factor to correct the difference of CMU muon reconstruction efficiency 
between MC and data has an uncertainty as described in Section 6.3. We vary 
the scaling factor ± 1 cr and calculate aR= 0.07, 0.05, and 0.07. 

XFT efficiency scaling factor 

We apply the XFT efficiency scaling factors data/MC in bins of inverse Pt from 
Herndon [80] to correct the signal and background efficiencies. The uncertainty on 
the kaon and pion scaling factors are varied ± 1 cr to evaluate a^. For the proton 
scaling factor, due to the limited statistics, we evaluate the systematic uncertainty 
following the suggestion in Herndon's analysis: we compare the difference by 
applying a constant efficiency scaling factor as shown in Figure 8.2, instead of the 
one based on the third order polynomial in Figure 6.15. For all the three modes, 
the systematic uncertainties are neghgible as expected, since the final states of 
our semileptonic and hadronic modes are almost identical and the difference in 
the ionization of the pion and muon is insignificant, is less than 0.01. 

Aft and Ac polarizations 

There is not yet a precise measurement of the production polarizations of Aj, and 
Ac, while the Standard Model predicts both particles are produced polarized. 
The angular distribution of of the A^ daughters is parameterized by 

dN 

-cxl + PscosO, (8.3) 

a cos 9 
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Figure 8.2: The relative proton XFT efficiency between MC and data in bins of 
1/Pt fit to a constant by Herndon [80] . 




Figure 8.3: Angle definition for the A5 production polarization, where the dashed 
line indicates the momentum of A^, in the rest frame of A^, and n is the polariza- 
tion axis normal to the beam proton-A^ production plane. 
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where Vb is the product of the polarization and the asymmetry parameter 
of the weak decay. is defined as the angle between the A+ momentum in the 
Af, rest frame and the axis normal to the beam proton-A^ production plane, n. 
Therefore, 

cose = P(Ac) -n, (8.4) 

where 

.^P(rtxP(A^ (8.5) 
|P(p)xP(At)|' 

Here "x" ("•") means vector (scalar) product of two vectors. See Figure 8.3 for 
the definition of and n. The angular distribution of Ac daughters is parame- 
terized in a similar way; 

dN 

-ocl + T^ccos^, (8.6) 

a cos 9 

where 9 is defined as the angle between the proton momentum in the rest frame 
of Ac and the Ac momentum in the lab frame, i.e. 

COS0 = P(Ac) -p, (8.7) 

Note that here P(Ac) is in a different frame from that in Equation 8.4. The value 
of Vj3 {Vc) is ± 1 for the polarized state and for the unpolarized state. 

The Bgenerator and EvtGen do not include the polarization of A^ and Ac. We 
use the default settings for the central value of R. We study the systematics due 
to a non-null polarization using the generator-level signal MC without the detec- 
tor and trigger simulation. We use the "acceptance-rejection (Von Neumann)" 
method [74] and reweight the MC according to: 

(l + 7'scos0) • (l + T'ccos^), 

where all combinations of Vb and Vc for values at -1, 0, 1 are used. Each MC 
starts from a different random seed. We apply generator-level analysis-like cuts 
and obtain the efficiency ratio for each combination. We compare these efficiency 
ratios with that from the MC generated with zero Vb and Vc- We find the 
efficiency ratios are mainly determined by the Vc as seen in Table 8.11, i.e. the 
efficiency ratios with the same Vc, but different Vb are consistent with each other. 
Therefore, we apply scaling factors from the two Vc values: -1 and 1 (± 1.017%) 
on all the relative efficiencies and re-calculate R. We find au — 0.37. 



Ac Dalitz structure 

The Ac from our A;, A+tt" and B —>■ A+/i~X signal, decays into p, K, and 
TT in the final state. However, any two Ac daughters could form an intermediate 

resonant state, see Table 8.12. The resonant structure is called the "Dalitz" 
structure in the literature, and is usually displayed with a Dalitz plot [105], where 
the invariant mass square of one pair of daughters is plotted versus another pair in 
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Table 8.11: e(A(, — >• A+tt )/e{A.f, —>■ A+/i u^) from each combination of Vb and 
Vc. 





Vc 




r(A,,^A+T 




s 


caliiig 


Factor 






~u ) 








3 


225 ±0 


010 








1 
1 





-1 


3 


193 ±0 


007 





990 ± 





004 





1 


3 


281 ±0 


006 


1 


017 ± 





004 


1 





3 


232 ±0 


006 


1 


002 ± 





004 


1 


-1 


3 


193 ±0 


009 





990 ± 





004 


1 


1 


3 


275 ±0 


009 


1 


016 ± 





004 


-1 





3 


235 ±0 


006 


1 


003 ± 





004 


-1 


1 


3 


274 ±0 


009 


1 


015 ± 





004 


-1 


-1 


3 


175 ±0 


Oil 





985 ± 





005 



the two-dimension. Figure 8.4 (left) shows the Dahtz plot from the B A'^fi^X 
data after sideband subtraction. If a resonance exists, a concentrated area near 
the mass of the resonant particle will be visible. The momenta of p, K and tt 
are affected by the Dahtz structure and A;, decays have different efficiencies for 
various structures. However, EvtGen does not take into account the interference 
of each resonant state. Each state is considered as an independent decay with a 
BTZ measured by E791 [106] and listed in Table 8.12. See Figure 8.4 (right) for 
the Ac Dalitz structure in the MC. 

Without a better model to describe the Ac Dalitz structure, we study the 
change in the efficiency ratio of A5 A+vr" to A5 ^ A+/x~I7^ by varying the BTZ 
in Table 8.12 ± la. We generate four sets of A5 — > A+7r~ and Af, A+/i~z7^ MC 
samples without detector and trigger simulation, where Ac decay is forced to one 
single mode. The efficiency of A(, decay with A+ — > pii'~7r"'" (total), is then the 
sum of the BTZ weighted efficiency of each individual Ac mode. 

R - ^c(A.^A+vr-) 

'''^ - s^iA,^ At f^-v,y ^^-^^ 

where Sc {Rc) is the total weighted efficiency (ratio) using the central value of 
each Ac BTZ. We re-calculate the absolute and relative efficiency, by varying BTZ 
of each Ac decay ± 1 cr; 

P ^ E!Bn,-e, + {Bn,+aj)-ej 
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Table 8.12: Ac decays with p, K, n in the final state. 



Decay Mode 


Bn{%) 


e(Ai,^A+7r ) 
e(Ai,^A+^-!7^) 


Rj 


p/s:*(890)° 


1.6 ± 0.5 


3.17± 0.01 


3.22 


A(1232)++X- 


0.86± 0.30 


3.24± 0.01 


3.23 


A(1520)7r+ 


0.59± 0.21 


3.30± 0.01 


3.23 


non-resonant 


2.8 ± 0.8 


3.24± 0.01 


3.23 


Rc 


3.23 ± 0.01 



_ £j(A„^A+;r , 



where Sj (Rj) is the total weighted efficiency (ratio) with BTZ of j*'* mode varied 
by ± 1 (T and other BTZ fixed. 

We find a fractional change of 0.3% after adding the difference of each Rj from 
Rc in quadrature. We apply this fractional change to the relative efficiencies of 
all the semileptonic backgrounds to the hadronic signal and calculate aR= 0.07. 



Afe Lifetime 

The world average A^ lifetime is lower than the theoretical prediction. A smaller 
Ab lifetime gives a smaller efficiency for reconstructing A;, decays. While we cut 
on the CT of A^ — > A+tt", the A;, — > A+/i~I7^ is not fully reconstructed and we 
actually cut on the pseudo-cr of the inclusive semileptonic decays. Therefore, 
systematics due to the uncertainty on the A;, lifetime may not cancel in our 
measurement. We study this effect by generating A;, — A+/i^l7^ and A;, — >• A+tt" 
MC without detector and trigger simulation. We vary the lifetime of A;, ± 15% 
around the central value: 1.229 ps. We compare the difference of the relative 
efficiency ratio from the central value. We then apply a scaling factor from the 
signals, on the efficiency ratios of the semileptonic backgrounds to the hadronic 
mode, and calculate aR=0.22. 



Semileptonic Af, decay model 

In Section 6.3, we introduce a scaling factor, fc, which accommodates the accep- 
tance difference between the flat phase space MC and the form factor weighted 
MC. We vary the /c ± 1 cr according to its statistical uncertainty and obtain 
an = ±0.57. 
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Figure 8.4: Dalitz structure in the sideband subtracted B A+yU~X data 
(left) and MC (right). The concentrated areas in the top figure indicate the 
existence of i^*(892)'' and A(1520). Clearly, the destructive interference between 
the resonant states are not simulated in the MC. 



8.1.2 Systematic Uncertainty for Each Mode 

Tables 8.13-8.15 list the result of systematic uncertainties as discussed above. 
The systematics from the external information are separated from the ones from 
the CDF MC and measurements. Table 8.16 summarizes the uncertainties from 
each category. The statistical uncertainties on the relative branching fractions 
are also listed for comparison. 



8.1.3 Consistency Check of R 

In order to detect any unexpected systematics in i?, we separate the data and 
MC into several groups of independent subsets according to the run number, 
vertex position, cr and Pt of the charm and B hadrons, and etc. We cross-check 
the consistency of the R within each group. Figure 8.5 displays the result of the 
cross-check, where the uncertainties in the figure are statistical only. The R from 
all the subsets are consistent with the other subsets in the same group. 

8.2 Measurement Result 

For the control modes, we measure the relative branching fractions to be: 



B*+7i- 



17.7 ± 2.3 (stat) ± 0.6 (syst) ± 0.4 (BR) ± 1.1 {UBR), 
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Table 8.13: Statistical and systematic uncertainties of '^^^—^^ ^'^^ . 



Source 




ctr 


Statistical 




± 2.3 


Measured SR, 


B" D*+TT- 




± 0.29 






± 0.31 






< 0.01 


± 0.43 


Unmeasured BTZ 


D*+Ti- 




± 0.05 


Df D*+TT- 




± 0.04 



Dt- 



± 0.03 
± 0.02 



B- - 




± 0.70 


B- - 




± 0.39 


^- 


-> D*+T-v^ 


± 0.31 


^- 




± 0.53 


^- 




± 0.34 


^- 




± 0.19 



± 1.09 



CDF Internal Systematics 
Fitting of B" D*+tt- < 0.01 

Fake estimate ±0.07 
bb and cc background ± 0.22 

MC sample size ± 0.28 

MC Pt{B°) ± 0.38 

TT interaction with the material ±0.22 
CMU reconstruction efficiency scahng factor ±0.07 
XFT efficieiic}' scaling factor < 0.01 

± 0.58 
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Table 8.14: Statistical and systematic uncertainties of 



Source 




ctr 


Statistical 




± 1.0 


Measured BTZ 






± 0.70 






± 0.22 


B- ^ D^ii-V^ 




± 0.08 






± 0.11 


T jlV^Vr 




< 0.01 


fs/fd 




± 0.01 


± 0.75 


Unmeasured BTZ 



D*+7r- 
■ D*+7r- 



± 0.01 
± 0.01 



B- - 




± 0.17 


B- - 




± 0.79 






± 0.39 






± 0.10 






± 0.09 



± 0.91 



CDF Internal Systematics 



Fitting of B D+n- 

Fitting of 5° ^ D+n'V^ 

Fake ^ estimate 

bb and cc background 

MC sample size 

MC Pt(5°) 

TT interaction with the material 

CMU reconstruction efficiency scaling factor 

XFT efficiency scaling factor 



± 0.38 
± 0.13 
± 0.07 
± 0.22 
± 0.18 
± 0.32 
± 0.17 
± 0.05 
< 0.01 
± 0.62 
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Table 8.15: Statistical and systematic uncertainties of 



g(Aj,^A+/j-i/;^) 
S(Ab-»A+7r-) ■ 



Source 




ctr 


Statistical 




± 3.0 


Measured BTZ 






+0.73 
-2.07 

< 0.01 


Total 




+0.73 
-2.07 


Unmeasured BTZ 



Aft 
Aft 
Aft 
Aft 
Aft 
Aft 
B- 



Ae(2593)+/x-l7^ 
Ae(2625)+//-I7^ 
E07r+//-I7^, Aft ^ 

A+7r°7rV"i^/x' ^6 

A+T-I7^ 

* A+p/x-i7^ 



E+ttV J^/i, Aft 



A+TT+TT I/^ 



± 0.21 
± 0.27 
± 0.24 
± 0.05 
± 0.20 
± 0.10 
± 0.11 
± 0.11 



Total 



± 0.50 



CDF Internal Systematics 



Fitting of Aft A+tt" 
Fake /j, estimate 
bb and cc background 
MC sample size 
MC PT{Ab) 

TT interaction with the material 
CMU reconstruction efficiency scaling factor 
XFT efficiency scaling factor 
Aft and A^ polarizations 
AjJ" Dalitz structure 
Aft lifetime 

Semileptonic Aft decay model 

+1.09 

-1.15 



± 0.63 
± 0.17 

± 0.04 
± 0.32 

+0.28 
-0.50 

± 0.22 
± 0.07 
< 0.01 
± 0.37 
± 0.07 
± 0.22 
± 0.57 
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which are consistent with the ratios obtained by the PDG, 19.7±1.7 and 7.8±1.0 
at the 0.7 and 1.1 a level. Finally, we measure the relative branching fraction 
to be: 

^^t^\ ^ ^l^^'^ = 20-0 ± 3.0 {stat) ± 1.2 {syst)^^^] {BR) ± 0.5 {UBR). 
o{J\h Ajn ) —I. I 

The uncertainties of the relative branching fractions are from statistics, CDF in- 
ternal systematics, external measured branching ratios and unmeasured branch- 
ing ratios, respectively. 

8.3 Estimate of the B{h ^ A+/i-F^) 

Wc have just presented the first measurement of the ratio of A^ exclusive scmilep- 
tonic to hadronic branching fractions. The ratio provides important input for the 
absolute branching fraction of A;, — >• A+/i~I7^ or A5 A+vr^. Leibovich et al. [33] 
predict B{Ab ^ A+tt') = 0.45% and ^(Aft A+fx-V^) = 6.6%, which gives a 
relative branching fraction of 14.7. However, the largest theoretical uncertainty 
from the functional form of the Isgiu-Wisc function is 30%, due to the assump- 
tion of the large Nc limit (see Section 2.3). Our measurement of the ratio has a 
19% uncertainty and may stimulate additional theoretical work. Multiplying our 
Aft relative branching fraction, with our derivation of B{Ah —>■ A+tt") from the 
CDF measurement of -Hi^TX^-^^^iM^^t^-) rgg. g^^^.^^ ^ ^. 

<Tgo{PT>6.0)B{B"^D+Tr-) ^ ^ 

B{Ab A+TT-) = (OAI ± 0.19 {stat © syst) ^JJ'JJg (Pr)) %, 

we obtain 

B{Ab ^ A+i^-V^) = (s.l ± 1.2 {stat) ^ {syst) ± 4.3 {B{A, ^ A+tt"))) %. 

which is also consistent with a recent DELPHI result derived from the A^ — > 
A+/i~I7^ form factor measurement [37], 

B{A, A+/.-I7,)™^ = {^■0^ll{stat)^];l{syst) 

Combining our and DELPHI'S numbers, we obtain {5.5 ± l.S{statQ) syst)) %. 
Our relative branching ratios and the derived B{Ab — > A+tt"), B{Af, — > A+//~I7^) 
are all in agreement with the predictions by Leibovich et al, within large uncer- 
tainties. Note that the dominant uncertainties of B{Ai, A+tt") arise from 
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and B(A+ — > pK 7r+). New CDF-II measurements of ^ are anticipated. How- 
ever, a better measurement of i3(A^ — > pK~7r~^) has only been proposed by Duni- 
etz [107] and Migliozzi [108]. Improvements in the B{Ai) A^vr") will reduce the 
uncertainties in our determination of the exclusive semileptonic branching ratio. 



8.4 Conclusion 

We analyze 171.5 pb~^ of data collected with the CDF-II detector in the pp 
collisions at ^/s — 1.96 TeV. Using a novel secondary vertex track trigger, we re- 
construct 1237 ±97 B A'^fi~X decays and 179 ±19 A5 ^ A+vr^ decays. This 
is the largest A^ sample in the world, which enables us to measure the relative A^ 
branching fractions and examine the Heavy Quark Effective Theory. We have also 
observed several A^ semileptonic decays which have never been seen in the other 
experiments: A^ Ac{2593)+ pL' X , Ai, Ac(2625)+/i-X, A^ E°7r+/i-X, and 
Aft — > E^+7r~//~X. In addition, we reconstruct the D*^ and 5^ -— > 

decays similar to our A;, decays and use them as the control samples to under- 
stand the issues associated with the A;, measurement. After the estimate and 
the subtraction of the background in the inclusive semileptonic signal, we correct 
the yield observed in the data with the trigger and reconstruction efficiencies 
obtained from the Monte Carlo. We find the relative branching fraction of the 
control modes in good agreement with the values obtained by the PDG. We 
measure the ratio of A^ branching fraction to be: 

^^^A ^ ^1+^^^ - 20.0 ± ?>.Q{stat) ± l.2{syst) ^^^]{BR) ± 0.5(C/Si?). 

The uncertainty is dominated by the size of the data sample and the branching 
ratio of A;, — > A+tt". More data and a more precise measurement of B{Ai, 
A+7r~) in the future will immediately improve our relative branching fraction 
measurement and our determination of B{Af, — > A+/i~I7^). 



8.5 Future 

In addition to the measurement presented in this dissertation, several measure- 
ments can be performed at CDF and increase our knowledge of the A^. 



Aft lifetime Previous analyses [7,109-111] used semileptonic decays to measure 
the Af, lifetime so the A;, was not fully reconstructed. A scaling factor based on 
the MC was applied to convert the total momenta of the observed daughters to 
the Afe parent momentum. On top of that, possible backgrounds from the other A;, 
decays might not have been included. Therefore, a A^ lifetime measurement using 
fully reconstructed decays, such as A^ — > A+tt" or Aj, — > J/^A, can provide a 
model-independent comparison with the current world average and the theoretical 
prediction and possibly resolve the discrepancy as noted in Chapter 1. 
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Aft production cross-section Instead of measuring the ratio of two branching 

fractions, we can measure an absohite branching ratio directly if the A5 production 
cross-section at the Tevatron is known. An ahernative way is to measure the ratio 
of the A;, to production cross-section. As noted in Section 8.3, an accurate 
measurement of the will reduce the uncertainty on the B{Af, A+tt") when 
normalizing the A;, decays to the B meson decays. 

Af) A+/i~I/^ form factor and \Vcb\ As the three-dimensional (3D) vertex 
reconstruction using the SVX information is not fully developed at CDF-II, yet, 
this analysis made requirements only on the two-dimensional vertex (L^y) when 
reconstructing Ag/x events. Once the 3D vertex reconstruction is mature, we can 
use kinematic constraints to obtain the neutrino momentum and calculate the 
scalar product of the A5 and the A^ four- velocities, w. The w distribution can be 
fit to the Isgur-Wise function to obtain the slope parameter, p^, in Equation 2.13. 
When an exact form of the Isgur-Wise function is established, the A;, semileptonic 
decays can give a \Vcb\ competitive with that from the B meson decays [112,113]. 

CP asymmetry in the A^ — > pn and A^ pK decays Bensalem, et al. [114] 

and Dunietz [115] predict a large T violation (or CP asymmetry) in the charmless 
Afe decays in the Standard Model. The asymmetry can be studied without a need 
for the flavor tagging since the baryon decays are not affected by the mixing. 
However, A;, charmless decays, such as A(, pir and A;, —>■ pK, have not been 
observed, yet. A recent search for these decays at CDF-II by Carosi, et al. [116] 
set an upper limit on the branching fraction B{Ab ~^ ph) < 23 ■ 10~^ at 90% 
confidence level, and improved the previous upper limit by ALEPH [117]. More 
data in the future can either improve the upper limit or result a first observation 
of the signal. Once enough signal is established, a CP asymmetry measurement 
can be made. 
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